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Abstract A study has been made on the influence changes in surface tensio 
of interfacial area in different types of distillation equiy it. It has been 
films are stabilized when the surface tension of the reth ere down tl 
with decreasing surface tension liquid films break up into narrow rivulets o 
area in spravs, however, is not much influenced by changes in surface tensi 
kor agiven liquid mixture ibuliza f filtos can be demonstrated t ersing 
the direction of mass transfer magnitud the efleets depends distinetly on tl 
of the surface tension gradient developing in the 
In apparatus where the interfacial area ts mai f the tilm tvim 
mixtures causing the surface tension of the retlus t increase may | ty 
as those in svsterms in which the surface tension of the reflux decreases. In the 


equipment with interfacial area mainly in the form of droy 


plet sprays, this influence 


on mass transfer is but small 


Résumé — Les auteurs ont ctudic Teffet de changements dans la tension 
formation de Vinterface dans differents types de colonnes de distillation. Ils o1 
films de liquide sonte stabilisés quand la tension superticielle du liquide ce 
tandis qu'une diminution de la tension superticielle cause la rupture de ces filn 
petits ruisseaux ou gouttes. Cependant, dans des emes de gouttes linterfac 
influencée beaucoup par des changements dans la tension supertictell 

Pour un mélange de liquides donne, on peut demontrer la rupture ou la stabilisation des 
films en changeant la direction du transfert de tmiasse la grandeur des effets depend nettement 
de la valeur du gradient de la tension superticielle qui se forme dans le liquide en reflux 

Dans les colonnes ot Tinterface a principalement la forme d'un film, les vitesses du transfert 
de masse, pour le cas de melanges qui provoquent une iwmentation de la tems superticiell 
du liquide de reflux, peuvent étre deux fois u plus) plus clevees que celles qui presentent 
lorsque la tension superticielle du liquide de reflux diminue Dans le cas d'installations indus 
trielles comportant des interfaces principalement en forme de gouttes, cette influence du systéme 


sur le transfert de masse nest que faible. 


Zusammenfassung Kine Untersuchung uber den Einfluss von Anderungen der Oberflachen 
spannung auf clic Grenzflichenbildung in Destillationskolonnen verschiedener Art hat ergeben 
dass Fliissigkeitsfilme stabilisiert werden, wenn die Oberflachenspannuny les Riickflu 

zunimmt, wihrend bei abnehmender Oberflachenspannung soleche Filme zu schmalen StrOmen 


oder Tropfen zerrissen werden. Jedoch wird die Grenzfliche in Tropfensystemen von Ander- 


ungen der Obertlichenspannung nicht stark beeintlusst 


*The contents of this paper were read at an “ Interne Arbeitssitzung des Fachausschusses Destillation, Rektifi 


kation und Extraktion der V.D.1. Fachgruppe Verfahrenstechnik,” 24 April, 1956, Bingen-on-Rhine (Germany) [13]. 
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Fiir cin bestimmtes Flissigkeitsgemisch lasst sich das Zerreissen oder Stabilisieren der Filme 


HARMENS 


nachweisen, indem die Richtung des Stoffaustausches umgekehrt wird: die Effekte werden 


deutlich von der Grosse des sich im Rickfluss entwickelnden Oberflachenspannungsgefilles 


bedingt 


In Apparaturen, wo die Grenzfliche hauptsiichlich in Form eines Films auftritt, kann bei 


Gemischen, die die Obertlichenspannung des Riickflusses steigern der Stoffaustausch zweimal 


so schnell (oder schneller) sein als der bei Systemen, wo die Oberflichenspannung des Rick flusses 


abnimmt. Beitechnischen \pparaturen mit einer Grenzflache hauptsiachlich in Form eines Sprudel- 


lL INTRODUCTION 


IN operations of mass transfer between a gas 
or vapour and a liquid the transfer rates obtained 
ae pn nd on the rates of diffusion in the two phrase s 
and the magnitude of the contact area between 
the phases. These quantities are influenced by the 
type of equipment in which the operation ts 
carried out, the physical properties of the phrase s 
and the operating « onditions appli d. Fre que ntly 
the effect of these factors is mutually de pend nt 
and cannot clearly be se parat d. As an xample, 
the change in op rating conditions n t\ affect the 
mass transfer rates on a bubble cap tray quite 
differently from those in a packed column. This 
is particularly true in connexion with the subject 
of the present study, which discusses the influence 
of surface forces on vapour-liquid interfacial 
area and the resulting mass transfer rates. 

In general, little attention has been paid in the 
past to the formation of interfacial area in actual 
vapour-liquid mass transfer operations. Most 
of the work done concerns the study of bubble 
sizes emerging from submerged orifices, the degree 
of wetting of packing units In a pae ked column 
or the characteristics of liquid sprays. The 
formation of interfacial area under conditions of 
mass and or heat transfer is not covered by thes 
previous studies. It is found that the interfacial 
area depends primarily on the flow conditions of 
the phases and mass forces. Surface forces 
influence the interfacial area to a minor degree, 
the general effect being a smaller area with 
mcreasing surface tension. 

This situation, however, may change consider- 
ably if the formation of interfacial area takes 
place under mass and/or heat transfer, particu- 
larly if owing to these phenomena surface tension 
gradients develop along the vapour—liquid inter- 


face. It has long been known that surface tension 


betts ist dieser Eintluss des Svstems auf den Stoffaustausch nur gering. 


gradients cause rapid movements in the surface 
and may lead to a spreading or a contraction of 
the surface [7]. Such contraction phenomena 
have indeed been observed under certaim con- 
ditions in the absorption of ethyl aleohol vapour 
[8], and ammonia [1] into water films, It is 
readily understandable that the contraction of 
the liquid surface into small rivulets will have a 
profound effect on the overall mass transfer 
rates. Presumably many hitherto unexplained 
anomalies in mass transfer rates reported in 
literature 3, 10) have to be ascribed to this 
type of surface effects. 

It is noteworthy that the relation between 
mass transfer and interfacial area has recently 
received attention in the case of liquid liquid 
contacting (9). The stability of dispersed droplets 
has been shown to depend on changes of surface 
tension caused by the transfer of a solute through 
the interface. 

Complex surface phenomena as quoted above 
are in general very difficult to analyse theoreti- 
cally. The accumulation of a general insight into 
these problems is a first objective that is best 
obtainable by carrying out adequate experiments. 
It has been found that in the case of vapour 
liquid contacting (distillation) the latter can be 


realized relatively easily, The first purpose of 


this paper is to describe such experiments, and 
to give a discussion and possible explanation 


of the results. 


2, Surrace TENSION CHANGES IN 
DISTILLATION 


The surface tension or interfacial tension of the 
reflux is subject to changes during the downward 
flow of the liquid in the fractionating column. 
This ts caused by the changes in composition and 


the increasing temperature, It is possible that 
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temperature and composition effects reinforce or 
neutralize each other in this respect. Since in 
vapour-liquid systems temperature and composi- 
tion are closely interrelated, the possible variations 
in surface tension can easily be calculated from 
the equilibrium compositions and temperatures 
and the surface tension of the pure components 
at the same temperatures, 

It is then found that for systems composed of 
compounds of homologous series, the possible 
changes in surface tension are about 2-3 dyn em. 
For mixtures of compounds of a less related nature 
much larger variations in surface tensions are, 
however, possible. Mixtures of hydrocarbons and 
alcohols, for instance, may show surface tension 
ranges of about 5dyn em at their boiling points. 
Pure hydrocarbon systems can also show large 
surface tension gradients, if they are composed 


For 


surface tensions in the svstem benzemn 


of paraflins and aromatics. instance, the 


plan 


range between 21 and I2dyvn em, those of the 


mixture heptane toluene between 12 and 


From these data it is seen that in 


dyn em. 
a normal distillation the reflux running down in a 
fractionating column decreases in surface tension 
whereas with the second mixture 


It will be 


tension only deve lop 


in the first case, 
an merease occurs. clear however, that 
these changes in surface 
if the operation of the column is such that a 
gradient in concentration and temperature occurs 
along the column. In the case of operation close 
to minimum reflux the changes in composition, 
and therefore also the gradient in surface tension 
in the reflux stream are small, even if there is a 
wide gap between the surface tensions of the pure 
components of the system. This applies in a 
to all mixtures where relative 
that total 


reflux, the concentration gradient in the reflux 


more general sense 
volatility is small. In case, even at 
\n example of such a system 
is the widely test 
methyleyclohexane. Though the surface tensions 
of the components (12 and 15 dyn/cm) differ 
considerably, the (1-07) 
allow of the building up of noticeable surface 
the normal 
conditions of distillation. 

From the above description it is seen that three 


may be negligible. 


used mixture n-heptane 


low x value does not 


tension gradients in reflux under 


types of systems can be distinguished with respect 
to the changes in surface tension de veloping in the 
reflux flow. For convenience these types will be 
denoted as negative, positive and neutral re- 
the 


in positive systems 


spectively. In negative systems, reflux 


decreases in surface tension ; 
it increases. In both cases the relative volatilities 
in the system have to be large enough to allow 
the development of appreciable gradients. The 
term neutral is given to those systems In which 


either the components do not show a difference 


in surface tension or in which the relative vola- 


tility is very low and the gradients in surface 


tension are consequently always small. 


The mixtures employed in the present study 


are listed in Table 1. 


Table 1. 


Survey of binary systems used 


Surface tension 


System Boiling point | at boiling point 
(©) (dyn/cm) 
n-heptane 12 
methvleyclohexane 101-8 15 
n-heptane 12 
toluene 110-7 18-5 
benzene 80-2 21 
n-heptane 12 
benzene SO-2 21 
toluene 110-7 
benzene 21 
cyclohexane SO-8 17-5 
(azeotrope) 77°5 20 
ethanol 78-3 18 
2.2.4-trimethylpentane 99-1 11-5 
(azeotrope) 71 16 
Most of the experiments were done with 


the negative benzene-n-heptane, the positive 


and n-heptane—methyleyelo- 
Because of the low 


taken as a 


n-heptane toluene 
hexane mixtures. x value, 
the latter mixture will, however, be 
representative of a neutral system. 

Two systems are of particular interest—viz. 
the mixtures ethanol-2.2.4—trimethylpentane and 
benzene—cyclohexane. These mixtures both show 
an azeotrope at about equimolar compositions. 
As a result 
different in sign, according to whether the mixture 


the changes in surface tension are 
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distilled has a low or a high content of the low- columns in which the liqguidl ts dispersed by a 


boiling constituent. mechanical atomizing device. 
It is found that, with respect to the effeet of 


changes in surface tension. supported interfacial 


3 Us) 


area behaves quite differs ntly from unsupported 


The formation of interfacial differs 


interfacial area. Therefore. representative types 


according to tha ol equipment used, 


ol equipment of both yroups detined above were 


Roughly, two vroups by distin uished, In 


used in the present study. Particulars about the 


the first. interfacial area is obtained by spreading 


equipment are listed in Table 2. 


the liquid itt thin films over area ilt wt\ present 
in the equipment. In this case therefore, one can Surrace Tension wirn 
Sy ik of supported mitertac ial area. \ prototype SUPPORTED int) RFACIAL Ans 
of thus erou thy wetted wall ceolun 1) Packed 
(1) Wetted wall tyjpn columns 


columns and moditied wetted wall co umns such 


as Viereux columns also belong to this group. \ most striking example of the influence of 


surface phenomena on the performance of dis 


In the second tyyn ol equipment the inter 


facial area is obtained by mixing th quid and tillation equipment is obtained when operating 9 
vapour phases. As a result, the formation of @ Wetted wall column, HETP data measured for 1952/ 
vapour bubbles and droplets occurs; the inter the small concentric tube column at total reflux 

facial area in this case is not direct supported are given mm kig l In the case of the n hept ne 

by the equipment | xamples typ of methvleyclohexane system, the liquid phrase 

tray columns which may operat ther witl 

predominantly liquid’ interfacia rea (foam 7 

and bubble systems) or droplet interfacial area af 


Pure droplet interfacial area ts obtained in spray 


Tabl 2. Surve equipmcu fine 


(1) Concentric tube colwme 


Length 200 mm, internal diameter 6 mm. width of P ya 


annular space 0-60 mm. Only the outer wall wetted 


inner wall dry 


columns 


Varving length, diameter 25 mm, distance between L& 


similar indentations 40 mm 


Packed column 


(3) Length 250mm, diameter 40 mm. Packing : pores 4 


lain Raschig rings, 6 mm diameter 


Packed column 


Length 100 mm, diameter 25 mm. Packing : metal 


Fenske helices, 38 mum. or Dixon wauze rigs, mm 


(5) Spray column 


ric. 1. Separating power of concentric tube column 


IS rotating discs, 40mm diameter, spaced 30 mm 


with different test mixtures 


A n-heptane methyl yclohexam 
Oldershan Siete plate column 


apart on axis. Column diameter 100 mr 


venzene n-heptane 


Diameter 25 mm. plate spacing modified to LOO tin, 


2 or 5 plates. 


(7) Perforated plate column showed compl te wetting of the column wall and 
Diameter 450 mm 1 per cent rea 


the separating power observed agreed very well 


10mm hole size Plate spacing 400 nn 


with the theoretically predicted values [12]. 


é 
10) 
predic tice 


However, with the negative benzene n-heptane 
system very bad separation was obtained, the 
HETP being about ten times the 


value. It was observed that in this case the column 


theoretical 


wall was poorly wetted : the liquid flowed down 


in narrow rivulets. Improvement could be 
obtained by fixing a wire coil of low pitch (about 
2mm) along the column wall. 
capillary forces, the liquid was now retained and 
spread out on the wall and nearly the same 
HETP was observed for the benzene—» heptane 
mixture as for 
(Fig. 1). 


similar 


n-heptane methyleyclohexan 


difference in separating power 
caused by a different degree in wetting is observed 
with the 


were obtained for a variety of svstems under a 


Vigreux columns. ‘Total reflux data 


wide range of load conditions. In the range of 


Vapour velocities of about 0-4 to 0O-6m see the 
separating power was approximately constant. 
HETP values observed in this 


plotted itt 


The average 
Fig. 2 against average 


data for the negative mixture ar nearly twice 


ranve are 


mixture composition. As shown, the 


500; 

300} 
uy 200 
I 

100} 

Mole /o volatile component in. mixture 
Fic. 2 Separating power of 25 mm Vigreux column with 


various test mixtures. Vapour velocity 0-6 m see, 


benzene n-heptane 
v nzene cyclohexane 


n-heptane toluene 


n-heptane methyvleyclohexane 


those observed for the n heptane methyleyelo- 
hexane and the n-he ptane toluene systems. The 
wetting of the column wall observed visually was 
in accordance with these results : 
benzene n-heptane partial wetting and the for- 
rivulets, as illustrated in’ Fie. 3, 


mation of 


Because of 


in the case of 
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occurred, whereas the other two mixtures showed 


almost complete wetting. The similar behaviour 
of the latter two explains the almost identical 
HETP values for these mixtures. 

The observation that particularly the benzen 
n-heptane mixture gave poor wetting only scems 
explicable by surface effects. Since there is no 


major difference in the static ~ surface tensions 
of the different systems employed, it was believed 
that the sign of the changes in the reflux surface 
tension might perhaps be responsible. In order 
to test this supposition, some more experiments 
First. 


carried out with 


were done with the Vigreux column. 
total reflux distillations wer 
tha mixture benzene cyclohexane. \t high 
concentrations, this constituent is the 


because ot the 


ln 
less volatile azeotrop inthe 


system. Therefore, reflux surface tension in- 
creases downwards in the column (** positive 

mixture), Complete wetting was observed in this 
case. With a low be nzene content of the mixture, 
cyclohexane is the less volatile component and 
the surface tension of — the reflux decreases 


negative” As with benzene ptane 


svstem). 


channelling of the reflux was now observed. 


The separating power measured for the two 
mixture compositions was in accordance with the 
wetting characteristics found (see Fig. 2). 

In a second series of experiments it was attempt- 
ed to measure mass transfer rates with th 
normally ** positive systems n-heptane- toluene 
and ptane—me thy le yclohexane under nega- 
tive” conditions. This cannot be realized in 
normal distillation ; however, by countercurrent 
absorption ot vapour rich in volatile component 
into a licguicd rich in the heavy component the 
direction of mass transfer may be reversed and 
thus also the sign of the surface tension changes 
in the downtlowing liquid. Moreover, by changing 
the length of the Vigreux column or by adjusting 
the compositions of entering liquid and vapour, 
the surface tension gradient in the liquid can 
easily be varied and its influence determined. 
Data on the absorption experiments with 
Vigreux columns of various lengths are collected 
in Table 3. 


and composition, — the 


Besides top and bottom tlow rates 


average driving force, 


ic. the average difference between the operating 
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line and equilibrium line is given (in terms of 
liquid composition ). The latter is re presentative 
for the average changes in concentration (and 
hence for changes in surface tension) in the reflux 
At a high 
force, mass transfer is rapid and a large gradient 
the \t 
driving forces, the reverse is true. 

As all 


reflux surface tension decreased in the absorption 


over one theoretical tray. driving 


in reflux surface tension results. low 


expected, under conditions where 
experiments, channelling of the retlux occurred, 
This became more severe the higher the driving 
the HETP 
driving force. This is clearly 
illustrated in Fig. 4. in which the HETP data are 
plotted. It that both the 
data fall on one singk Furthermore, the 
HETP value obtained in normal distillation with 
the 


correlates 


force applied. Consequently, also 


increased with the 
for svstems 


is sech 


line. 


muxture 
the 


n-heptane 


if plotte d at 


bee nzcne 
this 


* channelling 


well with 
average driving force applied in the distillation 
experiments. 

Fig. 4 also shows the data obtained under fully 
The 
value for the n-heptane-methyleyclohexane sys 
the 


wetting conditions with various mixtures. 


tem at high driving force was obtained by 


Table 3. 


Absorption riments in 


Concentrations im 


mole 


HARMENS 


absorption of methyleyclohexane into n-heptane 
(see Table 3). It 
of positive changes in surface tension there ts no 
of the the HETP, 


presumably because is almost 


is seen that under conditions 


dri force on 
the 


may 


influence 
wetting 
complete. Further, it be concluded from 
Fig. 4 that at low driving forces, the HETP of 
negative systems seems to approach that of the 
positive mixtures, te. wetting becomes complete 


and channelling disappears, 


O*40 


mole fraction 


effect of 


column 


driving force on separating power of 
tested test 


Vapour velocity mi 


Kia. 4 


Vivgreux with various mixtures. 


ln nzem 


" hept 


n-heptane 
ptane thyleyclohexane 


toluene 


columns. livrage velocities 0-3 m/sec. 


pe 


vapour 


cent heptane 


Feed rates Product rates 


Column 


length 


(mmole hr) hr) 


System 


bottom bottom 


lop 
16 

15: 

16 

17 


17 


n-heptane toluene 


n-heptane 


methyleyclohexane 


Compositions, (mole per cent) Average 


Feeds Products Theoretical HETP 


plate 


ing 


force 


(mem) 


top hottom top bottom (mole 


28 


16: 
17 


oc 
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9 
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Driving force, 
— 
top 
1-48 17-8 23-2 340 
= 1-06 17-6 22-3 > 2-3 
16-6 20-5 0 21 . 
0-65 18-1 14-7 20-8 0 10) 4h 52:55 37 13 
’ 21:9 155 118 | 25-6 100 30 | 145 30 
1-48 20-0 16-6 12-1 25-4 0 100 il iv 3-8 ooo 
141 16-5 11-6 20-0 0 100 23-5 él 26-5 2-3 
0-88 21-1 169 146 234 0 100 23 57-5 27-5 21 420 
0-65 15-7 18-6 0 45 +h 1-3 
O48 16-3 16-5 14-1 0 100 55 os 
0-48 19-7 18-1 162 «216 100 0 71 38 33 20 


hia 3 Photograph of channelling phenomenon in 


Vigreux column when distilling benzene—n heptane 
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There is one remarkable feature in the correla- 
tion of the HETP data in Fig. 4 which should be 
The 


between n-heptane and methyleyclohexane is less 


mentioned. difference in surface tension 
than half the difference in surface tension be- 
This 


means that since for all mixtures the range of 


tween n-heptane and toluene or benzene. 


liquid concentration under channelling conditions 
is about the same, the surface tension gradient in 
the reflux during the n-heptane aromatics tests 
was about double the gradient in the n-heptane 

methyleyelohexane experiments. In spite of this 
fact the intensity of channelling in the formation 
of interfacial area seems to be the same. 

The observations of channelling liquid) phases 
reported in literature are also to be interpreted 
on the basis of decreasing surface tensions. 
and Parsty [8] found that the contrae- 
tion of water films in grid type tower packings 
occurred when ethanol vapour was absorbed. In 
this case surface tensions of the liquid film decrease 
because of the increasing ethanol concentration. 

Boxp and Dona. [1] observed the break-up 
of water films when hydrogen chloride or ammonia 
was absorbed in a rippling film in a wetted wall 
column. Owing to the heat of absorption, a 
decrease in surface tension also develops in these 
cases, It is interesting that the break-up did not 
rippling of the liquid film was 


occur when 


prevented, 


Packed columns 


(2) 

In packed columns the interfacial area is sup- 
ported by randomly distributed packing elements. 
In view of the results obtained with the wetted 
wall columns, the effect of the direction of surface 
tension changes on the HETP was also deter- 
mined for this type of columns. Two types of 
packing were investigated in relatively small 
columns : 6mm porcelain Raschig rings (in a 
40mm diameter column) and two fine metal 
In the 


case of the Raschig rings, efliciency varied but 


packings (in a 25 mm diameter column), 


little with the boil-up rate at vapour velocities 
of 0-2 to 0-3 msec. HETP data measured under 
these conditions are plotted in Fig. 5. It is seen 
that just as with Vigreux columns, the HETP 


for the n-heptane methyleyclohexane and n hep- 
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tane-toluene mixtures are about the same, 


whereas for the negative benzene n-heptane 


mixture about the double HETP values are 
recorded. 
300, 
> 
I 
4 6C BO 100 
Mole volatile nponent in mixture 


Fic. 5. Separating power of 6mm Raschig rings, with 


various test mixtures, 
benzene n-heptane 


n-heptane toluene 


Vapour velocity 0-2-0°3 m /see. 


n-heptane methylcyclohexane 


HETP 


observed with positive and negative systems was 


\ similar difference between values 
found for the fine metal packings such as Fenske 
helices and Dixon gauze rings. Table 4 shows the 
data obtained : 

Table 4. HETP data for fine metal packings. 
Average liquid composition 65 mole per cent volatile 


component 


HETP 


Vapour velocity 


Packing System (m /sec) (mm) 
Fenske helices n-heptane 1S 
toluene 
benzene 0-06 30 
n-heptane 
Dixon rings n-heptane 0-12 20 
toluene 
benzene 35 


n-heptane 


Though in the cases of the packed columns the 
actual flow of liquid over the packing elements 
was dillicult to observe, it may be safely assumed 
that the differences in separating power with the 
different systems must again be ascribed entirely 
to a different degree in wetting. 
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5. SURFACI 


I NSU 
(1) Spray colui is 


In order to study the effeet of chanees mn 


surtace tension when mitertacial are formed by 


dispersed drops use was made of a spray colunm 


in Which the droplets were cvoncrats 1) rotating 
alises, The column diameter wis 
provided with i central rotatin \ cnn Which 
dises of were | The 
retlux was euided onto the topa spraved 
it to the wall. From the wall, tl und was 
collected in a per pheral outter. fro whieh at 
was fed back tor the ext rotatin Ih total 
the cohimn contamed ral 
eutters tin power « 
at 1000 of t} Was cl ned wit! 


data obt titi 


hia. 6 Separating power f spray colun ith 
fest muxtures 
heptane toluene 
benvzen carne 


difference in mass transfer rate between the two 
svstems is observed result entirety 
unexpected, since im this” particu type of 
apparatus droplet formation occu! external 
forces MOreovel onee they irc ed ther 
seems litth possibility for interaction and co 
alescence between thy droplets by i ot the 
short residence time. Furthermore, the droplets 
formed are of such small size that the scarcely 
ln exper ted to disintegrate still fur by the 
action of decreasin surtace tensor Since only 
such disintegration or coalescence car ilter the 
magnitude of imterfacial area thy bsence of 


effects of changes in surface tension on sé parating 
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power for this particular type of equipment seems 


sullicrently « plain a. 


Perforate plate columus 


In perforated plate columns where th plates 
have a low percentage of open areca and small 
rloratrons. the lisp rsion of thy 


leads to the 


Vapour into 
the liquid frequently formation of 


bubbles or foam \ specimen of equipment of 
this is the Oldershaw perforated plate 
Plate 


column were determined with a variety of svstems. 


Data obtained im 


laboratory column elicreneres for this 


recorded in Table 5. part of 


Which is yl bie. 


It is now seen that approximately the san 


eroney (about 50.55 per ce nt)is obse rved 


for neutra or neovative mixtures, whereas 


by 7 Separating power of Oldershaw sieve plates 


with varhous test mixtures 


Ir mole of 


volatile « omponent 


Vv benzene toluene MO 
benzene n-heptane 
n-heptane methvileyclohexane 

pot carne toluete 15 
n-heptane toluene 75 


for the positive carne toluene system nearly 


the double eflicieney of 90 per cent is reached, 


Phis result is obvious when the behaviour of 


the licyuaicl op the trays is observed: with 


the w-heptane toluene system severe foaming 


whilst the other systems, in contrast, 
high bed. 
clearly ck monstrated by the photographs inh Fig. Ss. 


ck Ve lops, 


show a nuxaderately spray This Is 

\s with the Vigreux column the relation of the 
observed poe nomena with the direction of changes 
in surface tension was investigated more closely 


by the distillation of an azeotropie mixture and 
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Pic, 8. Tray action on Oldershaw sieve plates with 


benzene n-heptane (left) and n-heptane toluene system 


(right) respectively. Vapour velocity about 0-3 m see 


7 


Fic. 9% Tray action on Oldershaw sieve plates with 


ethanol 2.2.4 trimethylpentane mixtures, Left 7°, etha- 
75° 


nol, trav eflicieney 30 percent. Right ethanol 


tray efficiency about 70 per cent 


i 
~ 
VOL, 
by 9 
ie 
ye =. 
a 
4 


The influence of surface phenomena on the performance of distillation columns 


Table 


Efficie ney of Oldershax 


sieve plate 


Irerage mole 


System volatile comp 


n-heptane methyvleyclohexane 


benzene toluene 


benzene-n-heptane 


n-heptane-toluene 


ethanol-2.2.4 trimethylpentane 


Vapour velocity Tray efficiency 


(m/sec) { ) 


Foam Height 


(mm) 


The 
azeotropic mixture ethanol 2.2.4 trimethylpen- 
tane was used. 


by carrying out absorption experiments, 


At ethanol concentrations below 
53 mole per cent, this component is the most 
volatile. Since the 2.2.4 


trimethylpentane is lowest, the system behaves 


surface tension of 


‘negatively “in this region of concentration. 
At ethanol concentrations above 53 mole per cent, 
the 


octane ts the most volatile component and 
the 


mixture shows “* positive = properties, This 


conclusion Is contirmed ma spectacular manner 


by the distillation tests. Fig. 9 shows photographs 
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of the liquid-vapour behaviour for low and high 
ethanol concentrations at about the same vapour 
velocities. The tray efliciencies (30-40 per cent 
and about 70 per cent respectively) are in accor- 
dance with the formation of interfacial area 
observed. 

Absorption experiments with the Oldershaw 
sieve plates were carried out with n-heptane 
methyleyelohexane and n-heptane toluene mix- 
tures under conditions of increasing surface 
tension changes in the reflux flow, Table 6 gives 
the collected data. Tray efliciencies and visually 
observed foam heights of the absorption and total 
reflux distillation experiments are plotted against 
liquid driving force in Figs. 10 and 11, which are 
hence similar to the HETP graph given in Fig. 4 
for the Vigreux column. 

distinct the 


dependence of tray efliciency and foam heights 


similarity is noted between 
on the average driv Ing forces applied. The degree 
of foaming and increase in tray efliciency with 
the two mixtures, however, is clearly different and 
the 


both systems. This is in contrast to the Vigreux 


in line with surface tension differences in 


column, in which — with decreasing surface ten- 
sions and using a channelling mixture—no such 
influences of the system could be detected. 
Figs. 10 and 11 show that at low driving force, 
the foaming effect seems to disappear and the 
50-55 


In normal distillation practice these low driving 


per cent tray efficiency ts approached, 


forces are obtained at low or high concentrations 


or under pinched conditions. In these cir- 


cumstances, a loss in tray efliciency is therefore 


Table 6. 


vapour velocity 0-3 m_ see. 


HLARMENS 


Driving force, mole fraction 


Average tray efficiency,” 


Fic, 10. Influence of driving force on tray efliciency of 
Oldershaw sieve plates for various test mixtures. Vapour 
about O3m 


velocity 


benzene n-heptane 
benzene toluene 

n-heptane toluene 

A n-heptane methyleyclohexane 
Open points : normal distillation 


Closed points : absorption experiments 


60;-—— 


O20 


Driving force, male fraction 


Fic. 11. Influence of driving force on foam height on 
Oldershaw sieve plates with increasing surface tensions 


of retlux. 


@ n-heptane toluene 
n-heptane methyleyclohexane 
Open points : normal distillation 


Closed points : absorption experiments 


Absorption experiments with Oldershaw sieve tray column. Column with 5 plates, average 
Concentrations in mole per cent n-heptane 


Feed rates 
(mole /hr) (mole hr) 
System 


top bottom top bottom top 


n-heptane 13-8 100 


methyleyclohexane 
49 
: 
100 


n-heptane toluene 18-8 


Product rates Compositions, (mole per cent) 


Feeds 


bottom 


lverage 
Foam 
height 


(mem) 


Products driving Theoretical Plate 
- force plates efficiency 
top _bottom|(mole %) (%) 


20 
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to be expected when a positive mixture ts distilled, 
For total reflux distillation this is shown explicitly 
in Fig. 12. 

The concentration effect on plate efliciency as 
demonstrated in Fig. 12 has also been reported 


in the literature ; in all cases in which noticeable 


2 
* 40 
2 
8 Oo 20 40 60 BO 100 
Mole °lo volotiie component in mixture 
Fig. 12. Effect of concentration on tray efliciency of 


Oldershaw sieve plates in distillation with n-heptane 


toluene. Average vapour velocity 0-3 m/sec 


effects were recorded the binary system used can 
be identified as positive, Fig. 13 shows data of 
[10, 11] for the distillation of methyl- 
cyclohexane toluene and n-heptane methyleyelo- 


hexane in a small perforated plate column (38 mm 


» 
ote 
> 
20) 
| 
0 20 40 60 80 100 
Mole “/o volatile component in muxture 


Fic. 13. Concentration effect on tray efliciency of 


perforated plates observed by Thiissen. Average vapour 


ity about O-2 se« 


diameter). The behaviour of the first system was 
to be expected because of its high driving forces 
surface tension 


and large difference — in 


(3-5 dyn cm). It is of interest to note that also 
the xn-heptane-methyleyclohexane mixture, in 
spite of its low relative volatility, shows the 
concentration effects under extreme compositions. 
This mixture should therefore be considered as 


slightly positive instead of neutral. 


The influence of surface phenomena on the performance of distillation columns 


Tray efliciency maxima at mid-concentrations 
as found by Tutsssen and by the present authors 
for hydrocarbon systems have also been observed 
for ethanol-water mixtures [4, 5]. Because of the 
low surface tension and high volatility of ethanol 
this system should also be identified as positive. 

Further, some interesting cases are reported 
for the distillation of nitrogen oxygen mixtures. 
Nitrogen is the most volatile component and has 
a low surface tension at the boiling point (about 
Ydyn/em), whereas the heavy component oxygen 
has a high surface tension of about 13-5 dyn em. 
Accordingly, GRrASSMANN and Frank [3] found 
that 


bubbles into oxygen nitrogen liquid mixtures, 


when transferring oxygen from vapour 
higher transfer rates were observed than for the 


opposite case in which nitrogen vapour was 


transferred to the liquid. From the surface 
tensions of the pure compounds it is clear that 
in the first case positive surface tension changes 
develop in the liquid, which changes stabilize 
the bubble formation and hence favour mass 
transfer rates. 

The positive characteristics of the nitrogen 
oxygen system can also be derived from the 
distillation experiments with perforated plates 
earried out by LINDE (6). In the middle concen- 
tration range foaming was observed, whereas at 
extreme mixture compositions only a low spray 
bed deve loped, The same was found for the system 
nitrogen-argon, but foaming was absent with the 
These 


accordance with the surface tensions and relative 


mixture argon-oxygen. results are in 
volatil ties in the respective systems (Nitrogen 


argon : boiling point difference 10°C, surface ten- 


sions 9 and 12-5 dyn cm respectively ; argon-oxy- 
ven : boiling point difference about 1°C, surface 
tensions 12-5 and 13-5 dyn em respectively ). 
The effects reported above with perforated 
plate columns in general pertain to plates with 
very small open area and small perforations, 
conditions which are favourable for foam forma- 
tion. In technical columns, however, a higher 
free area and large perforations are used. There- 
fore, and also because of higher vapour velocities, 
a relatively larger amount of spray is formed and 
foaming diminishes. The results obtained with the 


spray column give rise to the expectation that 
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erage troy eff 


Ay 


Var ty, ry 


Tray eflicieney of a perforated tray, 16 per cent 


open area, 10 mm holes 


n-heptane toluene about 30 


n-heptane methvleyclohexane, about 350 


benzene toluene, about 350 


there would therefore be less difference between 


positive and negative mixtures with rv spect to 
tray ellicieneyv. This is indeed the case and can be 
Fig, 14 in 


observed for a perforated plate with 16 per cent 


seen trom which trav elliciencies 


open area, 10 mm hole size, (diameter 0.45 m) are 


plotted for three different svstems 
DiscusSION 
The 


paragraphs have shown that, depending on the 


experiments reported preceding 


type of equipment used, decreasing or increasing 


surface tensions of the reflux may influence 


interfacial area and separating power quite 


marke dly. \ survey of thn main trends obs« rve dl 


= 
Fable 7. 


ties pertaming to systems on 


Is given in For convenience, the proper 


which no surface 


tension changes occur (neutral systems) have 


been denoted in this table as “ unity 


Table 7. 


Survey of 


HARMENS 


In order to eive an explanation of these phe 


nomena, tt should be recalled that liquid surfaces 
of high interfacial tension contract when contacted 
This 

Manac 


7]. however, did not consider the influences 


with a surface of lower surface tension. 


is known as the “ Maragnoni effect 
NONI 
of vapour liquid) interaction his) studies, 

\ very well-known phenomenon which 
liquid film contraction oceurs as a result of mass 
transfer effects is the 


effect 


because of its 


so-called wie 
wall of the 
he na 


has been forme d. aleohol pre ferentially ev aporate Ss, 


vlass 


eres ps up the 


wetting film 


content and 
When 


wine creeps up the wall, and comes into contact 


with the latter 


leaving a residue of higher water 


therefore higher surface tension. fresh 


aqueous film, the contracts into 


small rivulets. 
\ similar explanation can be given for the 
taking 


mixture ina 


channelling phenomenon when 


place 


distilling a newative wetted wall 


type of column. Liquid films descending in such 


a column are never of uniform thiekness. and 


hence, locally. thin places will become more 


rapidly saturated with heavy component than the 
remainder of the tilm. Since the he avy component 
tension, areas of lower 
the 


tensions 


has the lowest surface 


surface tension develop locally remaining 


parts with higher surface therefore 


contract and the film breaks up into 


rivulets.* 


*kessentially the same explanation has been given by 


Bonp and Dowano (1) for their observations of the break 


down of water films in the case of ammonia absorption 


surface phenomena in distillation equipment 


Supported inele ial area 
Surface tension 
changes Separating 


Wetting powe 

no changes complete 
(neutral) 
increasing complete 
(positive) 
decreasing partial 


(negative) 


Effect of 


drit ing force 


nsupported interfacial area 
Type of Separating Effect of 


interfacial area power* driving force 


spray 
fourm 


spray 


*Separating power denoted unity for “ no change 
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Once the rivulets are formed, the edges, because 


of their smaller volume, more 
with the 


tension is thus lower at the edges and this stabilizes 


will always be 


saturated heavy component, Surface 


the contraction of the licquaicl (see sketch in Fig. 
15). 

It is easily understood that the differences of 
concentration in the liquid surfaces will become 


higher and contraction will become more severe 
if the 


The latter ts produced either bry an 


mass transfer per unit: area is’ higher. 


dl 


et Ss 


Fie. 15. Sehematic Hlustration of break-up of a liquid 


in oa Shaded arenas denote 


* negative svstem 


liquid of lower surface tension 


driving foree or because of higher mass transfer 
smaller HETP values). For the 
effect has 


influence 


coethcirents (i.e. 


Vigreux column, the driving fore 


already been illustrated in Fig. 4. The 
of the HETP on the severity of the channelling 
phenomenon Is recognized when COM paring the 
column 


results from the small concentric tubs 


(Fig. 1) and the Vigreux column (Fig. 2). Fig. 16 


100 


HE TP, wetting system 


Relation between loss in separating power and 


when distilling benzene n-heptane in various 


types of equipment, 
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gives a graphical representation of the relation 
between the ratio of HETP’s with channelling 
and complete wetting and the HETP of the 
columns with complete wetting. Furthermore, 
the same HETP ratios observed for the packed 
columns are plotted in this graph. 

It is thus seen that at the same HETP, the 
packed columns suffer a smaller loss im efliciency 
than the wetted wall type of columns. This ts 
easily explained by the fact that in the packed 
columns liquid is retained by capillary forces in 
interstices between the packing elements. The 
contraction of the liqquid surface in this case 
less effective and the loss tn interfacial 
Reduction of the 


channelling effect in this sense has also been 


will be 


area will therefore be smaller. 


found when prov iding the concentric tube column 
with a thin wire spiral. 

lhe distillation of positive mixtures in the 
Vivreux column and the Raschig rings column 
did not differences with 


reveal any significant 


respect to the distillation of a“ neutral ~ mixture 


in the same equipment, Considering for this 


case again the non-uniformities in the liquid 
films, it is now seen that thin and weak spots 
become preferentially saturated with the com- 
with the 


weak spots are thus reinforced, and as a result 


ponent highest surface tension. The 
the film tends to become of more uniform thickness 
and break-up will be prevented (see the sketch 
in Fig. 17). Wetting of the supporting solid 
surface with a film of the positive mixture will 


therefore be approximately the same as with 


Fic. 17. Sehematic illustration of the stabilization of a 
liquid film in a * positive’ svstem. Shaded areas denote 


liquid of lower surface tension. 


mixtures in which surface tension gradients are 


absent. 


holds 


good for supported interfacial area. In a system 


The latter conclusion, however. only 


where the vapour is dispersed into a liquid, the 
lifetime of the 
because of the 


bubbles formed may increase 


considerably reinforcement of 
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weak spots, These weak spots occur between 
adjacent bubbles (see Fig. 18); the preferential 


saturation of these spots with heavy component 
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hic. 18. Schematic illustration of the stabilization of 
vapour bubbles in a positive svater Shaded areas 


denote of lower surface 


causes the surface tension to rise locally. There 
fore, liquid is drawn in between thy bubbles and 
coalescence of the adjacent bubbles is prevented, 
If the changes in surface tension ar: irge enough, 
even a foam may develop in this manner. 

Since the degree of bubble stabil zation will 


depend on the change in surface tension of the 


liquid on the trays, the relation of foam height 
and tray efliciency to driving force as shown in 
Figs. 10 and 11 seems quite lo i 

Finally, the behaviour of neutral swstems and 
negative mixtures m a perforated plate column 
will bx briefly considered, kor hvdrocarbon 
mixtures in which no appreciable changes in 
surface tension occur, the stabilizing factor 
tending to produc prolonged bubble life is 
absent. Accordingly, bubbling is poor with such 
systems and spray formation occurs prefers ntially. 
From what has been said about the instability 
of liquid films when distilling a nm gative mixture. 
it will be clear that in this case the bubble stabi 
lizing factor is also absent. Approximately, the 
same interfacial area can therefore by expected 
for the two types of systems ; accordingly about 
equal tray efliciencies are observed for the 


‘neutral’ and me gative systems 


7. CONCLUSIONS 


The A ntal facts obse rved and the quah 


tative interpretation prov ided give rise to some 


definite conclusions regarding the formation of 
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and A, Harwens 


interfacial area in vapour liquid contacting 
equipment, 

The main conclusion is that changes in surface 
tension occurring in the liquid phase can have a 
spectacular effect on the interfacial area. Break- 
up of liquid films results from decreasing surface 
tensions and is of particular importance for 
wetted wall columns and packed columns. 
Stabilization of liquid films occurs when the 
interfacial tension increases ; even foaming ray 
deve lop in this case. The latter effect is of par 
ticular interest for the operation of plate columns. 
Since mass transfer rates are directly proportional 
to interfacial area, the surface tension effects 
are reflected in the HETP values of the various 
types of equipment. 

The above phenomena become especially notice 
able when the differences in surface tension 
between the components of the mixtures at the 
botline pomnt are higher than about 2Zdyn em 
and at driving forces above 5 mok per cent, 

\ further conclusion is that equipment with 
high mass transfer efliciencies is relatively mor 
sensitive to the surface effects than less eflicient 
apparatus. From a practical point of view this 
means that the effects are of particular Importance 
for laboratory apparatus, whilst technical equip 
ment will be much less sensitive to the surface 
phenomena deseribed. 

Finally, it is of interest to note that th 
influence of the changes in surface tension on 
mass transfer in general is quite severe and 
appreciably more pronounced than the effect of 
Variations in the static properties of the liquid 
phase, such as density, Viscosity and diffusivity, 
Trssen [10, 11], also Mousrap [8] et al., have 
in fact argued that variations in thes properties 
could not account for their observations. In this 
connexion it may be pointed out that when 
studies are made in which the static properties 
of the phases are the main subject, (e.g. when 
the contribution of the vapour and liquid phases 
to the resistances against mass transfer ts assessed) 
care must be taken to select experimental con- 
ditions such that the surface effects do not 
affect the observations. From the results of the 
present study it may be coneluded that such 


conditions cannot easily realized, 
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Abstract This article ithematieal treatment of the batch rectification of a binary 
lighter component in the distillate, both for 


mixture in the range of lo oneentrations of the 
the case in which the amount of holdup in the column may be neglected (Part 1) and for the 


case in which this holdup has te be taken into consideration (Part I It is assumed that at 


phase equilibrium the concentration of this component in the vapour ts directly proportional 
to the concentration in the liquid 
The relations derived in Part I show the importance of choosing a high reflux ratio in the 


process of removing the lighter mponent A graph is given which facilitates the finding of 


suitable values of the reflux ratio and of a suitable number of © transfer units” or theoretical plates 


for the column \ comparison is ide between removal of the lighter component by batch and 


by continuous reetifieation It ws found that, economically, batch rectification ts inferwr to 


continuous rectification if a high yree of purity of the heavier component is to be reached 


The amount of holdup is taken into consideration by assuming that the column is con 


tinuously in a state characteriz \ constant concentration ratios at a constant ratio between 


the amount of holdup and the a nt of material contained in the still rhis state is termed 


the stationary state of the second order The influence of the holdup is illustrated in a number 


of graphs, which bring to light s« ral particularities. Finally 
wh re termed the accumulation, accommodation, pinch and 


i discussion is given of some effects 


connected with the holdup, 
dead space effects (the first tw rrespond respectively to the depletion and flywheel effects of 


Picrorp ef. al. An example is given of a numerical computation 


Résumé On donne un traitement mathematique de la rectification discontinue d° un melange 


bunaire dans le domame de basses neentrations du composant le plus volatil dans le distillat. 


aussi bien dans le cas of la quantite de la retenue de la colonne est negligeable (Premiere Partie) 
que dans le cas ot il est pas ainsi (Deuxiéme Partie) admet qu’en cas d'équilibre 


des phases la concentration de « omposant dans la vapeur est directement proportionnelle 


a la concentration dans 


Il avere des relations dérivees dans la Premiere Partie que pour Penlévement du composant 


le plus volatil, on a intérét A choisir un taux de reflux cleve, Un graphique est donne en vue de 


faciliter la recherche de valeurs nvenables pour le taux de reflux et du nombre de * transfer 
units ou du nombre che pl t he rhyues de la cobonme l enlevement du composant ke 
plus volatil par rectification discontinue est compare a celui realise par rectification continue 


Cette comparaison montre qu'au | de vue économique, la rectification discontinue le cédea 


la rectification continuc agit d’obtenir une purete clevee de la composante la moin s 
volatile 


Pour porter en caleul la quantite de la retenue on admet que la colonne est continuellement 


dans un état qui est, le rapport entre la quantite de la retenue et le continu de la chaudiére étant 


constant, caracterisé par des rapports de concentration constants et qui est nomme etat station 


naire du deuxiéme ordre. Loinfluence de la retenue est illustree par quelques graphiques, ot: se 


réevélent différentes particularites. On traite entin quelques effets qui se rattachent a la retenue, 
effets dits d'accumulation, dav modation (ceux-ci correspondent respectivement a leffet 
d'épuisement et a leffet-volant de Piororp et collaborateurs), de “ pinch” et despace mort. 


On donne lexemple d'un calcul numerique 
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Zusammenfassung 


Der Aufsatz gibt eine mathematische Behandlung der absatzweisen-Rekti 


The final stage in the removal of the lighter component of a binary mixture by batch rectification 


fikation einer bindren Mis« hung im Bereich kleiner Konzentrationen der leichteren Komponente 


im Destillat, und zwar sowohl fiir den Fall, dass die Grosse des Betriebsinhalts (holdup) vernach 


lissigt werden kann (Teil 1), als auch fiir den Fall, dass dieser Betriebsinhalt berticksichtigt 


werden muss 


Es wird angenommen, dass beim Phas nyleichgewicht die 


Konzentration dieser 


Komponente im Dampf direkt proportional der Konzentration in der Fliissigkeit ist 


Die in Teil I abgeleiteten Bezichungen zeigen die Wi htigkeit der Wahl eines hohen Riick 


flussverhiltnisses, um die leichteren Komponente zu entfernen 


Kin mitgeteiltes Diagramm 


erleichtert das Auflinden geeigneter Werte des Riickflussverhaltnisses und einer veeigneten Zahl 


von Transferunits 


Kinheiten oder theoretischer Boden fiir die Kolonn 


Kentfernuny der 


leichteren Komponente bei der absatzweisen-Rektifikation wird mit der bei kontinuierlichem Betrieb 


verglichen. Esz eigt sich, dass wirtschaftlich geschen die absatzweisen-Rektifikation derkont inuler 


lichen Rektifikation unterlegen ist. wenn ein hoher Reinheitsgrad der schwereren Komponente 


erreicht werden soll. 


Der Betrag des Betriebsinhalts wird durch die Annahme berticksichtigt 


dass clie Kolonne 


stindig in einem Zustand ist, der durch konstante Konze ntrationsverhiiltnisse und ein konstantes 


Verhiiltnis zwischen dem Betriebsinhalt und dem Blaseninhalt gvekennzeichnet ist 


wird stationirer Zustand zweiter Art yenannt 


gramme erliutert, in denen sich 


verschiedene 


Dieser Zustand 


Der Einfluss des Betriebsinhalts wird durch Dia 


Einzelheite erkennen lassen. Schliesslich 


werden cinige Wirkungen diskutiert, die mit dem Betriebsinhalt in Zusammenhang stehen 


nimlich die der 


ersten beiden stimmen  iiberein 


Schwungradeffekt nach Pigrorp u.a.). Eine zahlenmiissige Berechnung wird als Beispiel veveben, 
- 


INTRODUCTION 

As a consequence of the non-stationary character 
of batch rectification computation of the process 
difficult the the 
holdup may be neglected and, consequently, the 


is rather even if amount of 
concentrations existing at a given moment may 
be equalized to those existing in a stationary 
state [12]. To account for the amount of holdup 
the aid of big computing machines [10, 13, 14] 
has had to be called in. The finding of a solution 
in the form of ordinary algebraic expressions 
seems to be fraught with unsurmountable 
mathematical difficulties [10]. Zorperwee [15], 
did 


empirical lines an expression for the so-called 


however, succeed in finding along semi- 
pole height, a measure of the slope of the recti- 
fication curve at 0-5, introduced by Bowman 


and Cicnetit [2]. This equation is of great 
importance in the choice of the right conditions 
[6]. Zumperwec further gave a relation between 
the pole height and the amount of the inter- 
mediate fraction between two symmetrical values 
of 2p (e.g. 0-95 and 0-05). This relation by itself is 
of no avail if one wants to know the amount of 
the lighter component which is left in still and 
column after the take-off of the 


fraction has been stopped, and which conta- 


intermediate 


A\kkumulation, der Akkomodation, des “ 
unterschiedlich mit 
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pinch und des toten Raums (dix 


dem  Erschépfungseffekt und dem 


minates the product obtained by simply emptying 
still and column or by distillation. 

One cause of the mathematical difliculties is 
the the 


equation. We have therefore restricted our con 


non-linear character of equilibrium 
siderations to the range of low concentrations of 
the the 


(e.g. ap 0-1), so that is becomes permissible 


lighter component — in distillate 
at phase equilibrium to assume the concentration 
of this component in the vapour to be directly 


proportional to its concentration in the liquid : 


y* (1) 


This renders possible a certain assumption, to 
be named in Part IT, by means of which a set of 
ordinary algebraic equations can be found. 
describing batch rectification in the range of low 
concentrations of the lighter component when 
the holdup cannot be neglected. Although in this 
way no direct contribution is made to the solu- 
tion of the general problem, we see the following 
advantages : 

(a) Calculation of the ratio between x, and x, 
(the mole fraction of the lighter component 
in the still) or ap) and & (the mean mole 
fraction in still and column). This may 


be used as an addition to the results 


| 
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obtained by Zurperwec. The method is 
the 


known e.g. 


further useful for 


s or a, if ap is from a 
measurement), 
batch rectification 


low 


the 


range of 


(b) Computation of 
process in the Up values, 

This computation is important in. those 

eases where the lighter component has 

to be thoroughly removed or where the 

available column has a low separating 
power, 

(c) Getting an insight into the influence of 


the amount of holdup. 


a general understanding and to 
for Part Il, in Part I 


calculations and considerations will be given in 


In order to get 


provide a background 


which the amount of holdup is neglected. In 


Part II allowance will be made for the amount of 


holdup. 


Before proceeding to the derivation of equations 
something should be said about the measure that 
will be used for expressing the separating power 
of a column 

The separating power of a column can_ be 
expressed in a number of theoretical plates (N) 
or in a number of transfer units (1). The number 
of theoretical plates of a plate column can easily 
be determined if the * Murphree plate efliciency 
is known [8]. If a column is used in which vapour 
and liquid are not contacted stage-wise, but con 
tinuously, e.g. by the use of a packing or the wall 
units is the 


of a tube, the number of transfer 


obvious measure to be used for the separating 
power [3). 


For N and [4] has derived equations 


which apply if the column is in a stationary 
state, the line of equilibrium satisfies equation (1) 
and the operating line is straight. In these 


A are related to « VL. and 


By dividing one 


equations N and 
the change in concentration. 
equation by the other Colburn finds* 


* Equation (2) can easily be derived (also for an equili- 


brium line of the type y* = a 


b) by calculating the 
number of transfer units corresponding to the concentra- 


tion change over one theoretical plate 


Burres 


calculation of 


A InakR 
N aR | 


(2) 


This equation renders it possible to assign a 
number of transfer units to a column of which 
the number of theoretical plates is known, and 
vice versa. 

It is likely that a plate column and a packed 
the 


(evenly distributed) and having the same separat- 


column having same amount of holdup 


ing power in the stationary state, will not behave 
batch This 


imphes that, by means of equation (2), results of 


very differently mn rectifications. 
calculations for a packed column may be applied 


to a plate column and vice versa. For this 
reason our calculations have been carried out for 
that type of column which seemed to involve 
the 


type in which vapour and liquid are not contacted 


the least mathematical difliculty, namely 


stage-wise and for which the number of transfer 


units is characteristic of the separating power. 


PARI 


COLUMN 


RECTIFICATION IN A WITH 


NEGLIGIBLE 
Derwation of equations 
Before proceeding to the derivation of equations, 


the following assumptions must be mentioned : 


(a) The mixture to be rectified is a binary 
mixture. The lighter component con 
stitutes only a small percentage of the 
distillate, so that equation (1) may be 


used to express the equilibrium relation. 
(b) The vapour leaving the still is in equili- 
brium with the liquid in the still. 
(c) The reflux has the same composition as 
the distillate. 
(d) The flows of vapour and liquid are the 
same In any cross-section, 
(ec) Vapour and liquid are not contacted 
The rate of 
the 
Co_purn [3] 


stage-wise. mass transfer is 


described by expression given by 


and (right-hand 


member of equation (4)). 
(f) Neither in the vapour nor in the liquid is 
there diffusion or the 


any mixing in 


direction of flow. 
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(g) The number of transfer units of the 


column does not change with time. 
(h) The reflux ratio does not change with time. 


If the amount of holdup may be left out of 


consideration, the material balance for batch 


rectification is: 
d(agW)+ ap)pdD = 0 


Wi D W, 
from which can be derived: 


d Inag As 
d In W Xs 


(3) 


in which equation X, has been written for 

* the ratio of the mole fractions of the lighter 
Tp 

component in still and distillate. 
X, must now be related to the relative volatility 
(x), the the 


transfer units (4). For the reader's convenience 


The quantity 


reflux ratio (r), and number of 
this will be discussed in some detail. 
Drawing up a material balance for an imaginary 


disk of a column provides the equation : 


dy 
dl 


O(y* — y) (4) 
The right-hand member is the expression for the 


rate of mass transfer given by CuHrLTon and 
CoLBuRN [3]. 


sionless quantity 


After introduction of the dimen 
A, detined by : 


da KO 
dl 


and termed the number of transfer units (after 
and CoLtnurn), the equation changes to : 


dy 
or, in the form of an integral, 


dy 


so as 


A 
* Purposely this is not expressed as j 


to leave open the possibility that KO Vis a function of 1. 


A material balance for the part of a column 
above an arbitrary horizontal cross-section is 
the 
applying to stationary states : 


given by following well-known equation 


Up 


Using equations (1) and (5) the integral may be 


worked out, which results in: 


l (x 


This equation results also from a conversion of 
an equation derived by Co_purn [4], 


The equation given above may also be written 


1) ajexp[—(aR—1)A] +1 


xR 


(6) 


that 
it is independent of xg or rp. We 


From equation (6) it is seen Xo is a 
constant, 1.e. 
remark that this conclusion could already have 
been drawn after replacing x and y by the ratios 
NX and Y in the equations from which equation 
(6) has been derived (this replacement 1s possible 
because of the homogeneous form of the equa 


tions). Equation (3) may now be changed to: 


1—Xs, 
and W 


already been derived for simple distillation from 
a still by Rayiercnu [11]. 


An analogous relation between ag has 


(a 1) (r + 1) 
Insertion of this expression for XY, in equation 


For A + o, 


(7) results in: 


1) In 


(x 1) (r : 


0 
This equation has already been derived, in a 
different form, by CoLBuRN and STearns [5]. 
It appears that equation (8) is already a good 
approximation when A has a relatively small 
value. Because equation (8) represents a removal 
of the more volatile component which is as 
rapid as possible at a given reflux ratio, it is in 
which is 


practice desirable to give A a value 
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sufficiently great to get a near approximation 


of equation (8). To survey quickly the require 


ments as to 1 we have drawn up a diagram 


(Fig. 1) which relates to each other .1, x r. and 
a quantity « which represents the relative devia 


tion of XN, from 


This diagram is based on the following equation : 


1) exp | 1) (Ga) 


which is another form of equation (6) and in 


which m and = are defined as 


x Xs 
Other variables could have been chosen instead 


But the latter have the 


that. if « is left out of consideration, the one (m) 


of m and 3. advantage 


is a function of A only, the other (2) of 1 only, 


The diagram may also be used in the calculation 


\ \ \\ 
\ 
AN 


YX 
AN 


\\ \\\ 


\\ 


Fic. 1. Nomograph, representing equation (6a) and relating 
to each other the ratio of the mole fractions of the lighter 


component in still liquid and in distillate (.\.¢), the relative 


volatility (a), the reflux ratio (r), and the number of 


transfer units (4). Significance of symbols 


of the factor (1 NX) NXg in equation (7), which 

(a + 1) 

following discussion an 


1s equal to 


In the approximate 


is needed. To derive such an 
0 


equation, equation (7) is written as 


equation for 


(x ltr 1) 
lr 


) 
If and are small, 


may be replaced 


by « and In =. by 


. Equation (7a) then 


results in: 


D 


W, (« 1) (r 


Discussion 


In the the 


aspects of the removal of the lighter component 


following discussion economic 
by batch rectification will take the central place. 
As a measure of the labour which has to be 
spent on removing a volatile component by batch 
negligible 


rectification with 


taken 


holdup, may be 
(r+ 1) DW, 
1—DW, 


The quantity FE is the amount of vapour that 
has to be formed per unit amount of product 
(the least volatile component), and is.c onse- 
quently, a measure of the heat consumption and 
of the time needed (if the diameter of the column 
is given) or the diameter of the column (if the 
time is given). Equation (9) leads to the following 
expression for E: 


In Te) 
E ~ (1 (10) 


x l 


0 

Usually, « and $ D/W, can easily be kept 
small, in which case E is not much greater than 
the smallest possible value. This is determined 
by In (ago “r7g) and « — 1. 
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Equation (9) shows that, to keep DW, small, 


a sulliciently high r value should be chosen. 


An example will be given here which proves that 


a high value of r is soon suitable: if one puts 

ty = 100,24 = 2 and D/W, = 0-1, 
from equation (9), disregarding the terms « and 
a value of 45 for r. 


If r exceeds a certain optimum value, E shows 


one finds 


an increase instead of a decrease. This is caused 


by a shortage of transfer units, as a result of 


which the decrease of } dD Wi, is smaller than the 


increase of 


Next a 


rectiheation 


batch 
rectification as 


comparison will be given of 
continuous 
the 


assuming this component to be present in low 


and 


methods of removing lighter component, 
concentrations. 

In a continuous rectification process 2 under- 
goes in a stripping section - which will, for the 
sake of convenience, be considered as consisting 
of theoretical plates a change per theoretical 
plate which is approximately given by the factor 
1)} not 
differ too little from the fraction of the lighter 
the 
quently, the change over N theoretical plates is 

xp |* 
approximately given by a factor ee To 


x L} xp prov ided does 


component in bottom product. Conse- 


give this factor a desired value, N and p must be 
balanced against each other. The result will be 
e.g. that for p a value is chosen which is 1-5 
times the minimum value, which is equal to 
1). different the 


mentioned factor will have little influence upon the 


value of afore- 


choice of p. It is seen that in the case of removal 
by continuous rectification N or A is strongly 
influenced by the requirement as to purity made 
The 


following antiparallel between batch rectification 


with regard to the heavier component. 


and continuous rectification will now be clear: 


batch rectification: E is highly and A only 
slightly dependent on the degree of purity 


of the heavier component required* 


continuous rectification: A is highly and p (com- 
parable with E) only slightly dependent on 


the degree of purity of the heavier com 


ponent required. 


The increase of EF with the degree of purity 
in batch rectification causes this process to be 
economically clearly inferior to continuous 
rectification if a high degree of purity ts required. 
In addition, in the continuous rectification process 
the stream of vapour leaving the stripping section 
can be utilized in a rectifving section for obtaining 
the lighter component in a pure state. All this 
alter the fact that the 


light component by batch rectification is. well 


does not removal of a 
feasible (if only the value of r is high enough) 


and is not as inapplicable as Luoyp [9] suggests. 


RECTIFICATION IN A COLUMN WITH 


NON-NEGI HoLtpupe 
Introduction 
In Part I it has been found that if the amount 


of holdup may be neglected, the ratio between the 


IGIBLI 


mole fractions of the lighter component in still 
and distillate is constant. Of course, this implies 
These 


ratios have other values if the amount of holdup 


constant ratios throughout the column. 
is not negligible, but may then be constant too, 
a condition being, as will be shown later, that the 
ratio h between the amount of holdup and the 
amount of liquid in still and column remains 
constant. In practice this can be realized by 
the 


characterized 


constantly supplying the still with 


This 


by constant ratios of the concentrations and 


pure 


heavier component. state, 


approached asymptotically as time elapses, will 
be called the stationary state of the second order. 
In ordinary batch rectification this state will, 
in theory, not be reached, as during the rectifica- 


Yet 


tion process h does not remain constant. 


* If it is assumed that ¢ is, at given A, directly propor- 


tional to 1 1 (generally a rough approximation), D/W 
is inversely proportional to r | (good approximation), 


r, at given A, has its optimum value, and dE /dA has a 
the 


A and In 2g) 


constant value determined by economical factors, 
following simple relationship between 
can be derived : 

d 1) Aj 


d In (In agg /2g) 
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it will be assumed that, if A is small, the concen- 
tration will always be that characteristic of a 
state of the By this 


simplifying assumption the mathematical treat- 


stationary second order. 


ment of the problem of the influence of the 


amount of holdup in batch rectification ts greatly 


facilitated 


Derivation of equations 


To the assumptions made in Part | the follow- 


ing relating to the holdup — should be added ; 


(1 The amount of holdup per transfer unit 
is constant (both as regards place (A) and 


time 


The holdup caused by the vapour may be 


disregard 


The holdup of the condenser may be 


neglected 


If there is an appreciable amount of holdup, a 


material balance leads to 


ryll) 0 


For ordinary batch rectification one must add: 


dW dD 0 


It is difficult to derive from this, for the 


quantity S detined by 


not 


dina*y 


(dD) M 


case of rectification with constant W: 


the 


l 


(11) 


and in the case of ordinary batch rectification : 


h) dX, dh 
X 


dD dh 
M h 


Xe h?dXy, dh (12) 


{use is made of 


*One may also write dinzy, din M If h = 0, 


dinzp dinzg and M 
as the left member of equation 3 


W. Then this form is the same 


In ordinary batch rectification the concentra- 
tion distribution will (if A is sufliciently small) 
the for the 
stationary state of the second order, the existence 
For this 


lag behind distribution 


slightly 


of which will be assumed however. 
reason there is no point in keeping strictly to the 
differential quotients in equation (12) and trying 
to meet serious mathematical dilliculties. There- 
fore, an approximation will be made. 

In numerical calculations it was found that 
from h 0, V. and X\y decrease nearly linearly 
with A over a comparatively large interval of h. 
Consequently, hdX. dh can well be approxi- 
erator of the right-hand member of equation 
(12) by writing: 1 + (dXy/dh 
dX. dh). The latter term of this expression is 


quadratic in A and may therefore, for low values 


mated by writing YX, » and the num- 


of h, be neglected with respect to Vy... ,. which 


term itself is often small with respect to 1, 
Therefore the following equation, which gives a 
good approximation, may be written instead of 
equation (12): 


0 


13) 


Ss 


The value of Xs. » follows from equation (6). 
Equation (13) becomes exact when r approaches 
infinity. 

In the case of rectification with constant W, 
the quantities S, XN. and X are constant, so that 
integration of the equation defining S leads to 


In 2p) = In ty) — (14) 


In the case of ordinary batch rectification, the 
quantities 8S, X, and X are not constant, so that 


the appropriate equation then ts 


D 


s dD 


lin X X,) 

can with suflicient 
S,, In(1 — D/M,), in 
which expression S,, is the arithmetical mean of 


Often the integral form 


accuracy be replaced by 


the initial and the ultimate value of S. 
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The quantities Y,, X and S must be related to 
x, r, Aand h. To this end the differential equa- 
tions describing the process in a column have 
to be solved. These equations can be derived 
from «a material balance for an imaginary disk 


of a column and read : 
KO (y* y) 
KO (y"* y) 


As appears from an article by Bowman and 


Briant [1]* and from another by Jaswon and 
[7]. 


obtain a 


the treatment of these equations to 


numerical solution is 
the 


equilibrium line is assumed, This is mainly due 


a complicated 


business, even if existence of a straight 
to the fact that an initial disturbance (c.g. taking 
off distillate after a long time of total reflux) 
returns periodically as a discontinuity in the 
course of the concentrations with time. However, 
the discontinuities are smoothed out as distilla 
tion proceeds (see Fig. 9 of Jaswon and Sarru’s 
article). In this article it will be assumed that 
the 


initial state and has attained to a state which 


column has completely “forgotten” its 
has been called the stationary state of the second 
order because of the constancy of the concentra- 
tion ratios. A condition of this constancy is the 
passing of a straight equilibrium line through the 


More- 
over the holdup due to the vapour will be 


origin (y* xx), which will be assumed. 
neglected. In this way relatively simple equations 
result. 

variables 


the dimensionless 


A and + (which is not new), the latter of which ts 


By introducing 


detined by : 


and by putting v, 0, the equations given 


above are changed to: 


* Bowman and Briant use a different expression for 
the rate of mass transfer. 


R(y* — y) 


dA 
‘ di 
Elimination of y and 


4 


equation has been added by differentiation to 


is now possible after an 


A of the first equation. The result of this elimina- 
thon ms: 


dy* 


0 (16) 


dA2 dx dA or dAdr 


If the 
must be equal to f(A) 


constant, 2 
Now X 
for f (A) and x, for g (+). Substitution in equation 


concentration ratios are 


£ (rT). is chosen 


dy* 
(16) | with 


a} leads to an equation for ap 


and one for XY: 


dinap 

dr 
di? 


dX 
dX 


(a R A (18) 
(4 is a constant). 
It may be remarked that the derivation of the 
equation for the +-independent mole fraction 
ratio XY would not have been possible, if the 
differential equation (16) were not homogeneous. 
It will appear from the derivation and use of 
equation (21) that it is also necessary for the 
relation between y* and x (equation (1)) to be 
homogeneous. 

The (16) 
incorporated in equation (17) and equation (18) 


solution of equation which is 


is of such a nature that the terms of equation 


(16), multiplied by have values which are 


a 

independent of +. Therefore, it has to be con- 
sidered as an asymptotic solution. It describes 
reached after an 
infinite time or which are maintained if already 
present at ¢ 0. 


column conditions which are 


A in equations (17) and (18) is a constant, 
i.e. a quantity independent of A and +r. However, 
A and h. From this it 
follows that, strictly speaking, these equations 


it does depend on a, 7, 
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apply only if during the rectification A remains 


We shall 


unchanged, icc. if W is ke pt constant 


iim’. 


ordimary bateh reetification alt hous is not 
1 constant m this case. This is ess ntially what 
was meant by saving that also in such a case a 


column ts ina stationary stats of the second order. 


It must be remarked that for an ordinary batch 
rectification this stats corresponds to another 
value of A than it does for a rectification with 
onstant W oat the same h. 


Application of the boundary ¢ onditions 
l for A 1 A l 
for A A da 


(this from 2 


I) R 


follows 


leads to the following solution for equation (18) 


a l 
\ : exp | a, (A \) | 
a a, 
2 
a. 
| a, (Al (10 
a, ay 
where a, and a, are the roots of th quadratic 


R l a 0 


are some simple re lationships bet ween a, and ay, 


equation a® There 


which are useful in the reduction of expressions 


contaiming a, and a, 


ay 


Starting from the equation for Y we can easily 


tind the expressions for VY, VY. and \, 


1 
| AN 
R d \ 
(a 1) a 
exp |—a,(A \) 
a, a, 
1) ad 
exp | \) (20) 
a, ay 
y 
— 
x 
l l a 
ay, 
x a, 2 
l | 
exp | a, A) (21) 


Bu rres 


however, these equations also in the case of 


\y Vda 
rA Be 
ArA | a, a, 
rA 1 
“ exp ( a, A) (22) 
a, a, 

As regards the parameter 4 figuring in the 
above equations, it follows from equation (17) 
and the detinitions of S and + that 

ArA AS (23) 

By elimination of A from th equations (23) 

and (11) or (13) the following equations for 8 


are obtained : 


in the case of rectification with constant W: 


ArAX 
ArA+Q (26) 
in the case of ordinary batch rectification : 
Xs 
@Q SAwO (25) 
For the quantity @Q imtroduced into these 


equations a relatively simple equation holds : 


(rd — a,) exp (— a, A) + (a, — rA)exp (— ag A) 
(2-1) jexp( —a,).1+ {a,—(a— 1) }exp(—a, A) 
— tanhg A (26) 
(x 1)! q| tanhq A 

ip 4 (a, q (a, ay) VP A| 
By means of the equations given above 
numerical values of S, X,, X,, and X can be 
caleulated if values of x, r, A and h are given. 
For this purpose the parameter A has to be 


assigned a certain value and S must be calculated 
by means of equation (24) or (25). From equation 
23) the value of h corresponding to the chosen 
value of A is then found. Subsequently, another 
value is to be assigned to A, etc. In the range of 
low values of A numerical calculations can be 


J. 
VOI 
9 
1952/ 
a, as rR l 
(l+a,)(1 + a.) 
“2 
112 


shortened by assuming VY, and X,, to be linearly 


dependent on h. 


Special forms of the equations 


The forms which the equations assume in 
boundary cases (r +0, r + A +0, A w) 
will, for the sake of brevity, be omitted, except 
in cases where they will be needed in the discus 
sions. They can be derived through the introduc 
tion of new parameters, being the limits of Ar 
or AA: in those cases where A -- 0 and A -> &, 
lim A A has to be introduced to denote the place 
in the column. It will be seen that S can be 
expressed as an explicit function of A, except 
for the case where A 0. 

Also for some specific values of A the equations 


assume special forms. 


(1) By 


obtained which apply to the case where 


assuming A 0 equations are 


h is equal to zero, which means that the 


holdup has no influence. 
(2) With A the equations for XY, ¥, 


assume the following simple 


S and h 
forms: 


1) A] 


These equations occasion the following re 


l [x (a 


marks : 
(a) The reflux ratio r does not appear. 
(b) X and ¥ are equal and have the same 
value as with total reflux. Y can be 
shown to be smaller than YV at values of h 
exceeding the one mentioned above. 
(c) Elimination of A leads to the equation 
S which 
the relation between S and A if r= 0 
(W constant). 

(3) Both for A (1 y«R) and for 
A=(l , « R)*, q is equal to zero. The 
equations applying to these cases must be 


equation also gives 


found as limits. For values of A lying be- 


tween the two values mentioned above, q 
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11 


is imaginary. However, after replacement 
y (4 
goniometric 


of q by ig (g p*)) and intro- 


duction of functions the 
equations assume real forms again. The 
case in which A (1 \ « R)* coincides 
with the case mentioned im (2) if 


’ ; it occurs at other values of r 
x 
if A and r is not too great. 
The case in which A = (1 y « R)* can 


be encountered with constant W only if A 
is very small (4 (38+ for 


1). 


The influence of the amount of holdup 


The influence of the amount of holdup is of a 


complicated nature. Various effects may be 


distinguished. Before discussing them we will 
first explain some figures illustrating the influence 


of the holdup. W is assumed to be constant. 


Fig. 2 gives X as a function of A for the case 
r A 
The curve for A 


10) at various values of A 
0 (h 0) 
The influence of A can 


and henee of h. 


shows the general course. 


Fic. 2. XN as a function of A at different values of A 
(x 1-5, 7 Ww, A 10). The values 0, 0-004, 0-008 
and 0-012 of A correspond with values of h equal to 0, 
0-077, 0-158, 0-390 respectively, for rectification with 
constant W. The right-hand side of the figure shows the 
curve for A 0, the left-hand side shows the deviations 
from this curve on an enlarged scale. 
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be seen from the deviation curves plotted on an 
enlarged X scale. The influence of the holdup 
on X, 4 and hence on Ny proves to be strong: 
decreases strongly as A and rise. 

A remarkable aspect of Fig. 2 is that the 
deviation curves for a value of A, which is 1-518 
smaller than A, have a common point of inter 
section lying on the A axis. At this value of A, 
NX has the same value as if there were no influence 
of the amount of holdup and from there the 
influence of the amount of holdup on X ts reversed. 
In reality the curves do not intersect in on 
mathematical point. For, numerical calculations 
show that (dX dA), 4 and (dX dA), 
are equal to zero not at the same value of A but 
at values differing by 0.00044. 

As V is a function of A \ Fig. 2 after 


x 


Fic. 3. S as a function of A at different reflux ratios for 
rectification with constant W (a 15; A 10). 


shifting of the zero-point of the A axis also holds 
for other values of A. It should be borne in 
mind, however, that at other values of A the 


values of A correspond with other values of h. 


Fig. 3 shows, for the case (« 1-5, A 10) the 
ratios between S and hk for a number of values of r. 

From this figure it is seen that below a certain 
value of r—lving between 40 and 80 - maxima 
occur, which are shifted to higher values of 
h as r decreases. A similar picture ts shown in 
Fig. 2 of Hourman and Husatn’s article [6). 
This figure applies to ordinary batch rectification 
in the case where zp) — 0.5*. It is seen that an 
analogy exists between the situation for rp, = 0-5 
and the situation for small values of rp. A 
difference is in the fact that the curves of 
Houtman and Husain have a_ practically 
symmetrical course near a maximum. Another 
and moreover remarkable difference is that in 
Fig. 3 the curves have a common point of inter 
section (P) (see p. 118, (2)), and that the series of 
maxima continues to the right of P towards 
higher values of S. 

To the right of P the highest value of S, 
viz, «(1 A), is reached by the curve for r 0, 
However, this case the maximum has 
degenerated to a discontinuity. Moreover, a 
study of the equations shows that the highest 
point (Q) corresponds with a situation in the 
column at which a 0 for A < A. This also is a 
degeneration. The same is true for the branch 
QR (a hyperbole corresponding to the equation 
AS x A). Here the limit of the assumption of a 
stationary state of the second order has been 
reached. It is assumable that this is indicated 
already in the region below QR, by the length of 
the time needed in rectification with constant 
W for approaching to the stationary state of the 
second order if the initial state deviates from it. 
This implies that the assumption of a stationary 
state of the second order existing in ordinary 
batch rectification is inadequate at high values of 
h and low values ot r. 

* The quantities plotted on the axes of Hourman and 
Hr figure are (twice the pole height) and 
OS 
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The two dotted curves starting from P enclose an area 
R)* (see p. 113, (3). 
One curve connects P with Q and is intersected by the 


in which A is greater than (1 


curves for values of r ranging from 0 to 2 ( 1 (« 1)). 
The maxima of these curves do not lie on the curve PQ. 
but slightly to the left of it. The other curve is intersected 
at a small angle by curves (not shown in the figure) for 


values of r ranging from 2 to about 3 


Fig. 4 shows the relation between S and h for 
a number of values of r at a small value of 


A (0-5); here also, « 1-5. For practical purposes 


this figure has little importance, but it is interest- 


ing from a mathematical point of view, 


After passing through the common point of intersection 
P, the curves appear to intersect again, thir time not in one 
point, but, for not too low values of r, still in a very small 
area (G). It is not accidental that in this area S is only 
little greater than a. The curve for 7 0 lies rather far 
outside this area. At the curve for 0-1 a shifting in the 
direction of the curve for r 0 is already perceptible. 
This shifting is accompanied by the occurrence of minimum 


values of S which are reached (for given values of h) 


at certain values of + 


Fic. 4. S as a function of A at different reflux ratios for 
rectification with constant W (« 15; A 0-5). 


Several aspects of the influence of the amount 
of holdup come to light clearly in the extreme 


cases where r = 0, r > «©. A 0. A + ow, which 


will be discussed below. 


The final stage in the removal of the lighter component of a binary mixture by batch rectification 
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r=0. If r — 0 and h may be put equal to 
zero the (stagnant) liquid in a column is in 
equilibrium with the vapour and has the same 
composition throughout the column. This changes 
as soon as h becomes greater than zero. For a 


given cross-section of the column the liquid 
composition then lags behind the vapour composi- 
tion which changes with time. Now the liquid 
tries to accommodate itself to the vapour by 
giving off volatile matter to the passing vapour. 
As a consequence a concentration gradient is 
maintained in the longitudinal sense of the 
column, on account of which NX, is lower than 


(ef. equation (27)). 


This is manifest in an increase of S 


This effect of the holdup, lowering XV, and 
called the 
It oceurs under all cireum- 


consequently raising S, might be 
accommodation effect. 
stances, except when r approaches to infinity and, 
state of the first 


consequently, stationary 


order results, and also when A is equal to zero, 


which means that there is no possibility of 
accommodation. This effect agrees in nature with 
the * fly wheel effect of Picrorp et al. {10}. 
r—»>e. In the case of r > © one is concerned 
with a stationary state of the first order. The 
change in the concentration from the still to the 
top of the column is then independent of h. 
However, there is an influence of A on 


R(= X,+h(Xy 


X to rise and hence S to decrease. This effect of 


Xx)); an increase of h causes 


the holdup, lowering 8S, which proceeds from the 
accumulation of volatile matter in the holdup 
h(X, — Xs)), called — the 
accumulation effect. In nature, it corresponds 
“ depletion effect ” of Pigford et al. 

effect 


circumstances, as X,, is always greater than XQ, 


(term might be 
with the 


The accumulation occurs under all 


which is a result of : 


(1) the rectifying power of a column (the 
only cause if r > 0), 
(2) the accommodation effect, 


(3) a third cause, which is the only one if 
A=0 (see sub A = 0). 


If A+ © orr = 0, Xy, at not too high values 
of h, is only slightly greater than X,, so that the 
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accumulation effect then has little signifieance 
(see sub A — ow). 

The accumulation and accommodation effects 
oppose each other. This is manifest in the 
occurrence of maxima in most of the curves 
shown in Fig. 3 and in the fact that these curves 
have a common point of intersection (P). 

A=0. In this case the column serves only as 
a temporary storage for the liquid fed in as reflux 
at the top, there being no effective contact 
between Vapour and liquuid, It is remarkabk 
that the concentration of the lighter component 
in the liquid now increases in the downward 
direction (at a given moment). This is caused 
by the steady decrease of the concentration in the 
reflux and by the absence of mixing in the direc- 
tion of tlow (assumption (f) p. 106), It will be clear 
that in this cas Xy Consequently, the 
accumulation effect occurs here. In view of the 
special cause of this effect if A 0 it might 
here be termed the dead space effect. 

The dead space effect also occurs in normal 
cases, when it may be considered as that part of 
the accumulation effect which is due to the 
tendency of the liquid contained in the column 
to lag behind in concentration. In normal cases 
there is, admittedly, no question of an increase 
in the downward direction of the concentration 
of the lighter component, but there is as a rudi- 
ment a rise of VY with respect to A near the 
top of the column (se¢ Fig. 2). 

The dead space effect disappears if A or r 
approaches infinity, and is felt more strongly 
as /1 or r decreases. It is clearly observed in 
Fig. 4, in that the curve for r — o is intersected 
for the second tim by the other curves at G 
and in that at high values of h, S is highest if r 
approaches infinity. 

A->®. In this case XY and Y undergo a 
stepwise change at the top of the column. It 
to 
1 and that, in the column, for small values 


can be proved that Y jumps from As ane 
of h, X increases approximately linearly from 
(1 hy) o (A 0) to Xsa (A + A) (this 
increase is caused by the accommodation effect). 
This means that at small values of h the main 
change in the concentration is at the top of the 


column, so that X,, differs little from Xx. At low 
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values of A the accumulation effect is conse- 
quently felt only slightly and the column behaves, 
as far as this effect is concerned, almost as if the 
holdup were part of the contents of the still. 

Equating of X and X., which implies neglect 
of a term containing h? (accumulation effect) 
with respect to a term containing h (accommoda 
tion effect) leads to: 


Ss = (27) 


\ derivation in which no approximations are 
used leads to the same equation, 

The case with r 0 (A finite) is alike in 
character to the case where A - ». If r 0. 


there is no concentration jump at the top of the 


column; and , do differ, however. Equa 
tion (27) applic provided h<A (A 1). 
The cases A + © andr 0 also agree in that 


if h approaches zero, the liquid assumes the 
same concentration throughout the column 
(except the top in the case that A + o) and th 
column is in a_ so-called “pinch” state. The 
weakening of the accumulation effect as a result 
of a small difference between X,, and VN. might 
therefore be termed the pinch effect. 

It will be clear that this effect is stronger 


according as 1 is greater and r is smaller. 


Vume rical ram ple 


Given: a column with five transfer units. 
by means of which, in an ordinary batch rectifica- 
tion at a reflux ratio of 12, a substance is removed 
for which « 2; the mole fraction of this sub- 
stance in the distillate is equal to 0.02 at the 
moment when A 0-1, 


Find : 
(1) the value of & at this moment ; 


(2) the value of % reached when 5 per cent of 
the still contents has been distilled off. 


After 5 per cent of the contents of the still 
has been distilled off, A is equal to 0-1: O95 
0.105263, Consequently, calculations have to be 
made for h 0-1 and h 0.105263. Calculation 
shows (a choice of A values has been made with 
the aid of preceding calculations) : 


for A 0.022362 : for A = 0.028862 : 
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The final stage 


No = 0.052595 0-051722 (eq. (21)) 
X yy = 0-221725 


13.3569 


0.221349 (eq. (22)) 
38-3354 
Ss 13-3167 3318 (eq. (25)) 


h 0.100755 0-105141 (eq. (23)) 


whence it follows for dX 


toh 0-1: 


dh and dX,, dh in a 


range close 


dX, dX y 
0.1990 
dh dh 


0-O857 


From small extrapolations it now follows that: 


or ft 


Xy 21790 


J 
NX. = 0.051698 

for h 
= 0.221839 


This leads to: 


0.105268 


0.069649 7} 
>for h 
S= 183141 | 


0-105263 


X = 0.069555 
for h 
Ss 13-3321 | 


0-02 
is equal to 0.02 
0.0013930 (answer to point 1). 
Applying equation (15) and considering S to b 


At the moment when ap and kh = 0-1, 
X — 0.069649, so that 2 
0.069649 
value 13-3231, finds : 


a constant having the one 


0.50491 


so that at the moment when h 
(i.e. after 5 per cent of the 
been distilled off) z 
0-069555 : 0-069649 
point 2). 

Remarks. (1) From the values found for 
dX./dh and dX,, dh it follows that the value of 
the neglected term h? (dX, dh dX. dh) (see 
p. 110) is 0.001133. As this term is much smaller 
Xs 0 it 
neglect it. The 


0-1052638 
still contents has 
= 0.50491 0-0013930 


0.0007024 (answer to 


than was, indeed, permissible to 


numerical results show more 


decimal places than is permitted in view of the 
neglect of this term, however. 


in the removal of the lighter component of 


a binary mixture by batch rectification 


(2) That 
change with h from h 
the fact that the 
h Otoh 
tively, 


X, and X,, show almost linear 
Otoh 


decrease of 


0-1 is seen from 
X, and X,, from 
0-1 is 0-01994 and 0-00836 respec- 
numbers, after 
little 


and 


which 
deviate but 
dX 


also shows a linear change. 


division by 0-1, 
from the values found for 
dX,, dh. The parameter A 
From the 
of X it is clear that a similar remark does not 


definition 


apply to this quantity ; it tends to a minimum 


(and hence S tends to a maximum). 


NOTATION* 


roots of the quadratic equation 
a® —(aR | Aja 

amount of holdup as a fraction of the amount 
of liquid contained in still and column, 

HOM 
distance from the bottom end of a column 
(metres) 
(a 1) A 


\ (p® — A) 


(a, 
$ (a, Ay) 
L 


reflux ratio, 


time (seconds) 
amount of holdup formed by liquid per metre 
of column (moles, metre) 
amount of holdup formed by vapour per 
metre of column (moles /metre) 
mole fraction of the lighter component in the 
liquid 
mean mole fraction of the lighter component 
in the holdup formed by the liquid 
(1 h)agt+ hay 
mole fraction of the lighter component in the 
vapour 
mole fraction of the lighter component in the 
vapour that is in equilibrium with liquid in 
which the mole fraction of the lighter com- 
ponent is equal to a. 


* For the sake of convenience amounts are expressed in 
moles and consequently fractions are considered as being 
mole fractions. One is free to use other units for expressing 
amounts (the use of a volumetric unit requires the intro- 
Instead of 
unit length”’’ and 


duction of a proportionality factor). * meter” 


and second’ may be read “ unit 


time.” 
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parameter 

amount of distillate (moles) 
D D 

1) { 

We, 


amount of holdup (moles) 


mass transfer coeflicient (moles m? sec) 

rate of flow of liquid through horizontal cross- 
section of a column (moles sec) 

of 


column (mole Ss). 


contained still and 


i 


number of theoretical plates 


amount liquid 


area of interface between liquid and vapour 
per metre of column (metres) 
ArANXy, 


Indices 


D 
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Mass-transfer from single solid spheres—I 


Transfer at low Reynolds numbers 


F. H. Garner and R. B. Kerry 


(Received 13 January 1958) 


Abstract 
acid spheres to be measured. Results are obtained for the overall and local mass-transfer rates 


\ low-speed water-tunnel is described to enable dissolution rates of pelleted organic 


from j in. diameter benzoic acid spheres to water at 30-0°C between Reynolds numbers of 2-3 
and 255. Comparison with similar data enables the upper limit to be extended to 900. Results 
show that the bulk flow depresses the mass-transfer rate over a certain flow range and that 
free-convective effects are not entirely absent until Np, = 750. This effect is the more pronounced 
in downflow. At Reynolds numbers greater than 250 the overall transfer results tend towards 
the theoretical relationship 
N 2A ga! 

Résumé L auteur decrit un tunnel a eau, de faible vitesse, permettant la mesure des vitesses 
de dissolution de sphéres d’acide organique. Les vitesses de transfert massique, globales ct locales, 
entre des spheres dacide benzoique de de diamétre et a 30-0 C, ont été obtenues 


pour des nombres de Reynolds de 2-3 a 255. La comparaison de ces mesures avec d'autres 
résultats similaires a permis d’extrapoler la valeur limite du nombre de Reynolds jusqu’a 900. 
Les résultats montrent, que dans un certain domaine écoulement diminue la vitesse de transfert 
massique et que les influences de la convection libre n'ont enti¢rement disparues qu’au dessus 
dun nombre de Reynolds de 750. Cet effet est le plus important pour écoulement descendant 
Pour les nombres de Reynolds supérieurs a 250, les résultats du transfert global tendent a obéir 
a Téquation théorique 
Ng, = 0-04 Np,' ? 

Zusammenfassung Zur Messung der Auflosungsgeschwindigkeit von kompakten Kugeln 
aus oganischen Séuren wurde win Wasserkanal fiir geringe Geschwindigkeiten errichtet. Gemessen 
wurde die gesamte und die Ortliche Stoffiibertragungsgeschwindigkeit von Kugeln aus Benzoe- 
siure von 19mm Durchmesser an Wasser von 30°C im Bereich der Reynoldszahl von 2-3 bis 
255. Durch Vergleich mit Abniichen Messungen konnte die obere Grenze auf 900 ausgedehnt 
werden. Die Versuchsergebnisse zeigten, dass die HauptstrOmung in gewissen Geschwindig 
keitsbereichen die Stoffiibergangszahl herabsetzt und dass bis Re 750 auch der Einfluss der 
freien Konvektion nicht zu vernachlissigen ist, insbesondere bie abwirtsgerichteter StrOmung. 
Bei Reyvnoldszahlen tiber 250 erreicht die Stoffiibergangszahl die theoretische Beziehung 


INTRODUCTION direction of decreasing concentration. This is 


OPERATIONS in which one fluid phase is dispersed 
throughout another are common and the limiting 
case of infinite drop viscosity (a solid dispersed 
phase) is of some importance. This limit is 
approached, whenever droplets are contaminated 
with surface-active compounds, which include 
many greases. 

Whenever there is a concentration gradient, 
there is a migration of ions or molecules in the 


known as mass-transfer. If such a process takes 
place uninfluenced by extraneous forces, this 
movement is termed molecular diffusion. In 
practice the differences in concentration result in 
differences of density, which cause a flow due to 
gravitational forces. This motion convects 
material from the phase boundary and so pro- 
motes mass-transfer. Such a flow is known as 


free or natural convection. The solute laden 
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material at the interfacial boundary will either 
fall or rise, depending upon whether this solution 
is heavier or lighter than the solvent. There may 
be a relative motion between the drop and the 
continuous phase due to centripetal, gravitational 
or inertial forces on the drop and the inertia of 
the fluid. 
of motion can be simplified by rigidly supporting 


For solid drops, analysis of this type 


the drop and allowing a controlled flow to wash 
The latter called 
convection. 

While some attention has been given to these 


over it. motion 1s forced 


processes separately, there ts little information 


about transfer when all modes of transfer ar 
equally significant. In general, the 
material transferred by molecular diffusion alone 
is small. The present paper will deal with some 
aspects of the proble m of transfer when both free 
and forced convection occur. The work in this 
paper is confined to spheres, which approximat« 


to drop shapes. 


TuEeory 


An estimate of the rate of mass-transfer in 


gives a datum for com- 
Pre- 


rie thods are time consuming and 


forced convection alone 
parison with the total mass-transfer rate. 
vious “ exact 
demand tedious numerical integrations [6, 26). 
Others based on the concept of the boundary- 
layer (see below) are limited in SCO The only 
method u 


rates from the whol 


available to predict mass-transfer 
assumes that this 
layer is constant throughout the wake behind 
Revnolds 


consists of a 


surface 


the sphere. Since at low numbers 


(Np. $50) the 
vortex [8], this assumption cannot be sound. A 


wake toroidal 
new method is therefor proposed : 

\ mass balance over anv volume of an incom- 
pressible fluid, in which there are neither sources 


nor sinks, gives 


0. (1) 


Newton's law of fluid stress holds, the 


When 


equation of motion becomes for steady flows: 


U.VU=F U (2) 


Equations (1) and (2) have no known elegant 


Gagner and R. 


amount of 
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solution. These equations may be simplified by 
the concept of a boundary-layer over a surface. 
It is postulated that all 
confined to a thin zone near the surface and the 


effects are 


Viscous 
llow outside this layer is inviscid. It is assumed 
that the thickness of this laver, 4 is small com- 
pared to the characteristic dimension of the 
order 5 are thus 


surface and quantities of 


neglected. Equations (1) and (2) then become : 


oa oY 
ou vil 
oy* 


been solved by Tomorika [22] for 
that the 


boundary-layer is described 


These hay 


flow around a sph re. He assumes 
motion outsick the 
by inviscid flow over a spherical surface. The 


solution may be written 


Z 


where 
du 
Z My, My . ete 
da 


and g (A), & (A) and A* (A) are algebraic functions 
of A. 
before separation of the forward-flow, but begins 


Equation (4) can be used over the surface 


to break down when the flow in the boundary- 
layer decelerates. 

Certain theoretical considerations [12] indicate 
that the skin-friction in the back-flow region ts 
range: 40 < Np», 1000 


small over the flow 


Fic. 1. Variation of skin-friction over sphere surface 
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(Fig. 1). If the skin-friction in the wake is 


neglected, then, assuming a quartic velocity 


get: 


distribution in the boundary-layer, we 


2+ Alp 


It follows that A 12. (t) 


This is a sufficient solution of wake flow for the 
present purposes. 

For solute of low solubility it is presumed that 
effect on the 


hydrodynamics of the problem. Thus equations 


the diffusion has no significant 
(3) and that of mass-transfer can be regarded as 
independent equations. 
It is assumed there are no ionic effects and the 
diffusion is isothermal and isobaric. If the solvent 
ve locity is small, an analysis of the resultant flux 
due to diffusion over an elemental volume of fluid 


under steady conditions is given by : 


U.V. 


v? Ce 7 


To solve the problem, following the method 


suggested by Grarron [1]], we assume all 


a thin layer 
detine the 


concentration changes occur over 
near the surface. For convenience we 


variables : 


It is assumed that the concentration distribution 
in this laver can be given by a simple algebraic 


expansion, VIZ: 


the first five terms are taken. 


coefficients, d,..... a,, are 


In practice only 
The 


boundary 


found from 


conditions: four of these can be 


written down straightway : 


At » 
Ar 
At» 0. 
ub 


The tifth condition is found from equation (7). 


At the surface this becomes : 


Mass-transfer from single solid spheres 


I. Transfer at low Reynolds numbers 


25a) 
oY 0 oy” 


Replacing pv by m, the radial mass flux per unit 


u 0 


surface, we get 


Whence the 


unknown coetlicients become : 


a 
2 
Me 
12 27,4) 
a. 
Me Mo 
Mo) 
y= - 
6 Mo 
where i, 5 yy. 


The rate of dissolution, m, is given by 


nee 
Me 311 (wp (Ss 
3 


The positive root is extraneous, as this implies 
that mass is being transferred to the sphere. 


For solutes of low solubility equation (8) can be 


simplified by retaining the first two terms of 
the binomial expansion, this gives : 
Me 2 (vip y,). 


To solve equation (9) we must evaluate the 
thickness of the mass-transfer layer, 5,,. Formal 
calculations [14, 15] for a flat plate and a rotating 


with the condition that the mass- and momentum- 


disk give dy, This expression agrees 
transfer film thicknesses are equal at Ny. = 1. 
It is unlikely that such a relationship depends 
upon the geometry of the system and this is 
presumed to be valid in the present case. 

The hydrodynamic boundary-layer thickness 
can be expressed in terms of a single flow para- 
meter by: 
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(10) 


A Uy. 


Following Tomotika’s method, we find the varia- 
tion of uy with angular distance from the forward 
stagnation point from the velocity potential of 


a sphere in inviscid flow : 


U reos | 1 + 


whence sin @. 


and from equation (10) 


382 U 
A cos Np, Ny cos @. (11) 
¢ d 
By definition the mass-transfer coetlicient 
K, is m/p — and thus by combining 
equations (9) and (11), we obtain the local 


Sherwood number at an angle @ from the forward 


pole : 


To obtain the average mass-transfer rate over a 
0, to 


band from @ @,. equation (12) is 


integrated over that surface : 


[Nou] 


27 a* sin 6 dé 


12 \ 3 Nae’ } Os 
| 27a*sin dO 
de 
A 
é 3 N 2 N { 13) 


In this manner overall transfer rates and those 
for the 
calculated 


forward-flow and wake areas can be 


using values of A caleulated from 


equations (4) and (6). The integrands can be 
evaluated graphically and such values are com- 
pared to and 


those GARNER 


[9]. 


presented by 


These obtained 


SUCKLING 


were experi- 


GARNER and R 


12 


B. Keey 


mentally in a 3in. diameter, horizontal water- 
tunnel over the range: 80 < Np, < 900, 
the recalculated the 
constant Schmidt number as used in the experi 


when 


results are on basis of 


mental section. A rectilinear correlation between 
Ny, and 4/N,», is found, and thus forced convec 


tion is assumed to predominate in this case. 


Table 1. 


tical rates of mass-transfer in forced convection 


Comparison of experimental and theore 


Nsh Nie Nise 


theoretical experimental 
Area (T) 

overall 
forward-tlow 1-025 1-08 
wake TRS 0-67 it7 
forward stagnation 

point 1-59 1-68 65 
rear stagnation 

pont 13-0 
ratio, forward (wake 28 120 


Fic. 2. Theoretical local mass-transfer rates in forced 


convection, 


the values for 


fortuitous, 


The close agreement between 


overall transfer is somewhat as 


theoretical values compared to experimental 
ones are low in the forward-flow area and high 
in the wake. The distribution of local transfer 
rates is given in Fig. 2, Within the rear hemisphere 
there is a pronounced minimum, which corres 
Nevertheless the 


relationship provides an adequate 


ponds to the separation ring. 
theoretical 
basis for comparisons with further experimental 
work. 
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Mass-transfer from single solid spheres—I. 


APPARATUS 


Most water-tunnels have been designed pri- 
marily for hydraulic studies on cavitation [16]. 
Applications of such tunnels include free-stream 
propeller testing, the study of hydrofoils and the 
Recently 
tunnels to study the dissolution of solid spheres [7] 


hvdroballistic design of missile noses. 


and to investigate diffusion-controlled chemical 
reactions [12] have been reported. 


Arift 


Fic. 3. Simplified diagram of water tunnel. 


The plant consists of two water-tunnels placed 
endwise so that there would be similar hydro- 
dynamic conditions for both upflow and downtlow. 
It is 
designed for flow rates between 1-3 and 150 gal hr 


A flowsheet of the plant is given in Fig. 3. 


and working temperatures between 50° and 100°F, 

Water from the recirculation system discharges 
into one of the surge-tanks, fabricated from 
spun 18 $.W.G. copper sheet and holding about 
20 gal. This allows viscous effects to damp out 
eddies in the incoming fluid. The inlet water is 
distributed through concentric rows of orifices 
and a 200 Further- 


more by reducing the fluid velocity to nearly 


200 mesh bronze screen. 


zero, the increase in kinetic head athwart any 
cross-section is almost constant in the transition 
from the surge-tank to test-section diameter. 
The contraction is therefore smooth and of such 


Transfer at low Reynolds numbers 


a form that there is a uniform pressure distribu 
tion [15]. The 
honeycomb, consisting of 2 in. 


water then passes through a 


10mm _ hexa 
fluid 


A honeycomb 


gonal brass tubes. From these. filaments 
emerge parallel to the tube walls. 
like a screen of large free cross-sectional area, 
reduces the intensity of any upstream turbulence, 


but increases the scale of it. 


+ 


Half-section on XX 
Fic. 4. 


Test section. 


The test-section (Fig. 4) was designed for ease 
of dismantling for cleaning and for good observa- 
the field. This 
essentially of a 12 in. 6 in. Pyrex glass tube 


tion of flow section consists 
held between two end-plates by ¥ in. diameter 
tie-bars. These also secure 3 in. Perspex sheets, 
which form a jacket round the test-section so that 
refractive distortion can be eliminated. 

throttling the 


by-pass valve and, at low flow rates, by reducing 


Flow control is achieved by 


the voltage input to the pump motor with a 
For all the 
controllability, expressed as the average standard 
from the 


variable autotransformer. runs 


deviation mean flow rate. was 


5-1 per cent. 
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Due to the slow circulating rate temperature 
control was fairly difficult. In the lower surge 
tank, four 2 kW kettle heaters are let in symm 
trically together with a bimetallic temperature 
controller. The latter is connected in parallel to 
one of the heaters and the contacts of a hot-wire 
a singh 


vacuum-switch. Cooling is effected by 


pass, finned exchanger, having twenty-four } in. 
in. OLD. tube. the heat 
the column, all the 


pipes are cellular 


fins on a lo reduce 


flux through walls of the 


principal covered with a 


insulation of thermal con 
0-05 hrft F. tanks 


are lagged with a magnesia-asbestos cement and 


asbestos-aluminium 


ductivity surge 


coated with retarded hemihydrate gypsum plaster 


to give a smooth finish. It is evident that the 
temperature control depends upon the flow rate 
and ambient 


temperature. For example, the 


temperature range 1s 0-1 C at a circulation 
rate of 125 gal hr and a working temperature of 
30° The 
all runs with up-flow was 0-2 ( 


(86°F). average standard deviation for 


and for those 


with down-tlow 


EXPERIMENTAL 


Materials used in the making of solid spheres 


are usually crystalline. Compression of the solid 


is preferable to casting, as the latter process 
tends to give diffential orientation to the crvstals 
[17]. The material to be pelleted was « ompressed 
in a tly-pre ss, which gave a 


t 


pressure of 5 tons in®. Excess solid was trimmed 
from the ejected pellet and it was re¢ ompressed 
but 
later alone. The die and punch holders are shewn 
in Fig. 5 and the faces AA‘ are 


the punch and die fac 


several times, at first with extra material. 


machined so that 
parallel. \ 


wives a 


remain 


clearance of about 0-008 in. 


flash, 


machined surfaces. 


minimum 
the 
It was satisfactory to use the 
by the 


consistant with avoiding damaging 
materials as 
the 


pelleting properties [17]. 


supplied manufacturer, 


since random size distribution gave good 


In all cases anv diameter 
of the pellet was within 1 per cent of the mean 


diameter and, in general, pellets of uniformly 


high density could be obtained to within 1 per cent 
of the commonly accepted value of the solid density. 


To minimize disturbance to the flow. each 


Garner and R. B. 


Punch -hoider 


to fit into slide Allen bolt 


OO 


Ejector 


hic. 5. Die and punch to make half-in. spheres. 


sphere was supported at the 
The 
O-O1L5 in. 


rear Stagnation 


point | 10}, one chosen consisted of a 


1} im. diameter nichrome wire, 


soldered into a wv} ds in. diameter brass rod, 


which was bent into an ell and fitted into the 
support holder (Fig. 4). 

Before the spheres were mounted in the test- 
section, they were dipped into Teepol-B to ensure 
that the 
This procedure is claimed to have no apparent 
effect the (13, 19}. 


During each experiment the temperature and 


surface became completely wetted, 


upon rate of dissolution 


rotameter reading were noted at regular intervals 

and weighted mean values were used in evaluating 

the phy sical properties of the system. 
Dissolution of the spheres was studied photo- 


graphically using a double-extension bellows, 


plate camera of 170mm focal length in con- 


junction with a supplementary lens of — 8 diopt. 
A uniformly-lit, matt background was used. 

A number of jin. diameter benzoic acid 
spheres were prepared with a little red azo-dye, 


Rhodamine-B, the 
compression, As orthochromatic 


added in final stages of 


material is 
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Mass-transfer from single solid spheres —I. Transfer at low Reynolds numbers 


Reynolds No as parameter Reynolds No as parometer 


180 


Local mass-transfer rates. 
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blind to red, this was used to obtain photographs 
of the dye streaming from the sphere surface. 
Dye traces taken from such photographs are 
presented in Fig. 6. 

In the quantitative work one exposure was 
made at the beginning and end of the dissolution 
on panchromatic plates. The processed negative 
plates were projected by a diascope, giving a 
resultant magnification of 20 and the outlines 
of the pellets drawn out. These were superim- 
posed and the diminution of radii at successive 
then the 


silhouettes. average 


10° intervals was measured around 


periphery of the For an 


diminution of 0-020 in., this change could be 


measured to within 5 per cent, 


Mass-transfer coeflicients could then be calcu- 


lated readily from this diminution (/), since 


A, S 
t 
hop 
(14) 


Col 


and overall coeflicients are found by integrating 


over the surface : 


A, | Asin 0 tra 
0 


h sin dé 
2 Col 


The integrand can be evaluated graphically. 

To solve the proble m we need the following 
data : pellet density, the solubility and diffusivity 
of the solute in the solvent and the viscosity of 
the saturated solution. 

Benzoic acid of “ Analar”™ grade was used, 
as this was available commercially in high purity 


Pellet 


density was measured by Archimedes’ principle. 


and possesses good pelleting properties. 


Solubility of benzoic acid in water was calculated 
from Bourcorn’s data [3]. Due to the paucity 
of experimental data, liquid diffusivities are 
usually estimated from semi-empirical correla- 
tions. The one chosen is that proposed by WILKE 
and Pix CuanG [23], as this makes the most 
for solvent association. 


satisfactory allowance 


The viscosity of the solutions are assumed to be 
the same as that for pure water. 

The Reynolds number of the sphere, U,,d v 
varied from 2-3 to 255 and all results are expressed 
on a constant Schmidt number basis (Nx. = 788, 
T = 300°C). Other 
Reynolds numbers in horizontal flow are recal- 
Both sets of data are 


data at slightly higher 
culated on this basis [9]. 
plotted in Fig. 8; local mass-transfer rates are 
given in Figs. 9 (a, b and ¢). The position of the 
minimum mass-transfer rate is compared to the 
locus of separation, which was found by injecting 


red ink into the boundary-layer flow, (Fig. 10). 


DiscUSSION 


The pellet surface after each experiment was 
free from pitting, which suggests that the spheres 
dissolved rather than eroded away. This is due 
to the low voidage in the compressed spheres. 
The curves for the overall transfer rate (Fig. 8) 
show a high degree of correlation with low mean 
probable errors (downtlow 6 per cent; 
uptlow 7-6 per cent). These errors are of the same 
the method of 
estimating the diminution of radii. 

Below the Sherwood 
to be the same as that for molecular diffusion 
At higher Reynolds 


number the mass-transfer rate is depressed and 


order as those introduced in 


number seems 
and free convection alone. 
20 the curve for uptlow falls above 


Such a 


separation has been noticed previously in experi- 


above Np, 
those for horizontal- and down-flow. 
ments on packed beds [2a]. Since the data for 
horizontal flow fall between those for up- and 
down-flow, the latter can be justifiably extra- 
polated to their upper point of intersection, which 
750. this 


750 forced convection alone is 


appears to be about Np, Thus, in 
case, above Np, 
significant. 
Since the theoretical equation (13) agrees very 
well with the experimental curve for forced 
convection and that for horizontal flow, it seems 
that free-convective effects can promote or hinder 
Con 
versely a bulk flow can hinder mass-transfer by 
free This 
transfer rate will be termed “ interference.” 


mass-transfer rates in forced convection. 


convection. mutual depression of 
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It might that the convective 


velocity could be obtained by adding vectorially 


appear true 
the velocity associated with the free convection 
to that of the 
experiments with dye-containing pellets showed 


forced convection. However, the 
that this transitional region is one of metastability. 
Below Np, = 120 the steady wake behind the 
sphere began to break up and the dye effusing 
from the wake in upflow tended to drift down- 
stream at about one sphere diameter from th 
surface. Typical stages of the unsteady motion 
of the wake are given in Fig. 6 of a sphere of 
benzoic acid, in which a little dye had been 
impre ssed in the surface. As the caption indicates, 
these pictures are not strictly comparable, but 
they give a good indication of the changeover 
from a wake shewing a toroidal vortex to a 
It is thought the latter is 
flow. As 
the flow rate decreased the latter flow pattern 


tundish-shaped tail. 
more characteristic of free-convective 
was more persistant, but unsteadiness was still 
detected at a Reynolds number of 10, It is 
that 
capacity lag (23) and possibly the aforementioned 
the 


known free-convection 1s subject to a 


motion is due to inherent transients in 
transport process, 
the 


curve 


fall on a 


limits where forced 


transfer results 
the 


and free convection predominate, it follows that 


Since overall 


smooth between 
in this transitional region there are time-averaged 
values of the variables, which define the process. 
With this in mind we can examine the flow field 
and give a resultant velocity to each point within 
it. 

In understanding this 


the which 


interference occurs it is necessary to have some 


way in 


knowledge of the flow patterns in free convection. 


Since in the present case motion is caused by the 


solute laden fluid falling from the sphere under 
the influence of gravitational forces, from con- 
siderations of continuity there must be an upflow 
some distance from the sphere, as shown in 
Fig. 11. 
convection the maximum free-stream velocity is 


It is important to note that in forced 


reached only at some distance from the sphere. 

Now in upflow the motion in the wake near the 
the the 
currents, boundary 


same direction as free 


Outside’ the 


surface is in 


convective 


1 


Transfer at low Reynolds numbers 


infinity upstream 
(no disturbance) 


2aturated solution 
neovier than solvent 


nfinity downstream 
disturbance) 
riow pottern in free convection 


veloc distribution across section O.R. 


Fic. 11. Flow Pattern in free convection. 


layer, which encloses the wake, the free and 


forced convection streams aid another 


(see Fig. 11). 
area before separation of the flow. The reverse 


one 


They also aid each other over the 


is true in the case of downflow. It is thus expected 
that 


the lowering of the mass-transfer rate until the 


interference in mass-transfer will result in 


general flow has increased to a critical magnitude. 
the 
hemisphere, this opposition is the greater in 


Since ring of separation is over the rear 


downflow. This may be regarded primarily as 
an area effect. (The actual numerical values of 
this depression are 25 per cent in upflow and 
From Figs. 9 (a) and 


(b) it is seen that there is a considerable depression 


38 per cent in downflow). 


of the mass-transfer rate over the rear hemisphere 
When the 
depression is largest at the equator. This indicates 


in downtflow. flow is reversed the 
that the direction of circulation within the wake 
has an important bearing upon the extent of 
interference. 

The 


ime 
regimes is 


the 
Figs. 


convection 
9(c) and 10. 
It is seen that the pronounced minimum mass- 


transition between two 


illustrated in 


transfer rate, which occurs near the position of 
break-away of forward-flow, becomes less marked 
at lower Reynolds number. This change is also 
demonstrated by the silhouettes of dissolved 
pellets given in Fig. 6. The dotted curves give 
the original outline of the spheres. 
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The only data known to the authors on mass- 
transfer in this region refer to packed beds. 
WINTERKAMP [24] notes that at a critical Reynolds 
number of 1-7, dye from his spherules moved 
downwards against the bulk upflow. Garrney 
and Drew [7] show that the direction of flow is 
important at low fluid speed. Mass-transfer 
rates in upflow were found to be higher than those 
in downflow and this difference appeared to a 
function of the modified Grashof number. <A 
similar difference is noted by Drypen et al. [5). 
Recently Bar-ILaxn and Resnick [2] have 
presented data on the sublimation of naphthalene 
in fixed beds and show at Np, ~ 1 the j-factor is 
a maximum, whence the Sherwood number is a 
minimum. If it is assumed that, in general, the 
overall transfer curves form a family with a 
parameter depending upon the level of the 
Grashof and Schmidt number, then this result is 
consistant with the present work. 

Drorkin and Carni [4] find that heat-transfer 
coethcients for a rotating cylinder are independent 
of the speed of rotation until a certain critical 
value. At higher speeds both free convection and 
angular speed are important until eventually only 
the latter is important. Similar three regions are 
found in the present case. 

The problem of estimating mass-transfer rates 
in the transitional region remains. It is evident 
that methods such as treating the limiting 
values i.e. transfer by free or forced convection as 
vectors (16) and defining an equivalent free- 
convective velocity to use in a Re vnolds number 
[27] represent gross over-simplifications. While 
the present data are insuflicient to draw up any 
general correlation for miass-transfer in_ this 
transitional region, the overall transfer results 
suggest that forced convection can be assumed 
to be dominant until the Sherwood number is 
equal to that in free convection. Now we have, 


[N 2 A (N,z,)* N,.* 
free convection alon 
| 2 B Na? Ng, 
forced convection alone 


Equating [Ny], to [Nx], we find that forced- 
convective effects will be solely significant when 
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the Reynolds number exceeds a limiting value 
given by: 
Nu, > 0-4 (N;,)'Ng (16) 


when A 0-6 and B= 0-95 


CONCLUSIONS 


Minimum mass-transfer rates do not occur at 
zero flow rate, since the bulk flow can sometimes 
interfere with that caused by gravitational forces, 
The depression of the mass-transfer rate is 
particularly marked over the rear hemisphere, 
especially in down-flow. 

Over the range 900 Na. 250 the overall 
transfer results are in approximate agreement 


with the theoretical relationship : 


Nig, = 0-94 Np? Ne 
NOTATION 
a sphere radius / 
a constants, defined in text 
concentration VL 
d sphere diameter L 
D diffusivity L?T! 
extraneous force VLT? 
h diminution of radius L 
A, mass-transfer coeflicient TL 
m, Me transfer rates (defined in text) 
Ny Grashof number 
(modified for mass-transfer) 
Nin Reynolds number 
Schmidt number 
p — external pressure 
radial distance L 
ro distance from axis of revolution L 
surface area L* 
time interval 
velocities in w- and y directions LT 
velocity LT 
average velocity LT 
" total mass transferred M 
zy distances mutually at right-angles L 


Z Pohlhausen parameter (defined in text) 


boundary -lave r thickness 
(momentum-transfer) L 
boundary-layer thickness 
(mass-transfer) L 
y dimensionless distance normal to surface 


(defined in text) 
0 angle from forward stagnation point 
A Pohlhausen parameter, defined in text 
shear Viscosity WL 
kinematic viscosity 
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solid density VL? 
skin-friction intensit, 


wb dimensionless concentration 


(detined in text) 


Subscripts : 


0 at surface 


Bar-Inan M 


Bosworru R, ¢ 
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Some properties of the activity coefficient 


l. SORGATO 


Istituto di Impianti Industriali Chimici, Universita di Padova, via Marzolo, Padova 
( Received 2 May 1958) 


In previous papers the author has described the isobaric maximum of the activity 


Abstract 
ular parameters and interactions. He has demonstrated the 


coefficient by means of the mok 
universal character of the generalized function 
In this paper he shows that at this singular state the compressibility factor is numerically 


equal to the activity coeflicient. This empirical correlation has been found useful in the deter- 


mination of other properties of the activity function 


Résumé |. auteur a deécrit dans des articles antérieurs le maximum isobare du_ coefficient 


dactivite par les paramctres et les interactions moléculaires. Tl a démontré le caractére universel 


de la fonction généralisé 
Dans cet article fauteur montre que dans ce cas particulier le facteur de compressibilite 
est numeériquement égal au coefficient dactivite. Cette relation empirique est utile dans la 


détermination dautres proprictes de la fonction activite. 


Zusammenfassung — In friiheren Arbeiten hat der Autor das isobare Maximum des Aktivitits- 


koeflizienten mit Hilfe von mo 
universelle Charakter der generalisierten Funktion wurde betont 


lekularen Parametern und Wechselwirkungen beschrieben. Der 


dass bei diesem besonderen Zustand der Kompressi 


In dieser Arbeit zeigt der Verfasser 
Beziehuny 


bilitatsfaktor numerisch gleich det Aktivitatskoeflizienten ist Diese empirisehe 
erweist sich als brauchbar bei der Bestimmung anderer Kigenschaften der Aktivitétsfunktion. 


Using a statistical treatment we have studied 
the activity from the molecular point of view 
and written in differential form the equation of 


THe activity function describes in a simple way 
the deviations from the ideal state. The activity 
coellicient y represents the excess of the actual 
the line of maximum activity [4]. 

This line is represented in Fig. 1 in reduced 


; co-ordinates using the molecular parameters o 


Gibbs function G over the ideal value G* at the 
same te rature 
G.* 

P P (1) and detined by the Lennard-Jones 6-12 


m* 


In y 
equation, It is a straight line and it ts surprising 


to observe that the linearity extends into the 


In previous publications [1l, 2, 3) we have dense state, in the critical regions and at high 
I g 


pointed out that when plotted as a function of the 


temperature at a constant pressure, the activity gas N,. 
In addition the experimental values represented 


in Fig. 1 are the only ones available, but are not 


pressures, for example about 2000 atm for the 


coefficient shows a maximum: 


din 4 ii* H numerous ; nevertheless the agreement is good and 
| p RT? shows the universal character of this generalized 


function. 
The state of maximum activity as any other 


This state of maximum activity is a singular 
singular state may be used as the basis of impor- 


state; it may be considered to be a poimt of 
ideality, since the difference (/7/* —H) between tant correlations, and the property of linearity 
the enthalpies is zero. should be important in this respect. 
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6 7 6 
> KT/Ew 


1. The state of maximum activity in reduced 
co-ordinates. 
Nitrogen 
A Carbon monoxide 
O Argon 
DEFINING THE STATE 
OF Maximum AcTIVITY 


THe EQUATION 


We can also obtain the equation of the line 
using the virial theory, i.e. in terms of the inter- 
these forces determine the 
values of the second virial coeflicient B(7'), the 


third virial coetlicient C’ (7) and so on. 


molecular forces ; 


Since the coefficient may be conveniently 
calculated by the compressibility factor C with 
the equation (1), we shall employ the virial 


expansion in powers of P: 


Th 
RT 


D(T) 
RT RT 


(3) 
Combining with equation (1) we obtain the relative 
contribution of each virial coefficient to the 
magnitude of the function y: 

In y B(T) C’ — 
RT 2 


1 Ps 

D(T) —. (4) 
3 v7 
Now combining equations (2) and (4) we have 


the equation of the line of maximum activity : 


| B(T)—12C 4 
T 


First we can utilize it to obtain: 


The intercept of the line of maximum activity with 
the Taxis 


Putting P — 0: 


d 
T'p B(T) _B(T) (6) 


dT 
that is the curve B(T7') as a function of 7’ must 
have only one tangent, which passes through 
the origin. In order to obtain the temperature 
for P 


3rd virial coeflicients for the gas N, as power 


0, we may express the 2nd and the 


series : 
B(T) bT cT? 
19-2.10 4; b 0-081.10~4; ¢ 
5-2.10°° Amagat 
d eT + fT? 
77.10% — 20.10%; 
1-18.10 Amagat 
Then Tp 0 608°K. This is in 
reasonable agreement with the measurement 
of the Joule~Thomson effect (7 [1] 
and with the experimental values of B (7) [5] 
600 K). The 


slight difference between these experimental 


(a/c)? 


plotted in equation (6) (7 


values is attributal to the difficulty in measur- 
ing the 2nd virial coefficient with high accuracy. 
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The slope of the line at the origin 


dP T 
(a7)... (9) 


culated for 
vd of the 


co ordinate Ss we 


13-8 


In reduced 
the gas N, the 
15-4, 

$4 of the previous paper [3 Wi 


slope 


value obtainable graphically from 


verify 


cannot 


that the slope constant into 


the region of high pressures because we 


remains 
would 
require the 4th. the 5th .... virial coellicients. 


which are not well known. 

If we neglect the 3rd coetlicient ¢ T 
fall at the 

ture given by equation (6) and deting 


B(T) T line of 


function 


then 
all the values of y, tempera 
dl by thr 
maximum value of 


would be a straight line and 


independent of the 
that the 
to produce 


QP, dT); is 


max 


shows 


pressure | l 


influence of the higher hervents ts 
sloping te 


ie vative 


TO 


ACTIVITY 


COEFFICIENT 


COMPRESSIBILITY 
Let us 


in which we can 


first consider the low 
neglect he 

virial coeflicients. Since 


In 3 B(T)P RT is 


possible to write 


rh 


RT 
and using equation 3) we obtain 
Therefore the 
with the 
littl by litth 
true for the 
the C 


sures, 


isobars ol t he function te na to 


converge isobars ol thy function ¢ 
Is also 


that also 


. as the pressure P -) 
maximum values, knowin 
isobaric maxima at low 


curves show pres 


For pressures less than 200 atm the ompress! 
bility of N, shows an isobaric maximum (Table 1); 
the increase of temperature decreases the values 
of the maxima, owing to the 
the temperature. The value of the maxima tends 
to zero, as indicated by Table 1, in the range 
300 4100 C. We remember that 


experimental 


between and 


have shown the 


and 621 K. 


previously we 
values of 600 


idealizing effect of 


Table 


Compressibility coefficients for the gas N 


Pre ssure 


(atm) 


OO417T6 


OOS476 
O45074 


O1ZS17 
O44452 


O46710 


P 0. 


Further we 


We conclude that at the limit, when 


the two maxima have the same value. 


recognize this property for any other pressure 


along the line of maxima: 


PV RT (11 


To verify it let us consider the figures of Table 2. 
calculated tables of Micuens 
We have found the points at which (//* 
O and the 


from the 


corresponding temperatures and 


Table 2. 
calculated by the 


pressures are shown in alongside the 


values of C and —_— functions 


Pil and (G* 


(11) 
cal ulated 


Table Demonstration of the equation 
figures at 250 © are 


C. The 
from the tables of PERRY (9) 


calcul 
P (atm) by eq. (5) 
Gas CO 


2485 
2-060 
1-755 


1-531 
Gas N. 


1s2 
2-604 
2-354 
1-910 
1-437 
1-181 


335-0 1-170 


1 
| 
. 
VOL 
1958/ 
( 
— 
ssure revion, 
the 4th... 
equation 
( 
0 2-471 
100 020-2 1-722 
150 753-9 1-525 
wee 
1470-7 
25 1342-3 2-654 
0 1233-1 2-347 
MO 1025-1 1-050 
685-3 1475 1402 
132 


Some properties of the 


Table 2 indicates that equation (11) is applicabl 


in the whole of maximum activity, for 
The differences he 


are of the order of 1 


range 
pressures up to 1-500 atm. 
tween and 2 per 
cent and are due to the inevitable inaccuracy in 
reading in the Tables of Michels the values of the 
the values H* 


the values of G* 


functions. In addition, while 
and (PV)* are read at P = 0, 
read at P - 
This the 
The values of y calculated by equa 


must be being G m at 
= @. 
and 


tion (3) are in good agreement with the 


Increases differences between 


max* 
experi 
mental ones at not so high a pressure. 

We were 
retical explanation for equation (11) 


not able to tind a satisfactory theo 
Anvhow, 
For 


instance, let us separate the function (¢ 1) as 


it may be useful to know some properties. 


given by equation (3) in two parts: one, which 


represents In y, the other which re presents the 
Helmholtz function A: 

RT 3 R17 

At the point ; 

function C is expressed only by the second part : 


(12) 


- the overall behaviour of the 


Ap* — A 


d 


since the first part. which refers to In y. does not 
feel at this point the influence of 7. Therefor 
the variation of on part describes the variation 


of all the As the 
eflicients are related to the intermolecular forces, 


virial coetlicients. virial co 


this point may be of interest. 


OrTrueR PROPERTIES AT THE SINGULAR 


The re lationship Se C gives further possi 
bilities of illustrating the function activity. 

As we stated, at the state of maximum activity 

V) of N 

G,*) of Gibbs function over the 

Vp*), the 


Gibbs function G. In 


every system (N, P. molecules shows 
the excess (G 
ideal system (N, P, 


the 


which at same 


pressure has samc 


addition we may apply the relationship of G. N, 


Lewis : 
G G (14) 


P = fugacity. 


RT f 


/ 


activity coefficient 


We also may refer the real system (N, P, V) 
to a second ideal system (N. P,*, l’) with the 
same Gibbs function G and the same volume V 


Now let us consider the reciprocal position of 


(instead of the same pressure 


the two ideal states. Ws easily obtain : 


pP 
G,* RT \n RT in C (15) 
activity, when 


the point of maximum 


equation (1 


wives 
(16 


that is the point of is characterized by the 


max 
equivalence of the real 
(Fig. 2). The 


G Gp*) IS equal to the energy 


state over the two ideal 


states excess of free energy. 

G,*). 
that the gas receives during the ideal compression 
state (N, P,*. V) to the state (N, P. 
In other words, the compression of ideal 
ratio ( (P P,*) gives the 
the ideal state 


owing to the appearance 


from the 
V,*). 
molecules the 
when 


magnitude of the gain, 


becomes a real state, 

of co-volume and molecular interactions. We have 

elsewhere illustrated the mechanism of this path 
2). Here we point out that we have to speak 

only of gain and not of loss of ene rev, because it 
l and hence ( l. 

mention another application of 


take the 


which refers to the con- 


Is alwavs 
/ max 
Let us now 


equation (11) and into consideration 


activity coetlicient 
centrations as stated in the statistical theory by 
MAYER Mayer [10]. We the 


situation at V = cost and define ; 


and consider 


N* /N N* 


(17) 


N* number of ideal molecules with the same 
value of the Gibbs function G, in the same volume, 
as the real molecules N. 

Substantially 
the ideality condition (N*, f, V). 
the 
lig. 2 ; taking in consideration that we are dealing 


with ideal gases, it is easy to obtain from Fig. 2: 


we refer the actual situation to 
The concen- 


trations for two ideal states are shown in 


J 
pP 


(18) 


G —G,* 
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In y, Iny + (19) 


At 


the point of maximum activity — y, 


(Yean)” 3 in all other cases y is different from C 


and the values of In y and In C could also be 


zero or negative. 


NOTATION 


B(T) CCT), DCT) 2nd, 3rd, 4th... virial co- 
efficient 


Helmholtz function 


compressibility factor 


fugacity 


Gibbs function 


= enthalpy 


number of molecules 


Conc maximum energy of inter 9 
J action 
Fic. 2. Gibbs function of a real gas (curve OM) referred to activity coeflicient 
the ideal states Point M state of maximum activity, 


molecular diameter 


We have deduced equation (18) using a statis- “7 ersertpt the ideal state 


tical treatment [4]. From this simple ratio The index for quantity in the ideal state. for instance 


between the two activity coeflicients we obtain Vp*. shows the function, in this case the pressure P 


the general correlation : which is the same for the ideal and the real state 
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Kinetics of water-gas conversion reaction 


Institute of Industrial Chemical Plants-—University of Padua. Italy 


BortTOLIN! 


( Received 28 March 1958) 


Abstract—A description is given of experimental investigations into the kinetics of the water-gas 


conversion reaction as a function of the gas composition and its specific flow-rate. The results 


obtained are compared with the calculated reaction rate, assuming that the mechanism deter- 
mining the kinetics of the process as a whole is diffusion of the gas reagents in the catalyst surface 


and of the catalysis products in the gas. It is found that, at least at temperatures of about 500 ¢ 


and higher, this diffusion mechanism actually does determine the kinetics of the entire process, 


For experimental purposes the reactions of oxidation and reduction occurring between the 


gas and the catalyst were also studied as a function of the temperature and the H,O : CO ratio. 


Useful data are derived for the design and the running of commercial reactors. 


Résumé L ‘auteur décrit son travail experimental sur la cinétique de la réaction de conversion 


du gaz a lVeau, en fonction de la composition du gaz et de sa vitesse d’écoulement specifique. 


Les résultats obtenus sont comparés avec la vitesse de réaction calculée, en supposant que le 


mecanisme déeterminant la cinétique du procédé, est la diffusion des gaz réagissant a la surface du 


catalyseur, et des produits de la catalyse dans le gaz. Ila trouvé, au moins pour les températures 


de 500°C et plus, que le mécanisme de diffusion détermine la cinétique du procedeé entier. 


Dans un but expérimental, les réactions d’oxyvdation et de réduction s produisant entre le 


gaz et le catalvseur ont aussi été étudiées en fonction de la température et du rapport H,O : CO 


L’auteur en tire des données utiles pour le calcul et le pra tique des réeacteurs industriels. 


Zusammenfassung— Der Verfasser beschreibt eine Experimentaluntersuchung iiber dic 


Reaktionskinetik der Wasser-Gas-Konvertierung als Funktion der Gaszusammensetzung und 


des spezitischen Mengenstroms. Die erhaltenen Ergebnisse werden mit der berechneten Reak 


tionsgeschwindigkeit verglichen, wobei angenommen ist, dass der Mechanismus. der die Kinetik 


des Prozesses im ganzen bestimmt, die Diffusion der reagierenden Gase an der Obertliche des 


Katalysators und der Produkte im Gas ist. Es wird gefunden, dass mindestens bei Temperaturen 


in der Gegend von 500°C und mehr der Diffusionsmechanismus tatsichli h die Kinetik des ge 


samten Vorgangs bestimmt Fir experimentelle Zwecke werden auch die Reaktionen der 


Oxydation und Reduktion zwischen dem Gas und dem Katalvsator untersucht. und zwar als Funk 


tion der Temperatur und dem Verhiltnis H,O : CO. 


Fiir Entwurf und Betrieb iiblicher Reaktoren werden brauchbare Werte abgeleitet 


INTRODUCTION variation of the type and size of the catalysts, 


IN a previous paper [1] a study was made of a and the wide variety of compositions of the 


commercial conversion unit from the thermal and reacting gases. 


thermodynamic viewpoints ; it now became con- The reaction was studied at atmospheric pres- 


venient to continue the research into this reaction sure only and the tests were limited to gases 


and to examine the kinetics and study the be- without impurities and to a single type of catalyst. 


haviour of the relevant catalyst. Work was also confined to a relatively high tem- 


A complete study of the kinetics would be long perature, mainly in order to have a sufficiently 


and difficult, in view of the very different physical high residual CO content at equilibrium for the 


conditions of pressure and temperature accom- converted gases to be analysed with sufficient 


panying this reaction in commercial practice, the accuracy. 
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These restrictions to the experiments did not, 


however, prevent us from drawing certain im 


portant conclusions. 


EXPERIMENTAL SECTION 


Method adopt d 


Firstly we will indicate the experimental method 


adopted, employing the usual formula for the 


reaction rate for reactions on a solid catalyst 


da r dW (material balance 


ada 


where G indicates the mole of gas hr. W the ¢ of 
the catalyst, x the 
(in this case the CO, or H, partial pressure 
yr the 


In order to employ equation (1), we 


concentration ol i product 
and 


reaction rate in mole hr. o of catalvst. 
should 
our experimental reactor with a very 


bed dW of 


certain conditions of flow G the very small varia 


charge 


small catalyst and measure under 


tion dz in composition. But the analytical methods 
at our disposal do not enable us to obtain sufli 
ciently accurate data from these tests 

of the 


reactor 


Equation (1) also illustrates the slop 


rvs. WG curve we can then charge the 


with a given bed W of catalvst and analvse the 
composition 2 of the converted gas under various 
conditions of G in order to draw the required 
find the 


rin ntally, instead of tinding the 


curve and hence reaction rate 

curve of 
r (CO, or H, partial pressure) directly is expe 
dient to find the curve AR, where A is the equili- 
brium constant at the test temperature and R is 
the ratio between the partial pressures, calculated 
from analysis of the converted gas. The magnituck 
of A, R is a unique function of 2: in fact, indi 
cating by 4 the initial H,O CO ratio 


that 


of the gas 


is converted 
Poo Puw 
Poo Py 
(1 A) r} {A A) 


R 


(2) 


Knowing KR is therefore equival nt to know- 


ing z, but the function A. R has certain important 
Whereas for W 


towards infinity tends towards the equilibrium 


advantages over x. G tending 
partial pressure, that is closely linked with the 
initial composition of the gas and with the tem 
perature, the function A R always tends towards 
unity, irrespective of inevitable variations in the 
rimental Moreover, the 
asymptote of AR. unlike that of 2, remains the 
same for all the various test series, it can be stated 
that a of A Rvs. WG 
curves down under different conditions is entirely 
due to the different kinetic 


reaction. 


conditions. sinec 


Variation im thy cours 


onditions of tha 


1h seriplion of apparatus 

The apparatus used is simple and very similar 
to that of Barkiry et al [2), the 
kinetics of the 

The CO, 
acid by concentrated sulphuric acid at 
60°C, is collected in a large 401. flask discharging 
The 


con 


who studi al 


reverse reaction. 


obtained = by decomposing formic 


about 


from the bottom, acting as a vasometer. 


flow of and can by close ly 
trolled 


meter, 


gas 1s constant 


the water flow is measured by a rota- 


the the CO 


saturator completely immersed in a large thermo 


From gasometer passes into a 
stat. The thermal inertia of the thermostat ke« ps 
the fluctuations in temperature of the saturator to 
within several tenths of a degree C ntigrade, and 
the height of the 
it guarantees almost complete saturation of the 
gas, so that the value of the HO CQ) ratio can 


be found simply by reading off the temperature. 


saturation column is such that 


The outlet tube passes through the side wall of 
the thermostat and then enters immediately an 
electrically-heated tubular metal container. The 
gas thus arrives in the reactor with a constant 
water content and ts alre ady heated. The reactor 
is a glass tube with an inside diameter of 24 mm. 
at the 
two sets of electrical resistances that are controlled 


with a vertical furnace centre containing 


separately. A large metal tube inside the reactor 
makes the heating uniform. The lower half of the 
reactor is filled with sandstone rings in order to 
aid in heating the gas up to the reaction tempera 
ture ; above, separated 


immediately by an 
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asbestos cap, is the bed of catalyst. The BAMAG 
catalyst (about 85 per cent Fe,O,, 15 per cent 
Cr,O;) consists of uniform particles with a mean 
diameter of about 3mm. In order to reduce the 
wall effect, asbestos fibre is placed around the 
inside of the reactor near the catalyst. 

Despite the exothermic nature of the reaction. 
the gas must not be heated beyond the predeter- 
mined temperature ; thus heating in the upper 
half of the furnace is restricted, so that the reaction 
heat is large ly drawn off from the gas by thermal 
dispersion effects. The temperature along the 
bed of the catalyst 1S kept constant to within 
about 10°C. This liability to variation in tem- 
perature is the weak point of the apparatus. 

The outlet temperature of the gas from the 
catalyst is kept at the predetermined level to 
within 83°C. The head of the reaction tube is 
kept very hot by means of electrical resistances 
in order to prevent any condensation; the gas 
thus reaches a cooling condenser or a tube con- 
taining anhydride for water-vapour analysis, and 
from this the gas is drawn off for analysis. 


The gas balances on entry and on exit agree 
a 


to within 1 per cent. 
In the tests we did not find noticeable quantities 


of CH,. 


Reduction and oxidation of the catalyst 


We shall consider briefly the possibility of 
interference between gas and catalyst. 

We know that the catalyst, marketed as 
Fe,O, (we shall ignore the other components, that 
under our test conditions cannot contribute to 
oxidation or reduction), is reduced under normal 
conditions to Fes0,. Moreover, depending on the 
oxidizing potential of the gas in equilibrium with 
the catalyst, it is possible for all the oxide forms 
and metallic iron to appear and attain stability. 
We will now study this oxidation potential in 
relation to the gas composition and determine the 


form of iron oxide in equilibrium with a gas of 


any composition. 

If we assume that the equilibrium conditions 
for the formation of H,O from H, and of CO, 
from CO are satisfied, and indicate by Kio 
and Keo, the equilibrium constants of these 


two reactions, we can write : 


Kinetics of water-gas 


conversion reaction 


H,O 


P 05 


CO 
( (2a) 


(co) 


where the chemical symbols are taken as repre- 
senting concentrations. 
Multiplying these together, we obtain : 


H,O CO, 
ro lk Kiyo A COs CO H, (9) 


In the Cus being considered, all the CO, and H, 
come from the reaction, through which their 
ratio is unity. At 500°C for example, we can 
write [3] 
log 28.308 log | 1,0) 
CO 
and draw a logarithmic diagram with a straight 
line of unit ratio that gives the oxidizing potential 
as a function of the H 20 CO ratio (Fig. 1), 
In order to e alculate the oxidizing potentials 
of the y arious tron oxides, we resort to the classical 
methods. For example, for the reaction in which 


we are mainly interested : 


3 FeO + 0-5 O, = Fe,0,(K = 1/0,*) 


2 

a 

FeO Feeto0, 


Fe FeQ 030 


2 4 


| 

04 01 a 4 1 
Fic. 1. Oxidizing potentials of gaseous reaction mixtures 
and of the various oxide forms of the catalyst. Po,, is 
given in atm and the ratio H,O,CO is moles moles. 
Temperature 500°C. 
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from Perry's [4] data we find at 500 C: from 84C (H,O,/CO 1-21) to about 60°C 
CO 0-25) and back to 84°C. 
o K 219 uP 28.438, 
log K 14-31 log Po, : From 84° to about 77° the analyses of the 
converted gas were regular; below 77° there 
was a continuously increasing surplus of CO, 


compared with Hy, indicating that the CO reduced 
is log Po. - 28-0, 1.€. considerably higher. the catalyst . 


The diagram shows clearly the limit between the 
regions of stability of Fe,O, and of FeO, and 


In addition to this value of P,, in Fig. 1 we 
shall mark the pressure corresponding to the 
reduction of FeO. The oxygen potent: il of FeO, 


above 79° again, the situation was 
reversed and the part of the catalyst reduced 
earlier was reoxidized, analysis clearly showing 


the presence of a surplus of H, compared with 
in equilibrium. CO, 


indicates how oxide mixtures can no longer occur 


For experimental! verification a test lasting A temperature of 78 C on the saturator corres- 


many hours was carried out, in which the tem ponds to a ratio H,O CO 0-76; 


this experi- 
perature of the saturator was varied in stages, 


mental value agrees very well with the value 


0 7s 12.8 1 
Fic. 2. Experimental diagrams of the ratio KR (KR is a function of the 


degree of conversion) as a function of the reciprocal of the space velocity 


W/G. The curves are based on the quantity of catalyst Wing. G is the 


flow of gas ft® hr calculated at 20°C. The reaction temperature is 500°C, 
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0-9 calculated thermodynamically and confirms 
the accuracy of the original thermodynamic data. 

It is worth while to consider some important 
points concerning the running of industrial plant. 
It is obvious that the transformation from Fe,0, 
to FeO and vice versa is harmful to the catalyst ; 
according to our experience the activity of the 
catalyst remains practically unchanged after a 
reduction to FeO and a reoxidation to Fe,0, ; 
but 


large, as it is in commercial reactors 


especially if the quantity of the catalyst is 
it loses its 
mechanical properties as a result of the chemical 
and physical transformations, and breaks easily. 

It should be pointed out that in calculating 
H,, but 
often, especially in cracking gas, there is a large 


equation (3) it is supposed that CO, 


surplus of H, that causes the stability region 
of Fe,O, to rise to far higher values of the H,O, CO 
ratio. In such cases particular care must be taken 
not to reduce the catalyst and deteriorate it. It is 
a good rule to conduct operations with a surplus 
of steam and at a rather low temperature. 
the the 
between the stability regions of FeO and Fe,O, to 


fall to lower values of the H,O CO ratio. 


Lowering temperature causes limits 


EXPERIMENTAL RESULTS AND THEIR 


ELABORATION 
Jt is clear that for our kinetic study the tests 
must be limited to the stability region of the 
Fe,0O, form and thus be confined to a gas with 
a sufliciently high H,O CO ratio; 
various test were carried out with 
having H,O CO ratios of 0-84, 1-18 and 2-04. 
Since it was seen that the degree of conversion, 


in fact the 


series yas 


and hence the reaction rate, does not depend on 
the ratio W G only, as had been expected, but 
also on the absolute value of W (or of G) at the 


same space velocity GW, the conversion condi- 


tions were examined experimentally for each of 
the compositions given above, and with W = 6, 
12 and 16-7 g. 

Figure 2 gives experimental curves showing 
K/R as a function of W/G (g catalyst,1 moist 
gas, hr) with various initial gas compositions 
(H,O,/CO ratio) and different amounts of catalyst. 
Each test 
saturator temperature, and hence the H,O/CO 


run was carried out keeping the 


ratio, constant, and thus it was possible to relate 
the volume of CO recorded on the flowmeter to 
the volume of moist gas entering into reaction. 
which is the volume referred to in the diagrams. 

All these diagrams relate to tests conducted 
with a reaction furnace at 500°C, 

From the KR vs. WG curves it is easy, using 
equation (2) to obtain the corresponding wv vs. 
W/G curves. Using equation (1) all curves of the 
reaction rate r can finally be expressed as a 
function of x; these are given as full-line curves 
in Fig. 3. 


Theoretical calculation of the reaction rate 


If the kinetics of the process were controlled 
by the chemical reaction or by the adsorption or 
desorption mechanism of any one chemical 
component in the catalyst, the rate, defined as 
the quantity of CO, produced per unit of time per 
unit of catalyst (see equation (1)), would depend 
entirely on the degree of conversion or the com- 
position of the gas, and of course on the tempera 
ture and the pressure, but should not be directly 
affected by the flow-dynamic conditions of the 
gas stream and hence by the quantity of gas 
treated in the infinitesimal reactor dW. Similarly 
the degree of conversion (or the ratio KR in 
Fig. 2) obtained in a reactor of finite dimensions 
W would depend entirely on the space velocity 
G/W and would not be affected by the absolute 
values of G and W. 

The experimental diagrams in Figs. 2 and 3 
clearly show on the contrary that in the reaction 
process as a whole, the flow-dynamic conditions 
of the gas current in the reactor are of funda- 
mental importance. It is thus logical to assume 
that the diffusion mechanism of the reagents in 
the gaseous mass in the catalyst surface and of 
the surface products in the gas current do not 
play a negligible part in the kinetics of the process. 

In order to study the correctness of this hy po- 
thesis we propose to calculate what would be 
the reaction rate if the diffusion processes alone 
controlled the process. Of the authors who have 
studied and developed methods of calculating 
the effect of diffusion on the kinetics of a process 
based on a catalyst, we should mention CoLBuRN 


[5], Curmron [6, 7], Gamson [8], Witke [9], 
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Fic. 3. Diagrams of reaction rate r, expressed in mole /hr, g of catalyst, as a 


function of the degree of conversion or partial pressure Poo, Py, The 
reaction temperature is 500°C, 


Linton [10], Hopson [11] and Leva [12, 13]. ever obviously linked with the composition on 
We shall carry out our calculation according to the catalyst surface, a composition differing from 
the method developed by Houcen and Watson — the mean value for the gas stream in view of the 
in their well-known book [14]. difficulty of diffusion of the various individual 
The chemical kinetics of the process is how- gaseous elements between the centre of the stream 
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and the catalyst surface. It is then necessary to 
calculate the difference in partial pressure of the 
various components of the mixture between the 
gas stream and the reaction surface as a function 
of the following variables: the gas velocity, its 
density, viscosity and diffusivity; the catalyst 
particle dimensions and the reaction rate, that 
influence the diffusion mechanism. 

Extending the analysis and employing the 
non-dimensional method, we obtain the following 
relationships : 


P 
AP, 


Pp M,, 
Hy, (4) 


2/3 


H 
da a. Ya F D 


(5) 


Lim 
Suitable tabulation of the available data gives : 


Yq = 1:82 (D, G/py?™ (6) 


If the rate of chemical reaction on the catalyst 
surface were such that there were practically 
equilibrium conditions on the catalyst, the rate 
of the process as a whole would be controlled 
solely by the diffusion mechanism. Knowing 
A P,, we should then be able to calculate the 
reaction rate, substituting (4) for r: 


G 
“Mo Pala,” 


r 


In order to obtain AP ,, we must calculate the 
partial pressures of the different components on 


the catalyst surface. Indicating these by a dash 
), letting indicate the pressure in 
the centre of the gas and indicating by A and 
B the reactants and by R and S the products, 


we can write: 


(7) 
where K is the equilibrium constant of the reac- 
tion. Moreover, according to (4): 

AP ,: AP,g: APp: APs 


Ha: Hay: Hay: as (8) 


1 1 1 


The four equations (7) and (8) are soluble in the 
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four unknowns AP ,, ... . and thus AP,, and 
finally by means of formulae (4'), (5) and (6) the 
reaction rate r should be calculable. 

The importance of this value of r is clear; it 
represents the maximum rate of the process, 
which cannot be exceeded at the predetermined 
temperature no matter how high the rate of the 
actual chemical reaction may be, unless of course 
the flow-dynamic system of the reactor is changed. 

The diffusivity was calculated by the method of 
GILLILAND [15], according to whom D,,,, . . 
can be calculated as the mean by weight of the 
coeflicients relating to the possible pairs of pure 
components. These binary coefficients may be 
calculated with the aid of the following semi- 


empirical relationship : 


P(v,"" vy ) My My 
where D ,, is in em? sec, T in °K, P in atm, V, 
and J, are the molecular volumes of A and B, 
M the corresponding molecular weights. 

The fact that these diffusion coefficients vary 
with the composition of the mixture somewhat 
complicates our calculations, but does not intro- 
duce any difliculty regarding concepts. 

Writing equation (7) we wish to take into 
account our hypothesis, according to which the 
gas is in equilibrium at the catalyst surface. It 
must be noted, however, that the stabilized 
equilibrium will be that relating to the temperature 
of the catalyst surface, and this is higher than the 
temperature of the gas stream (and of the control 
thermometer), in view of the exothermic nature 
of the reaction. In order to calculate this difference 
in temperature we must write : 


AT=Q/h (10) 


where Q is in cal hr, em? of catalyst surface ; 
hin cal /hr, °C, em*. Calculation of Q is easy when 
the reaction rate r (mole hr, g of catalyst) and 
the heat of reaction H, (cal/mole) are known: 


r AH, p, 
a 


(11) 


In order to calculate the heat transfer coefficient 
h, just as was done to calculate the diffusion or 


vOL. 
9 
58/59 
Pp’, 
| 
| 


mass transfer of material we can derive the 


following relationships : 


(12) 


1-05 Re 


(13) 


where Cp is in cal mole. °C: G is the volume 
specified in mole hr, Cp A ts Prandtl’s 
non-dimensional number, 

To use equations (12) and (13) presupposes a 
knowledge of the reaction rate and, hence, the 
surface temperature: we must therefore proceed 
by successive approximations, In order to lower 
the work of calculation we note that the diffusion 
the heat 


transfer 


rate) and 
the heat 


coetlicient of convection h) depend on the same 


factor yg (or the reaction 


transmission factor y, (or 
power of the Reynold’s number (cf. equations 
(6) and (13)), an algebraic consequence of the 
fact that the mechanisms of the two processes 
are identical, the temperature gradient from the 
catalyst surface and the centre of the gas being 
independent of the gas tlow rate, and hence it 
will suffice to perform the calculation for any 


value of the flow rate G. 


The characteristics of the catalyst are: 
a, Tem”? cm® (outer surface of the particles) ; 
Pr 1-3 cm* (apparent density); D, = 03em 


(mean equivalent diameter of the particles). For 
70 


are 


the gas mixture we shall suppose that Cp 
and Prandtl (Cp y,A) 0-74 
invariable with the temperature and the com- 


cal/mole, ~C 


position. For calculating the diffusivity we shall 
use the molecular volumes given by ArNovp [16]. 

The results of these calculations are given In 
Fig. 4, where the curves of r as a function of 
Po, for a constant specific flow rate G are shown 
as full lines. This means (according to our initial 
hypothesis) that, for a certain pressure Poo, and 
a certain flow rate G at a given point in the reactor 
the reaction rate at this point is r. The broken-line 
curves on the other hand are for constant W 
(g of catalyst). These curves have the same sig 
nificance as the continuous experimental curves 
in Fig. 3. In order to draw these curves for 
constant W, we have to integrate equation (1). 
Expressing G in g/hr, em*, equation (1) becomes ; 
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where S indicates the cross-section of the reaction 
tube in cm?. Substituting equation (4’) in this 
relationship we obtain an equation by means of 
which can be calculated the necessary amount of 
catalyst W for obtaining under certain conditions 
(c.g. of flow rate) a given degree of conversion &, 
or looking at it in a different way the degree of 
conversion obtainable from a bed W of catalyst 
This 


second point of v iew enables us to draw the curves 


as a function of the flow rate of gas G. 


for constant W. 


CONCLUSIONS 

The full-line curves in Fig. 3 give the expert- 
mental results of our work and the broken-line 
curves give the theoretical results, based on the 
hypothesis that the reaction rate is controlled 
by the diffusion processes between the catalyst 
surface and the surrounding gases. 

The reaction rate depends very closely on the 
flow-dynamic conditions in the reactor, hence on 
the gas 


volume per unit of reactor cross-section 
and the particle-size of the catalyst. 
These 


can be explained quantitatively if we suppose that 


differences in the rate with the flow rate 


the kinetics of the process are controlled by diffu- 
sion of the reactants and the products between 
the catalyst surface and the gases. 

The experimental curves of the reaction rate 
as a function of the degree of conversion are 
correspondingly parallel to the curves calculated 
taking the process of diffusion only into account. 

All this enables us to state with certainty that 
in the reaction process as a whole, under the 
investigation conditions, the diffusion processes 
If in addition to 
the 
reaction had had any appreciable influence, we 


are of overriding importance. 
diffusion the purely chemical kinetics of 
should not have been able to explain the position 
in the above manner. 

We still have to explain however why the 
experimental and theoretical values of the reaction 
rate differ by a factor 1-7. This may be due to 


various factors—-viz. the approximation of the 


W G S | da 
MJ r 
s/s 
VOL. 
9 
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INITIOL CONDITIONS Ry 


INITIAL CONDITIONS 


We 6.7 gr. 
W:12 gr. 
6 gr. 


a function of the degree of conversion or partial pressure Poy, 
W is the weight of catalyst in g ; 


the 
wall effect, in view of the small diameter of the 


methods of calculation ; very extensive 


reaction tube, and the relatively low thickness 
of the bed of catalyst compared with the size 


of the catalyst particles and the consequent 


heterogeneity in the bed cause a fairly irregular 


Diagrams of the theoretical reaction rate r (mole 


hr, g of catalyst) as 
Py . 
G is the flow rate in g/hr, em? 


that 
any irregularity in the tlow of gas leads to an 


distribution of gas. It is easy to foresee 

apparent reduction in the reaction rate. 
From the industrial point of view our conclusions 

lead us to aflirm that, in order to improve the 


activity of conversion catalysts, great care must 
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be given to studying their composition, their 


D4, ~ diffusivity of component A of mixture, em? ‘sec 


macroscopic structure, and he nee their specific 


; D, = equivalent diameter of grains of catalyst, em 
weight, the size and shape of the particles, th ote of 
gas, 
outer surface, etc. In addition the shape of the onail 
reactor is very important; provided that this h — heat transfer coefficient, cal /hr, em*, « 
does not unduly restrict the amount of the charge. Ha, = size of a unit of diffusion for component A, em 
long reaction tubes are always pre fe rable. MH, molar heat of reaction. cal mole 
The results of this investigation must not be 4 = diffusion factor 
extended indiscriminately to other cases. It , heat transfer factor 
can however be stated that diffusion is of over k equilibrium constant 
riding importance for temperatures of 500° and V,, ~ mean molecular weight of gas, ¢ mole 
over. At lower te mperatures the chemical reaction \P, difference in partial pressure of component A 
proper has such a retarding action that it in turn between catalyst surface and gas stream, atm 
predominates. The number of experiments made Q = heat released by the reaction, cal hr, em?® 
by us at relatively low temperatures is insuflicient y = reaction rate, mole ‘hr, gr 
to enable us to draw any definite conclusions on R — ratio between the partial pressures of the gases 
this point, but it is significant, that at a te mpera after reaction 
ture of 400° it is found that the reaction rate Y = section of reaction tube, em? 
diminishes, the fall not being du entirely to a \T — temperature difference between catalyst surface 
decrease in diffusivity or gas velo ity. and gas stream, “ 
molar concentration of CO, or H, 
NOTATION \ = thermal conductivity of gas, cal hr, em, “« 
i initial molar H,O CO ratio He Viscosity of gas, g/cm, sec 
a, catalyst surface p r unit of volume, cm*® em?*® p density of gas, g cm 
Cp — molar heat of gas, cal mole, density of catalyst, g 
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The thermal conductivity of catalyst particles 


Research and Development Laboratory, Socony Mobil Oil ¢ ompany, Ine,, Paulsboro, New Jersey 
( Received 17 July 1957) 


Abstract —Three different methods for the measurement of thermal conductivity of catalyst 


A. 


particles have been developed and are compared using silica-alumina cracking catalyst. Values 


of the thermal conductivity of a number of practical catalyst particles are reported also. 


Résumé 


Lauteur dévetoppe et compare trois méthodes différentes de mesure de la conducti- 


vite thermique de particules de catalyseurs, en utilisant un catalyseur de cracking du type 


silice-alumine. Des valeurs de la conductivité thermique d'un certain nombre de particules 


iumsilikat miteinander verglichen. 
die Wirmeleitfihigkeiten angegeben 


thermal conductivity of porous catalyst 


INTRODUC 


Tu 
particles can play an important role in reactions 
This 


occurring within the particle. situation 
arises because the heat of reaction can change 
the temperature within the particle to an extent 
determined in part by the thermal conductivity 
of the catalyst. In turn, this temperature change 
the within the 


catalyst in an further- 


may modify reaction kinetics 


undesirable manner ; 
more, highly exothermic reactions could lead to 
to such high internal temperatures that the 
catalyst itself would be damaged. 

There are two classical methods for measure- 
ment of thermal conductivity on single particles : 
(1) the steady state method in which the heat 
flux is measured through a body of well-defined 
geometry, and, (2) the temperature transient 
method in which the time course of the tempera- 
ture is measured when the boundary conditions 
on the object are suddenly changed. Because 
of the poor geometrical form and small size of the 
usual catalyst particle these methods cannot be 
applied without some modification. 

Conceptually, the steady state is the simpler 
methods. The the 


of two poor geometry of 


catalyst particles can be overcome in some cases 


de catalvseurs pratiques sont également presentees 


Zusammenfassung Messung der Wiirmeleitfihigkeit von Kontaktkérnern wurden drei 
verschiedene Methoden entwickelt und unter Verwendung eines Crack-Kontaktes aus Alumin 
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Fiir mehrere, in der Praxis verwendete Kontakte werden 


This 
made 


by grinding the particles into cylinders. 
technique, along with other modifications, 
it possible to apply the steady state method to a 
single particle of one catalyst. 

As for the temperature transient method, the 
the 


catalyst particles usually encountered for this 


transient times are too short with size of 
method to be applied easily to individual particles. 
To show this, consider an experiment in which the 
particle is heated to a temperature T, and then 
immersed in a well-stirred fluid bath at a tempera- 
ture 7. The temperature rise in the bath as a 
function of time is suflicient to determine the 
For the 


theoretical treatment of this problem see CARSLAW 


thermal conductivity of the particle. 


and Jarcer [1]. The larger the diameter and the 
smaller the thermal conductivity, the slower 
will be the temperature change. Consider a 
catalyst particle having a radius of 0-4.em and 
a thermal conductivity of 10% cal em sec °C. 
This 


equation 5 


represents an optimum case. Applying 
and the graph given in Fig. 23 on 
p. 206 of reference [1] to this data shows that 
the transient is essentially over in about 10 sec. 
Such short times make it difficult to satisfy the 
obtain 


that 


conditions 
data. It 


well-stirred necessary to 


meaningful becomes apparent 


| 
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RoBERT 


temperature transient experiments on the particle 
itself cannot be easily emploved to determine its 
conductivity. 

Both the steady state and temperature transient 
methods, however, may be used to measure the 
of 
consisting of catalyst dispersed in a 


conductivity large, heterogeneous bodies 
suitable 
matrix. The continuous phase should be a solid 
in order to prevent heat transfer by gross convec- 


tion. For example, it may consist of a castable 
plastic or a matrix of finely powdered material. 
Thus, by fabricating large pieces of cast material 
or enclosing large volumes of porous materials 


of 
possible to obtain the conductivity of the mixture 


in a vessel convenient shape, it becomes 
by conventional methods, From the conductivity 
of the mixture it is possible to deduce the con- 


ductivity of the components. 


Il. THRORETICAL 


Thermal conductivity measurements on single 


particle 


The thermal conductivity of a single cylindrical 
particle can be determined by measuring the heat 
flux through it when its ends are maintained at 
T, and le If 


heat loss through the sides is allowed, the heat 


temperatures respective ly. no 


flux, dQ dt, across the particle is then given by ; 
= KA 


i (1) 
dt Aa 


Here A is the thermal conductivity (cal sec-cm- 
C), 
A wx is the length (em), 
A is the cross-section (em*), 
A T is the temperature difference (7, 


B. Thermal conductivity using heterogeneous 
mixtures 


of 


conductivity 


The thermal conductivity mixtures will 
the of the 
components in the same way as the dielectric 


depend on thermal 
properties of a mixture depend on the dielectric 
properties of the individual components. In a 
W. PF. Jr. [2] this 


problem is formulated rigorously for dispersed 


recent paper by Brown, 


\. Sener 


particles of any shape and solved in series form. 


The solution is given as follows : 


kK 1 


5 BA) (2) 


the volume fraction of the embedded 
particles, 
the volume fraction of the 
continuous phase, 

the 


mixture, 


l 


observed conductivity of the 


A’ and 64 are detined by 


where Ky is the conductivity of the dispersed 
phase and K, is the conductivity of the con- 
tinuous phase. The quantity A appearing in the 
third and higher order term is a particle shape 
factor of the dispersed phase. 

In deriving this equation the following assump- 
tions were made : 

(a) The particle size is small compared to 
the dimensions of the heterogeneous 
body, 

K,, and K, are not too different, 
Uniform dispersion exists, and, 


(b) 
(c) 
(d) 


at the 


interface between the two phases. 


no thermal resistance exists 


These assumptions have to be borne in mind for 
experimental application. 

The first term of equation (2), which is indepen- 
dent of particle shape, gives a simple averaging 
formula. This can be improved by taking higher 
order terms into account. The second term is 
still independent of the particle geometry. The 
shape factor contributes relatively little to the 
correction even in the third term. 

Several special cases of (2) have been proposed 
in earlier literature. Maxwell derived an approxi- 
mate expression for the dielectric constant of 
mixtures on the assumption that the dispersed 
With the symbols 
used in this paper, Maxwell's equation becomes : 


particles are small spheres. 
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2K + K, 2K, + K, 


Using the definitions (3) and expanding equation 
(4) as series, one finds agreement between (2) 
and (4) up to the third order term if A is assigned 
the value zero. 

For arbitrarily shaped particles the geometry 
factor in (2) remains unknown. If it is neglected, 
equation (4) may be considered an approximation 
for randomly shaped particles also, In this paper, 
equation (4) will be used to determine A, from 
the measured values of A and K,. 

Both transient and steady state methods will 
be used to determine catalyst conductivity via a 
heterogeneous system. The transient method is 
applied to spheres of a plastic material containing 
the embedded catalyst particles. A thermocouple 
is placed in the centre of the sphere and the sphere 
0 this 


sphere begins to cool off by radiation into a 


heated to a temperature T,. At time ¢ 


medium at a temperature T — 0. The tempera- 
ture T (f) at time ¢ is given by : 


sin Ra — RacosRa 


T(t) =2T, ‘xp(—ka® t) (5) 
Ra — sin Racos 
where R is the radius of the sphere, 
it the time. 
k the thermal diffusivity. 


and aa constant. 


In order to evaluate this expression, it Is 
written : 
T(t) T(0) 2 N exp(— nt) 
The value of n is obtained from two different 
temperatures 7(t,) and T(t,) at the times ¢, 
and ¢, using : 
T(t,) 
= = exp nit t.) 
T(t.) I | 1 - 


Using this » from above, N is obtained from 
T(t,) 2 T(0) N exp (— nt) 


With known N.@ is obtained by means of a 
table of the function 


sina —2cosa 
N (5a) 


a sin @ cos x 
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which were computed by means of a digital 
computer. Finally, k is given by the equation : 


Nn k a 


For the steady state method, the continuous 
phase consists of finely powdered catalyst for 
which the conductivity is to be measured along 
with the air spaces between the particles. This 
choice of continuous phase was made so that 
condition (b) could be more easily satisfied. 
The dispersed phase then consists of large 
particles of the same material uniformly 
distributed throughout the finely divided material. 
This mixture is placed in the annular ring formed 
by two concentric cylinders of radius r, and ry. 
A radial temperature gradient is established by 
maintaining r, at a temperature 7, and r, at T,,. 
The ratio of 1, the length of the cylinder, to the 
radius r, is chosen greater than 5 so that the heat 
flow in axial direction is much less than that in 
radial direction and may be neglected. 

Assuming then that only radial heat flow 
occurs, the basic differential equation of heat 
conductions kA?T )T Mt reduces to: 


0 (6) 


d{ aT 


dr dr 


r, <r<r, when a steady state temperature 


gradient established. Upon integration, 


equation (6) vields 


r (dT dr) constant (7) 
i.e.. the heat flow in radial direction F 
2arlkK dT dr is constant within r, <r < ry. 


Let H,, be the heat per unit length produced by 
the heater. Since all the heat originates from 
the central heating rod, equation (7) can then 


he written 
il, (dT dr) 0 (8) 


Denoting the temperatures at r, and r, by 
T, and T,, respectively, it follows from (8) : 


H,, (re) 27Kk(T, T.) (9) 


If the heat is produced electrically in the heating 
rod, H, simply becomes y J?R (cal sec” ') where 
R is the resistance of the rod per unit length, 
J the current through it and y = 0-239 (cal W~'). 
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Thus, the final expression for the thermal con- 
ductivity is given by : 

J?R 
A Y In 

(cal see °C (10) 


Ill. ExPERIMENTAI 
A. Single partic 
It was found that by skilful grinding on a 
lathe, cylinders 0-4em in diameter and 0-4em 
long could be formed from the larger particles 
of a commercial silica-alumina cracking catalyst. 


These cylinders were used to determine the 


conductivity in the apparatus shown in Fig. 1. 


; 
=P, 


mai 
|| 
Freon 1) 
Catalyst 
sample 
al 
0 


Fig.l Apparatus for measuring the thermal conductivity 
of single cylindrical catalyst particles. 


The hot plate is a hollow copper rod which 
contains a tungsten heater wound on a piece of 
glass tubing. A thermocouple measures the 
temperature just below the contacting surface. 
The cold plate, also made of copper, is shaped to 


Ropert A, Sener 
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form a receptacle which holds a low boiling 
liquid; the latter maintains the temperature 
constant at this plate. The temperature is 
monitored by a thermocouple near the contacting 
plane. 

The heat flux through the cylindrical particle 
is measured by the rate of evaporation of Freon 11 
from the receptacle. This liquid was chosen 
because of its boiling point of 16°C and its low 
heat of evaporation of 43-8 cal/g. The vapour 
escapes through tube 7 (Fig. 1) and is collected 
above water in cylinder Z The volume of 
water displaced is measured by the pipette P. 
The heat flux is given by 


F = H p(dv at) 


where H is the heat of evaporation of Freon 11 
p is the vapour density 


and dV dt is the rate of evaporation of Freon 11. 


To secure good mechanical contact, the two 
plates are pressed onto the end planes of the 
cylindrical catalyst particle by means of a 
helical spring S and clamp C. The thermal 
contact between the cylinder and plates is 
improved by a thin film of mercury. In order to 
prevent heat transfer from one plate to the 
other by air conduction and convection, the 
entire device is suspended in the jar J which is 
evacuated by an oil diffusion pump at the outlet O. 
The pressure is measured by an ionization gauge. 

Since the thermal interface resistance between 
the catalyst particle and the holding plates may 
well be of the same order as the thermal resistance 
R of the bead itself, as defined by R = 1/ KA, 
the two must be separated. This can be done 
by measuring the overall resistance for pieces 
of different length and extrapolating to zero 
length. If a residual resistance is found it may be 
ascribed to the interface. 

Six different cylinders varying in length from 
0-79 to 0-43cm were measured. The applied 
temperature difference across the specimen was 
between 90 and 100°C, The evaporation rate 
due to the heat flux through the particle ranged 
from 3-6 to 7-8cm* per minute while the back- 
ground rate was 1-1 em*® min”. Background and 
effect measurements were made alternately for 
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several hours in order to level out the fluctuation. 

Figure 2 shows the results of the measurements. 
In this figure the thermal resistance R is plotted 
against the length of the specimen. According to 
equation (1) the measured points should fit a 
straight line through the origin if no interface 
resistance is present. The actual points are seen 
to scatter. The least square straight line shown 
in Fig. 2 intercepts the ordinate at a negative 


4° 7 
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Fic.2. The total thermal resistance of cylindrical catalyst 
particles of different length measured in the apparatus 
shown in Fig. 1 


value, namely 0.36. However, this value is 
of no aid in determining the presence or absence 
of contact resistance since a determination of 
its 95 per cent confidence limit by means of 


Student’s distribution vields 0-36 1-77. 
The possible sources of experimental errors 


which cause this scatter are : 

1. Additional heat transfer from the hot to 
cold plate by radiation. 

2. Additional heat transfer due to conduction 
and convection of the surrounding medium. 

3. Finite and fluctuating thermal interface 
resistance between bead and plates. 

An upper limit of the first effect can easily be 
obtained. The total energy per second radiated 


from a body of temperature 7’, and surface S$ 


into a surrounding of temperature 7, is given 
by 
E = Ze(Tf 


where o = 5- 
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surface coefficient which is about 0-3 for machined 
copper. With the actual values 7, 4x 10° °K 
and T, = 8 10° °K the above equation yields 
a total energy radiated by the hot plate 


10° cal sec™ 


Considering that only a fraction of this energy 
reaches the cold plate and that the heat flux 
through the bead is of the order of 0-5 cal see ' 
it is evident that the radiation effect is negligible. 

At the maintained pressure of 10° mm Hg 
heat transfer by convection may also be neglected. 
As is known from molecular theory, heat conduc- 
tion of gases decreases linearly with pressure from 
about 1mm Hg on downward. Thus, according 
to a rough calculation. air conduction at 
10° mm Hg is less than 1 per cent of the con- 
duction through the specimen, 


Fic. 3. The thermal conductivity of cylindrical catalyst 
particles of different length. 


Having ruled out the first two sources of error. 
the interface resistance between plates and 
particle remains to be considered. Variation in 
thermal contact for different specimens could 
account for the fluctuation of the measured 
values. In Fig. 3 the thermal conductivity vs. 
length is plotted. There, instead of fitting a 
straight horizontal line, the points are scattered 
at random within a strip defining the range of 
fluctuation. If the interface resistance is the 
only source of error, the upper limit of such a 
scattered array should represent the value 
closest to the true conductivity of the particle. 
This upper limit seems to be about 0-75 10° 
cal/sec-cm-°C, However, the statistics of six 
points is not suflicient to established this upper 
limit with satisfactory accuracy. 


< 4 
U Or! O-2 O-3 
Length /, cm 
} 
% 
- 
i > 3 
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B. Heterogeneous mixtures 
1. Transient measurements. This method was 


only applied to a silica-alumina cracking catalyst. 


Spheres 1 in. in diameter were used with a hole 


drilled in the centre to receive a thermocouple. 
These spheres consist of a lucite continuous 
phase containing from 0 to 80 per cent of 16 to 
1) mesh crushed catalyst particles. The lucite 
was used in the form of 1385 AB Transoptic 
mounting powder supplied by Buehler, Ltd., for 
preparing metallurgical specimens, 

The lucite powder was thoroughly mixed with 
crushed catalyst in the proper proportions and 
moulded into cylinders one inch in diameter at a 
temperature of 230°F under 4,000 to 6,000 p.s.i. 
The cylinders were then machined into spheres 
one inch in diameter. No lucite penetrated into 
the pores of the catalyst particles, as could be 
shown by density measurement. 

A sphere, with a thermocouple in place to 
the the 
heated to a uniform temperature of about 85°C 


measure temperature at centre, was 


in a small oven; it was then quickly lowered 


into a metal container completely surrounded 
by ice except for the hole in the top for the ball 
This can was in the form of a 
The 
time course of the temperature at the centre of 
the 
Brown recorder. The thermal conductivity of the 
the 


to pass through. 
eyvlinder 5in. high and 3} in. in diameter. 
recorded continuously on a 


sphere was 


sphere determined from theoretical 
treatment described in II B. 
Figure 4 shows the thermal diffusivity for 


was 


varying fractional volumes resulting from the 


——_ 


6) 


Diffusivity, cm sec" 


06 


Fractional volume, 


Fic. 4. 
consisting of catalyst granules and bonding lucite, as a 
function of the fractional volume x of the catalyst particles. 


The thermal diffusivity of heterogeneous spheres, 


A. Sener 


measurements. Each point is obtained from two 
or three different spheres with the same <. 
0 to x 0-8. Spheres 
of higher «x cannot be made because (1) the 


The curve ranges from « 


continuous phase serves as adhesive for the 


granule, and (2) arbitrarily shaped particles 
cannot be packed to occupy a volume to more 
than about 80 per cent with a reasonable range 
of particle size. 

Strictly speaking. it is necessary to make 
measurements for only two different values of 
0 and for 
any « > 0. But, measurement of several values 
for different 


which is especially desirable in the present case, 


the volume fraction, namely, for 
xs reduces the probable error, 


since the function N(x) given by equation (5a) 
has a small slope. 

After determining the conductivity of the 
the the 
diffusivity of the embedded catalyst granules is 


mixture and matrix phase, thermal 
then obtained by means of the Maxwell formula. 
It is found to be 3-5 « 10%em* sec! for the 
silica alumina catalyst used. With the measured 
0-19 cal and density 
the thermal conductivity is 
K, — 0-82 10* cal sec! 


C'. The error is estimated to be not larger 


specific heat of C 
of p 1-25 


then found to be 


em! 
than 10 per cent. 

2. Steady state method. Experimentally, the 
the 


eylinder is contained in a glass cylinder 32 em 


heterogeneous mixture forming hollow 


long and 2-5 cm in diameter. The heater consists 
of a uniformly and densely wound tungsten coil 
axially supported by a glass rod. To provide 
the 
heater, the coil is covered by a thin copper 
tubing having a radius of r, The 


temperature 7’, at r, is measured by a thermo- 


uniform temperature distribution along 


0-15 em. 
couple. A second thermocouple is placed at a 
point on a line parallel to the axis at r, = 2-0 em. 

In order to render axial heat losses of the 
heater negligible, the supporting glass rod is 
pencil-shaped at the end to give a point-like 
bearing. Thus, the temperature at the ends did 
not differ by more than 1-5 per cent from that 
in the middle of the heater. 
allowed from the time the inner cylinder was 
heated to the time the temperatures were taken. 


Two hours were 


150 


VOL 
9g 
1958/ 
5) i 
3} 
| d 
Zz 


The thermal conductivity of catalyst particles 


Table 


K partict 
em! 


Catalyst (cal seem 


NiUW 1-12 
Co/Mo (Dehydrogenation* catalyst) 
Cr/Al (Reforming catalyst) 70 
Co/Mo (Dehydrogenationt catalyst) 
Si/Al (Cracking catalyst) 
Pt /Al,O, (Reforming catalyst) O53 
Activated carbon O-G4 


K powder Density 
Particle Powder 
) (eal seem 3 3 
(g cnr”) (gcm™) 
1-83 1-48 
OSI 1-63 1-36 
O42 1-4 
0-33 1-54 
1-25 Os 
1-15 O-SS 
Ow O65 O52 


The resistance of the heater and the current 
through it were continuously controlled. 

The pulverized particles forming the continuous 
phase were of mesh size smaller than 80 while 
the dispersed particles were between mesh 6 and 
10. In these experiments the fractional volume 
ranged only between 50 and 68 in all experiments. 
Lower values were avoided because of sedimenta- 
tion between particles and powder which would 
disturb the homogeneity of the mixture. The 
upper limit was set by the densest possible packing 
of the particles used. In order to determine <x, 
the mixture was separated by a screen after the 
experiment. The fractional volume 8 = 1 — x 
of the powder was then measured in a graduated 
cylinder. To assure the same powder density 
in the volume measurement as in the experiment 
the containers were vibrated for a certain time 
by an electric shaker. 


A typical set of measured data for a mixture of 


Si O,-Al, O, beads with its powder is given 


below : 


Temp. of the heating rod (r, — 0-16 em), 
T, = 127°C 

Temp. at r, = 20cm, T,, — 48°C 

Electric input per cm heating rod, 


PR—20W 


Volume fraction, « 0-51 


After computing AK («) and A, from such data 


by means of equation (10) the thermal conduc- 


* 36 per cent CoO and 7-1 per cent MoO, supported on alpha-alumina, 180 m?/g 


+ 3-4 per cent CoO and 11-3 per cent MoO, supported on beta-alumina, 128 m?/g 


tivity of the embedded granules is obtained from 
equation (4). The values thus obtained for 


several different catalysts are given in Table 1. 


IV. DISCUSSION 


From the three methods described the steady 
state method with heterogeneous mixtures was 
chosen for routine measurements. In this method 
the theoretical analysis is straightforward and 
contains no assumptions which are not fulfilled 
experimentally. Furthermore, the quantities 
occurring in the final formula from which K is 
calculated are either electrical or mechanical 
quantities and thus accurately measurable. 

A comparison of the values of K for silica 
alumina cracking catalyst as measured by the 
three methods are given in Table 2. 

As is well known the thermal conductivity is a 
function of temperature. The mean temperature 
of each experiment is given in the last column 
of Table 2 and is seen to be slightly different in 
the last two experiments. CRraANDAL [3] has re- 
cently reported some new values of the conduc- 
tivity of sintered alumina for the tewperature 
range in question. His data show that the 
conductivity approximately doubles on decreas 
ing the temperature by 100°C. If one assumes 
that the alumina catalyst sample behaves in the 
same manner as the sintered alumina, this in- 
troduces approximately a 6 per cent difference 
between the results of the last two experiments. 

The value we obtained for the catalyst particles 
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Table 2. Comparison of thermal conductivity measurements obtained by the three methods. using 


a silica-alumina cracking catalyst 


Method 


Steady state single particle 
Transient heterogeneous mixture 
Steady state heterogeneous mixture 


is almost an order of magnitude lower in value 
than that obtained by CRANDAL for his alumina: 
that is, the cataly st gave 8&8 10° * cal cm see “( 


compared to 6 10 * cal/em see “C for sintered 
alumina, 
The the 


method may be significantly lower in value than 


value obtained by single particle 


the other two. However. the experiments using 
the single cylindrical particle were made under 
the the 
pores of the catalyst were filled with gas. The 


lower value obtained when no gases were present 


vacuum, while in other two methods 


suggests that the gas filling contributes appreciably 


to the conductivity. This could be due to the 


fact that the pore diameter (about 107 em) is 


small compared to the mean free path of the gas 
(about 10° 
rapid exchange of energy from wall to wall by 
The 
on the conductivity of porous solids has been 
(See 


molecules em), which allows more« 


gas molecules collisions. influence of gases 


discussed in the literature. 


Kuine 


The use of this data to determine the tempera- 


for example 


ture difference between the inside and outside 


of a catalyst particle during the course of a 


catalytic reaction 1s presented clsewhere 
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Abstract The theoretical plate concept in « hromatography has been extended to the treatment 


of the non-homogeneous chromatographic columin where the exchange coefficient k is not the 


same for all the plates in the column. The treatment is subject to the two basic assumptions of the 


theory, namely that the column is equivalent to a certain number of theoretical plates and that 


the eluent passes continuously through these plates. General ce position and elution equations 


have been derived and formulas for some special cases of practical interest were deduced \ 


numerical example illustrates the use of the derived equations, 


Résumé 


colonne de chromatographie non homogéne, ot le coetlicient d’échange & n'est pas le méme pour 


Le concept du plateau théorique en chromatographie a été étandu a l'étude d'un 


tous les plateaux de la colonne L étude est assujettie aux deux hyvpothéses fondamentales de la 


théorie que la colonne est équivalente a un certain nombre de plateaux théeoriques et que l'éluant 


passe continuellement a4 travers ces plateaux. Des équations générales de dépét et d’Alution 


ont été établies et des formules, pour certains cas spéciaux dintérét pratique, ont été déduites 


tn exemple numérique illustre Temploi des équations établies 


Zusammenfassung — Die Vorstellung des theoretischen Bodens in der ¢ hromatographie wurde 
auf die Behandlung der nicht-homogenen chromatographischen Sdéule ausgedehnt, bei welcher 


der Austauschkoeflizient 4 fiir alle Boden der Saéule nicht gleich ist. Die Behandlung geht aus 


von zwei Grundannahmen der Theorie, nimlich dass die Sdule einer bestinumten Anzahl theoreti 


scher Boden diquivalent ist und dass das Elutionsmittel kontinuierlich durch diese Boden strémt. 


Die allgemeine Abscheidung und die Elutionsyleichungen werden abgeleitet und Formeln fiir 


einige Sonderfalle von praktischen Interesse mitgeteilt. Ein numerisches Be ispiel verdeutlicht 


die Anwendung der abgeleiteten Gleichungen 


Tue differential equation which expresses the y, represents the concentration of solute on 


concentration on any plate in a homogeneous _ plate ”,k,, is the adsorption or exchange coeflicient 


chromatographic column as a function of the at plate n, N is the total number of plates in the 


column parameter » and the solution parameter column, S is the total weight of adsorbent. w 1s 


u was derived in a previous paper |] by means the weight of eluent which has crossed plate | 


of a differential material balance around plate n. and a Nw S., 


In deriving that equation it was assumed that When k,, and k,_, are constants or functions of 


the adsorption coeflicient k is the same for all 4« only as in the case of gradient elution, equation 


the plates in the column. On the other hand, (1) can be solved as a first order linear differentia! 


when the column is not homogeneous, as in the equation and the following relation is obtained 


case of a column made up of two or more layers ; 

of different adsorbents having different exchange = | k,,da 

coeflicients, then a differential material balance 

around plate gives | ks Ys exp | | dx ++ | (2) 
— Kya de dy, which is the recursion formula for plate in 


the non-homogeneous chromatographic column 


dy,, When k, and k,_, are constants equation (2 


day 


k,, Yn (1) 


reduces to 
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= exp (— 


| CXP (k,2) da 4 (3) 
The recursion formula for plate n in the homo- 


geneous column can be deduced from equation 
(3) by substituting k,, k 1 k and one gets 


y, = exp (— kr) 


| exp (Ker) da | (4) 


In this paper equation (3) will be applied to 
deduce the deposition and elution relations for 
some simple cases of the non-homogeneous 


column. 


Tue Deposition Cas} 


A solute is deposited on a column containing 
N theoretical plates from a solvent where the 
adsorption coetlicients between the solvent and 


respectively. 

Let us assume that the concentration of solute 
in the solvent entering the column is ¥, and the 
column is free of solute at the be ginning of the 
deposition process. 


By a material balance around Plate 1 


Ss 
(Yo — du dy, 


ay, 


(% — da dy, dx 


Yo — ky, 


Integrating and substituting the boundary con- 


dition y, = 0 when x = 0, one gets 


kyy 
— exp (— 
Yo 

or 


R, = 1 — exp (—k,2) 


In order to obtain the formula for Plate 2. one 
applies the recursion formula (3) 


Ye = exp (— | dx « | 
= exp (— k,2) 


A. Samp 
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exp (— kyr) 
E exp (k,r) — k, k, (k,—k,) a} 


when 0, = 0 
substituting in the preceding equation 


Yok, 
(ky — ky) ky 
and 
= 1 — —— exp (— k, j 
R, ie k, 2) 
exp (— kyr) 
ky — ky 


going to Plate 3 one can show that 


R 1 + ( k Tr) 
(k, —k,) —k) 
= kyks exp (- 
(ky k,) (ky k,) 


- ‘xp (— kar) 


and generally for Plate n 
= 1 
= exp (— k 
— k,) (ky — ky) (Kn — ky) 


> k 
(ky k,,) (ky k,,) (k,, 1 k,,) nv) 


Symbols. In order to simplify the writing of the 
above formula the following two symbols are 


introduced. 
= ky = the geometrical mean 
raised to the nth power and 
Kn (ky k,) (ky k,) (ky, k,) 
== ai (k, k,) 
exp (—k,2) exp ( —k,r) 
exp ( ~ 
(5) 
rel 
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The non-homogeneous chromatographic column 


One deduces from this formula that the order of — y k k 
: — exp (— k,x) + exp (— 
the plates is immaterial since by interchanging y° k, —k, . k, —k, 
any two plates we wind up with the same equation. , , ; 
. By applying the recursion formula again one can 


Special case. k,, varies linearly with n show that 
k, =k, +(n—I1)a a is a constant 
denoting k, — a by ky .*.k, = ky + na 
ke exp (— hk, 2) ; exp ( ky x) 
exp (— k, 2) /k, j exp | 
-(n 2)a a) 
exp (— hk. and generally for plate 
-l)a (n—2)a---- Yn 
VOL, 
k n 
9 exp ( — kyr) f 
52/59 a”! (n 1)! 


exp (ax) exp (2 ar) 


1) 
k, 


CXP (naz) 


—1) 
k 


n 

— k, 

a"-1 (n 1)! 
ex kx) /k 
exp ( rl ) mr 

(6) 
, One can also see that as in the case of deposition, 


m-l exp (rar) 


the order of the plates is immaterial. 
ELUTION oF A Zone OccurpyinGe ONE Equation (7) is similar to the Bateman equation 
Pirate ONLY derived by Bateman [3] for the concentrations 
¥, = yee of different species during a chain radioactive 
; decay process when only one species is present 


Applying recursion formula (3) at the beginning of the process. 


Ye = exp (— ky x) | | ks exp (kyr) dx + ¢ Special cases 


(1) &, varies linearly with n 


= exp (— k, x) k k, +(n —l1)a 


Substituting in equation (7) 


| | exp( = exp (ka) dx L¢ 


= exp (— k,2) 
exp (— kyr) ax2a----(n lja 


{ (ky a} 4 exp (— 


0 
i 
a a----(n 2)a 


when = 0, y, = 0 


exp (— 
(n—1l)a x (n—2)a---a 


= 
kx) 

ka 
ni Ann 
0 

‘c= 

and 
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” 
exp | ky 
1)! 
' going to plate », — 2 and applying equation (3) 
exp (— ax) again 
exp ( 
2 exp (— kyr) yf exp (— a) 
1)"-! exp } (7 lhar, 
Py ay’ PX exp da 
k 
a"— (n 1) exp (— ket) | yaya, Py du ( 
[exp (ar) (8 where 


The homogeneous column formula R é,", can hh 


be deduced as a limit to this formula as a tends 


to zero as follows: 
Integrating and substituting y,. . = 0 when 
exp (az) l ae 0 we get 
i 
a 2! 3! 
Yn, 2 exp (— kyr) af" a, 
As a tends to 0, the night hand side tends to a . 
exp (az) and for plate nm, nN, one can show that 
and lim — 
a +0 a 
mi ru 
also k,->k and k,” omg Ug) ag Pr, (9) 
k (ery where 


R,, exp ker) - 


1)! m [mt 


. 


It was shown by this author [2] that 


(2) The column is made up of tu 0” lawe rs l and 2. 


ym 


pu 


containing n, and ny plates with adsorption 


coefficients k, and ky, respectively 


where 


m 


p* ( y’ “) 


mw 
2 P“ can be evaluated directly from the Poisson 


. 
Tables as shown in Table 1 and hence re 
where 
can be obtained. 


a, Sand w* (ky kee) 2 uy “. In constructing Table 1, use was made of the 
k, k. [m 
u™ 
relation wy which was also deduced b 
Integrating and substituting the boundary con m! 
dition y,.., = O whens — 0 this author [2]. 
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Table 1. Values of YP," for different values of u 


1 [2] 3 1 2) {3 
*p 3 pe »4 vpa 


£-O00000 S-O00000 10-666667 


2-000000 1-333333 0 
LOS668 l S-O1S316 5684983 


100000 


2-000000 


135895 


2 145544 2 908422 2871790 2-812193 
3 105307 040037 3 1347997 1523793 1-289400 
‘ 142877 O7T5140 030167 366530 TS1467 742526 SATOT4A 
052653 (22487 OOTHSO OORZTSO 332022 215052 


136588 O78464 


OOLTST 


O16564 


6. (u 


1 i2 3 1 2 3 


” 


+ 400000 9-680000 14-197333 


2400000 2-S80000 2 30-4000 0 


1 1400718 1380282 ‘914718 3412277 6267723 7929610 
2 THOLSY 2 2-478575 3-789148 140462 
SOSSOS 221255 103340 3 LH63717 2-125431 2-015051 
221277 147591 OT 3664 020676 640552 1023165 1102226 O-912765 
95869 O51722 2019042 OOT7T3 5 LESS16 527917 0-384847 

016049 OOLS41 6 280088 294261 233656 O-151191 


(u 


(u 


il) \3 l i2) \3 
»>3 3 3 »5 ps & s 
I I I os I n I " I on I I 


0 1-000000 3-000000 $- 500000 500000 0 1000000 000000 12-500000 20-833333 
l 2040787 2450218 2 040787 +-006738 493262 12-34007 1 
2 12480935 1-201278 848500 2 9595ST2 3-047 166 5 4460906 6893975 
3 672125 520153 3 STOS4S 2-171818 3-274278 3-619697 
810357 209796 109560 T34074 14356844 1837454 1-782268 
5 1S4737 154620 O75176 OS4584 5 550507 ST7T337 H60097 S22 166 


BS4030 £95208 355367 


O24474 


O5S0TO2 


3 2 3) 


0 1-000000 3-G00000 6480000 7: 776000 0 1000000 5400000 14-580000 26-244000 
O72676 2627524 5923524 HO4517 10-175483 16-068517 
2 S74511 1758013 2-0990665 18256061 2 71004 345542: 6-742060 9326457 
3 1055760 1-O43903 3 2528181 5°112578 
570976 172027 BOGRS1 1-741472 2472407 2640171 
5 203562 277414 ‘195513 111318 5 “626689 1114783 1357624 1-282547 
155881 121533 073980 037338 H60915 ‘696709 ‘S85838 


. 
1. (u 2) 5. (u 1-0) 
7 
2. 2-4) 
9 
58/59 
3. @ = 3-0) 7. @ = 5) 
3-6) 8. (u 5-4) 
1 (2 
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Table 1 (contd.) 


13. (u 


6-6) 


»> 58 
Ps 


1000000 
T0413 
O20480 


521685 


Pa 


5800000 
2: 803028 
2805126 
2-065089 
1377807 


856122 


(2 
Vp oes 
2 P, 


16-820000 
12 -0O10072 
8-193357 
5208251 
3235142 
1855555 


58 


$2-518667 
20- 501695 
12-308358 
TOLOLOT 

3776065 


922417 


6 


1-000000 
OS0661 


2 P,, 6 


6600000 
5 601360 
3651667 
2756818 
069522 


1324105 


p_66 
21-780000 
16-178640 
11-566041 
7015274 
5158456 
188054 
1864739 


(3) 66 
2 P, 


47-916000 
31-737360 
20-170419 
12-255145 
7-0906689 


2-045016 


10, (u 


6) 


‘TRZO49 
958031 
714945 
554320 


6 


6-000000 
5002479 
$-O190830 
S-O8 1700 
1518060 

V63740 


(2) 
ps 
1S 000000 
12-907521 
8-O77601 
5-895802 
3-662889 
2°144820 


1181089 


23-002479 
14-024788 
8128806 
+-466007 
2-322178 


1141080 


918235 
827008 


7000000 
6000912 
5 O08 207 
3119608 
2-202600 
1593308 


24- 500000 
18-409088 
13-490881 
9-455028 
6335430 
4040830 


2447522 


57-166667 
38: 667579 
25-176698 
15-723660 
9-390230 
5349400 
2-001878 


15. (u 


8) 


2 


997971 
046382 
865771 
740823 
‘585887 


6- 200000 
5- 202029 
4°210641 
8-270259 
2 404455 
1-663665 
1077778 


19-220000 
14017971 
6-531071 
4°126583 
24620158 
1-385140 


62 


39-721333 
5- 703362 
5 9020382 
370061 
52443578 
2781460 

396320 


0000 
9004 


‘ p,® 


80000 
70003 
6-0033 
50171 
40595 
31591 
2-3503 


\2) 


82-0000 
24-0007 
18-0064 
13-9793 
6-7607 
44104 


(3) 


‘ 
P,, 


85-3333 
60-3336 
41-3372 
27-3579 
17-4381 
10-6774 

6-2670 


6-4) 


9) 


» 64 


987704 

953676 
“SS81081 
-764930 
616256 


| 


(a) 


6400000 
5 401662 
4413958 
3-460282 
2579201 
1814271 
1198015 


(3) P 6a 
” 
20-480000 
15-078338 
10-664380 
7204098 
4-624897 
2-810626 
1-612611 


\3) 
éa 
P, 


45-690667 
28-612329 
17-047949 
10-743851 
6118954 
3-308328 
1-695717 


0000 
‘9999 

‘9938 
‘9788 
“9450 


P,° 
98-0000 
80001 
700138 
6-0075 
5-0287 
4-O837 


(2) 


40-5000 
32-4999 
25-4986 
19-4911 
14-4624 
10-3787 

71793 


2 P, 


121-5000 
89-0001 
64-5015 
45-0104 
30-5480 
30-1693 
12-9900 


0 0 
1 | 
2 | 
3 
5 5 
| 14. (u 7) 
3 (2) \3 
9 
0 0 1958/5 
1 | 
2 2 
3 | 3 
4 | 4 
5 5 
11. (u = 6-2) 
(2) 
»| mm | test | 
0 1 0 l 
2 
3 
4 
5 | 5 
| | | 
| | 
12. (u= 16. (u= 
ka 
1 0 1 
| 
2 2 
3 3 | 
| 
6 6 | | 
158 


Table 1 (contd.) 


17. (u 


10) 


The non-homogeneous chromatographic column 


1-000 
1-0000 


12-0000 
11-0000 
10-0001 
98-0006 
8-0029 
70105 
6-0308 


72-0000 
61-0000 
50-9999 
41-9005 
33-9964 
26-9859 


20-9551 


j i2) 13) 
>» 10 lo p 10 p 10 
P,, « P, P,, 
0 10000 10-0000 SO-0000 166-6667 
1-0000 98-0000 41-0000 125-6667 
2 9905 80005 32-9995 02 
3 ‘9972 70033 25-9962 56-6710 
66-0136 19-0826 BO-G884 
5 50429 14-0307 21-7487 
6 +1100 10-8207 10-0190 
18. (w= 12) 
1 3 
12 12 p 12 p i2 


288-0000 
227-0000 
176-0001 
134-0008 
100-0044 
73-0185 
52-0634 


19. (u= 


14) 


P. 


1-0000 
1-0000 
T0000 
“9905 
9982 
9945 


il 
P,, 


14-0000 
13-0000 
12-0000 
11-0001 
10-0006 
9-024 
80079 


2 
2 


14 
98-0000 
85-0000 
73-0000 
61-9999 
51-9993 
42-9969 
34-9890 


372-3333 
2090-3333 
237-3334 
185-3341 
142-3372 
107-3482 


20. (u 


16) 


1-0000 
1-0000 
1-0000 
1-0000 
9986 


il) 
P,, 


16-0000 
15-0000 
14-0000 
13-0000 
12-0001 
11-0005 
10-0019 


16 
& Fe 


128-0000 
113-0000 
99-0000 
86-0000 
73-9999 
62-9994 
52-9975 


3) 


682-6667 
569-6667 
470-6667 
389-6667 
310-6668 
247-6674 
194-6699 


ILLUSTRATIVE EXAMPLE 

To illustrate the applications of these equations 
let us consider the simple case of a column made 
up of two layers ; layer 1 equivalent to 6 theore- 
tical plates and having an exchange coeff. k,, and 
layer 2 equivalent to 4 theoretical plates and 
an exchange coeff. k, such that k,/k, = 2. At 
the start of the elution process, the solute is 
adsorbed on the first plate only and its concen- 
tration is y,. It is required to determine : 


a—The concentration distribution of solute on 
the plates after an amount of eluent has passed 
through the plates such that u, = 8. 

b—The shape of the elution curve. 

c—The value of u, corresponding to maximum 


effluent concentration. 


(a) From Plate 1 to Plate 6, the concentra- 
tion of solute is determined from the equation 


R, = dy = 0-000835, R, = df = 0-002684 


0-010735, 


= 0-028626 


= 0-091604 


R, = of = 0-05725, R, = ¢5 


On layer 2; Ry = 0-091604 « k,/k, 


= 0183208 


From Plate 7 to Plate 10, the concentration is 
determined according to the equation 


R 
ny 
ni nq-1 u® 
e ay 
ke k, — ke 
Uy uy u Uy —= 4 
ky ky 
k 


[0] 


R, = Ress e~4 26 1° 


= 00183816 x 64 0-214870 
= 0-25187 


9 
58/59 
— 
R, ds R, ds 4 
— 
-_ 
0 457-3333 
1 
3 
5 
4 
0 a —— 4 a l 
1 2 
1 k, k, k, — 
» 
3 
5 
| 
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Table 2 


no 1 


me u® 
ny 1 


ay 


l 

OOO425 135335 00-0037 
OOTS41 IS 
Le O4UTST 
OST O27324 
OTS 44 OISS16 -O918 
O12277 1187 
855307 1533 

SS 5858 (MOTT 
022417 OOS02S O- 1788 
141080 1809 
O-18138 
Oo 
2-O45016 001360 1778 
2001878 120 » O-1675 
62670 3355 » O-1329 
10-0190 +461 >» 0-0619 
52-0054 6145 O-O195 
107-3482 8-317 » 00-0057 
194-6009 1-125 00-0014 


A. 3. Satp 
028 
o26 
o24 
0-22 
020 
o-12 
o10 
0-08 
0-06 VOL 
O04 9 
1958/ 
002 
000 
2 3 4 6 ? 8 9 io 
ria. 1 Concentration distribution of solute on a two-layer non-homogeneous column during elution. 
“2 
0 
2:4 24 
3 3 
36 
‘ 
‘4 
+4 54 
6 
“62 
8 
4 
Ww 10 
12 12 
14 14 
16 16 
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The non-homogeneous chromatographic column 
4) values of R,, corresponding to different values 

1! | ay of u are calculated from equation (9). Ry 

. represents the ratio between the concentration 

0-018316 64 O- 195434 on Plate 10 and the initial concentration on 
0-22905 Plate 1. It represents also the ratio between the 


rg concentration of solute in the effluent and 
, 4 
| the 


equilibrium with Plate 1 when the concentration 


concentration of solute in a_ solution in 


O-OLS316 64 0-136588 


on the plate is equal to y,. One can also plot 
0- 16008 


. Ry, vs. u. The curve will be the same as Ryo 


a 13 ' ps vs. u except that the ordinates will be multiplied 
6 


0-O18316 64 0-078464 Values of Ry, for different values of u* were 


by the constant value ky. 


0-09238 calculated according to equation (9) and are 
tabulated in Table 2. Fig. (2) is a plot of Ra. 


The different values of R, are plotted vs. n in 


Fig. (1). vs. u*. Re" ie equal to YP; and its values 


(b) In order to plot the elution curve, the were obtained from Table 1. 


Elution curve for a two-layer non-homogeneous column. 
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| 
= 
R, R, +3 26 4 
| 
VOL, 
o24 | 
0-22 | | 
020 : | 
4 
O16 | 
o-14 
R,, o12 
O10 
0-08 
0-02 
000 
2 4 6 2 16 18 
U 


Graphical solution of the equation 


For all the points 6. 


tand u, 


t4 exp ( 


(ec) Maximum concentration. In order to find dR, 


differentiates Ryo with respect to u* du* 2 
and equates to zero. 


i“ one 


max’ 


64 exp (— u*) ry - 64 exp 
R,, = exp(—u,) x wy 


[31 and for maximum R,, 


64 exp | Us) dR, 0 
du* 


Table 3 


60 


696709 999213 1-181089 1-385140 1612611 1864795 


355367 585855 1141089 1-396320 1-695717 2°045016 2: 001878 


q 
A. S. Sarp 
. 
os 
O25 
0-20 
J 
005 VOL 
9 
19052/ 
9-00 
5-2 $4 56 5-8 60 6-2 64 68 7-0 
U 
: 
k 
u : u* u* 
ky — ky Ry = — u*) | PS 
* >u® 
u*) 
since 
| 
u* 5 5+ 5-8 62 6-6 7 
— 
| Pe 2-447522 
8 
u® 
| Ps 
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non-homogeneous chromatographic co umn 


Therefore ky X ky X kgs ++ 
geometrical mean raised to the nth power 
(k, — k,) 
Plotting total number of plates in column 
{2} { 
= PY vs. u* and y = vs. u* 


. 


The two curves intersect at w=. In 
Fig. (3), u 


is found to be equal to 6-138. / Yo 


max 

total weight of adsorbent in the column 
Table 3 lists the values of | P® and | P® for hea 


(ky Ky) lbs 


2 
different values of u*. These values were obtained i weight of eluent or solvent that has passed 
3 through any plate in the column 
wNn/S 


concentration of solute on plate n, solute /y 


from Table 1 and were used in preparing Fig. 


NOTATION adsorbent during elution 


concentration of solute in eluent in equilibrium 

with plate » during elution 

concentration of solute on plate n at the beginn- 

ing of the elution process 

concentration of solute on adsorbent if in 
constants equilibrium with solvent containing the solute at 
adsorption or exchange coefficient for the homo a concentration > 
geneous column concentration of solute in solvent before entering 
adsorption or exchange coefficient at plate n for Plate 1 


the non-homogeneous column ! Poisson exponential function 


REFERENCES 
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Abstract 


was constructed at 50 mm Hy total pressure 


Nowember 


An enthalpy-concentration diagram of the binary system hydrogen peroxide 


Waltair, India 


water 


Lack of information on heat capacities of the 


solutions at different temperatures necessitated the author assuming that the excess heat capacity 


of mixing remains constant in the temperature range required in this diagram, 


vapour equilibrium data were also « 
KavaNacu and Ticknor [1 The 
Résumé — Lo auteur a prepare 


oxygenee-eau pour une pression totale ce 


un diagram enthalpi concentration du syvsteme 


mm iy 


Isobaric liquid 


tleulated at 50 mm Hy from the isothermal data of SCATCHARD 


results of these calculations are presented in a Figure 


binaire eau 


Un manque de données sur les capacités 


caloritiques des solutions a différentes température a obligé Tauteur a supposer que Texces de 


capacite calorifique du mélange reste constant dans le domaine de température de ce diagramme 


\ partir des données isothermes ce 
caleulé a SO mm ce 
ealeuls sont reproduits sur une figure 


Zusammenfassung —Fiir das biniire 


Gesamtdruck von 50 Torr ein Enthalpie-Konzentrationsdiagramm ent worfen 


SCATCHARD 


Svstem 


KAVANAGH et Ticknor I, lauteur aaussi 


Ily les donnees Is »bares equilibre vapeur Les re sultats cle cos 


Wasserstoffperoxvd-Wasser wurde bei einem 


Fehlende Kenntnis 


iiber die Warmekapazititen der Losungen bei verschiedenen Te mperaturen zwang zu der Annahme 


Ube rschuss- Wiirme kap witht der 
konstant b. 


class clic 


menwerten von ScaTcHuarRD, Kavanacu und Trckwnor (1) berechnet 


nungen sind in einem Bild wiedergegeben 


HyproGcen peroxide finds extensive use as an 


oxidizing agent in the bleaching of cellulosic 


and wool materials, as a foaming agent in the 


manufacture of porous substances, as a source of 


free radicals to initiate polymerization Processes, 
as a starting material for the preparation of the 
inorganic and organic peroxy compounds, and 
as a liquid propellant. Although only 47 per cent 
by weight of hydrogen peroxide (100 per cent) 
is available as oxygen for the combustion of a 
fuel, a substantial quantity of heat will be re 

leased when it is decomposed to water and oxygen. 
It can, therefore, be used as a monopropellant in 
which it is decomposed under pressure to yield a 
gaseous mixture of oxygen and superheated 
steam, or as a bipropellant in which the hydrogen 


peroxide is first decomposed and then the fuel, 


Mis huny 
Isobare Flissigkeit-Dampf-Gleichgewichte wurden bei 50 Torr aus den Isother 


im Temperaturbereich des Diagramme 


Die Ergebnisse dieser Rech 


such as hydrazine, burned in the hot decom- 
posit ion products. 
The this 


chemical is by the action of acids upon inorganic 


earlier commercial production of 
peroxides such as barium peroxide, This process, 
however, yields low strength hydrogen peroxide 
(6-8 per cent by weight). Over 80 per cent of the 
present world production involves the manu- 
facture, in the first place, of peroxy-disulphate 
by an electrochemical process followed by hydro- 
lysis. The concentration of hydrogen peroxide 
obtained in the above process (nearly 36 per cent 
by weight) is adequate for some of the industrial 
applications other than propellants which require, 
however, a concentration of around 90 per cent. 
Various methods of production and purification 


have been clearly outlined in the monograph 
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brought out by the American Chemical Society 
[2] and these generally involve the absorption of 
hydrogen peroxide vapours and_ fractionation 
of the liquid at low pressures. Since hydrogen 
peroxide decomposes at hight emperatures, the 
processes are usually carried out at 40-60 mm Hg 
total pressure 

The author found a great need for the enthalpy- 
concentration diagram for the hydrogen peroxide 
water system to base the design of the absorbers 
and the fractionators in the above process on a 
more fundamental and accurate basis. The 
application of these Merkel charts has been shown 
to the design of fractionators by Poncnon [4] 
and Savarir [5] and to the design of single and 
multiple effect evaporators by Hirscn 
McCaBpe [7] reviewed the use of such a chart for 
sodium hydroxide-water system in heat balance 
equations. Bosnxsakovic [8] has given a com- 
plete development of its use and applied to heat 
exchange processes and chemical processes. 
Although, the importance of using such a chart in 
solving a problem on a binary solution involving 
the first law of thermodynamics is very well 
known, not very many charts are available for 
some of the industrially important systems to the 
designers. This is primarily due to the lack of 
availability of the data on physical properties 
such as heat capacities of the mixtures at various 
temperatures and heats of mixing. 

In the following pages, a detailed account of 
the method used in preparing the Merkel chart is 


presented, 


Vapour EQUILIBRIUM 
50 mm Hg. 


LiqutIp 


SCATCHARD, Kavanacu and Ticknor [1] have 
recently reported isothermal phase equilibrium 
data, over wide temperature and composition 
ranges, obtained in a recirculation still. The data 
are believed to be the most accurate available to 
the author. Because of the uncertainity in the 
composition of the vapour, Scarcuarp et al., 
relied on the total pressure measurements and 
expressed total pressure in terms of vapour 
pressure, composition and chemical potential of 


the components as : 


Enthalpy-concentration diagram for the system hydrogen peroxide-water at 50 mm Hg total pressure 
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Py exp 


— (Bow — Vow) (Pr — 
P, exp 


— (Bin — Van) (Pr — (1) 


where 
(1 |B, B, (1 
BL 2z,,) (1 6 (2a) 


Hy, r* B, B, (3 
B,(1 — 22,)(5 — 6%,)] (2b) 
Substituting equations (2) into (1) and fitting the 


isothermal data to equation (1), the authors 


evaluated the constants as : 


B, 752 + 0.97 (3) 
B, = 85 
B, 138 
Remembering that 
RT In 5 (4 


the activity coeflicients were expressed as 


[— 1017 + 0-97 T + 85(1 — 42,) 


0-97 T + 85(3 tr.) 


u 


[— 1017 

13 (1 (9 62,,)| (5b) 
At even values of hydrogen peroxide composition, 
the bubble point temperature at 50 mm Hg was 
found out by trial and error procedure bearing in 


mind that : 
Py 50 Pa Vh P, Up, (6) 


The vapour pressures for hydrogen peroxide were 


obtained through the equation of ScarcHarp 
et al. : 


and 
VOL, 
9 
58/59 
\2 
In Ye 
RT 
+ 13(1 — 2z7,) (1 61,,) | (5a 
\ 
| 
= 


P. S. Murti 


log P, = 44.5760 - 


4025.3/T — 12-996 log T + 0-0046055 T (7) 


and for water, the equation of Kryrs [9] was 
employed. The calculated phase equilibrium data 
are presented in Table 1 together with dew point 
temperatures read from a large scale plot of 


Table 1. 


equilibrium at 50mm Hg total pressure 


Calculated values of liquid-vapour 


Mole Hydrogen | Bubbli Dew 
fraction peroxide point point 


liquid vapour (°C) 


0-000 0-000 38-11 38-11 
0-100 00-0041 ‘ 548 
0-200 0-0123 3. 60-8 
0-300 0-0280 7 65-2 
0-400 O-O0575 68-5 
0-500 00-1071 55°! 71-3 
0-600 0-1900 740 
0-700 00-3100 76-3 
800 O-4824 7 
0-900 O-7125 5 70-4 
1-000 1-000 80-42 


B. Data 1n ture Lieuip 


REGION 

The reference temperature is taken as O'C at 
which the enthalpies of the pure components, as 
liquids, are taken to be equal to zero. Integral 
heats of solution data were reported by 
GicurerRE and Moxisserre [10] at 26-9°C. How- 
ever, data at 25°C were calculated by Scuums, 
SATTERFIELD and Wenrworrn [2] using the 
excess heat capacity data and the equation (8) : 


(AH,)r, = (AH,) 


where (AH,)7,, (AH,)7, are the molal integral 
heats of solution, cal/g mole of HO, at temper- 
ture T, and T, respectively and AC, 


ACp,,(T, — T,) (8) 


om the 
difference between the heat capacities of the 
components of the solution and the solution 
itself per mole of hydrogen peroxide contained. 
The excess heat capacity data were listed in the 
A.C.S. monograph [2] and were assumed inde- 
pendent of temperature; this assumption was 


shown to hold good by the authors by a com- 
parison of the data on heat of dilution to infinite 
dilution at 14°, 20° with those reported by other 
investigators [11-14], and with those calculated 
from phase measurement of 
Scarcuarn et al. at 56°C. 

Heat capacity data for liquid hydrogen peroxide 
at different temperatures were not available in 
the literature. The data, therefore, were pre- 


equilibrium 


dicted from the ideal vapour heat capacities as 
outlined by Houcen and Warson [15]. It has 
been found that the reduced temperature, even 
at the highest bubble point temperature, is less 
than 0-55 and hence it was assumed that 
Cpt) — Cp) = constant and this constant was 
evaluated using the liquid heat capacity at 
25°C. The calculated heat capacities were found 
to follow a linear variation with temperature and 
the following equation was fitted to the data by 
the method of least squares : 


Cp, = 0-6220 


0.000235 °C. cal, g °C 


The enthalpy of a liquid mixture at any composi- 
tion, x, is then given by : 


(H; = 0.000235 dt 


. 


(0.6220 


®(Hy,o) t 2, (As), (0) 


The enthalpy of liquid water Hj. at any 
specified temperature is obtained through inter- 
polation of the tables of Keenan and Keyes [16], 
The integral heats of solution (AH’,), were 
calculated using the experimental data at 26.9°C, 
and equation (8) and then reducing to a basis of 
one pound of hydrogen peroxide. The points on 
the saturated liquid curve were similarly deter- 
mined using the data given in Table 1 and cal- 
culating the enthalpies of the liquid mixtures by 
the method outlined above. 


C. Enruatrepy Data at Dew Point 
TEMPERATURES AND IN SUPERHEATED 
Varpour REGION 


The terminal points on the saturated vapour 
curve are located by adding the latent heat of 
vaporization of the pure components to the values 
of the enthalpy at the two ends of the saturated 
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Enthalpy-concentration diagram for the system hydrogen peroxide-water at 50 mm Hg total pressure 


liquid curve. The other points on the saturated where 

vapour curve were located by the intersection of enthalpy of the superheated vapour 
the isotherm and the corresponding concentra- at a temperature ¢ °C, C.H.U Ib. 
tion of a point on the dew point — composition enthalpy of the saturated vapour at 
curve, the data of which are presented in Table 1. t, °C, C.H.U Ib. 

The heats of vaporization of hydrogen peroxide 
were obtained by using the equation (10) [2] 
which was derived from the vapour pressure 
equation and the application of Clapeyron 


[ ideal heat capacity of the vapour, 
cal/g mole °K. 


Mol. wt. Molecular weight 


piv) 


equation : The ideal heat capacity data were taken from 


AH, = 18412 + 0.021066 7? — 25.817 T (10) reference [2] and are given below. 


The data for water were taken from the tables Hydrogen peroxide 
presented by Keenan and Keyes [16]. » 
Cr, 8-514 + (0-948) 10° T — (0-3) 10° 7 
rhe enthalpies of the pure superheated vapours ae (T in °K) 
were obtained by using the following equation : 
Water 
H’ == 7-256 +- (0-2298) 10° T 


(0-0283) 10° 7° (Tin K) (12) 


Table 2. Calculated values of enthalpies in liquid region : reference enthalpies are zero for liquid 


hydrogen peroxide and water at OTC 


Liquid enthalpies 


/(t'C) 20 2! 30 At bubble point 


0-000 20-02 24-99 30-56 38-09 
(38-11°C, 
0-100 . 17-02 21-69 26-92 35-09 
(39- 2°C. 
0-200 5-38 14-10 23-46 31-96 32-62 
(40-75°C, 
0-300 : 2-97 11-29 5-5! 20-24 28-42 30-78 
(42-75°C. 
0-400 0-91 0-04 2-97 ‘ 25-31 29-63 
(45°25 
0-500 0-66 7-24 10-06 22-80 29-51 
(48° 5°C, 
0-600 1-43 9-65 21-11 30-67 
(52- 6°C. 
0-700 1-30 20-30 33-36 
(58-0 °C. 
0-800 0-064 “5: 38°24 
5°C, 
0-900 2-489 5 22-32 43-65 
1°C. 
50-57 


(80-42°C. 


1-000 


*Bubble point temperatures are given in parentheses. 
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(11) and 


integrating one can get the enthalpies in the 


Substituting these equations into 


vapour region. It was assumed that the heats of 


solution of vapours were equal to zero, and that 


the enthalpy of the vapours change very little 
with pressure, particularly at low pressures. By 
these two assumptions, it is evident that the 
mixture of linear 
the weight hydrogen 
Hence, the the 


region can be located by drawing straight lines 


enthalpy of a vapours is a 


function of fraction of 


peroxide. isotherms in vapour 


through those points on the two pure com- 


ponents. 


All the calculated values of the enthalpies are 
and 38 at 
fractions of hydrogen peroxide, and the accom- 


presented in Tables 2 even weight 
panying Figure indicates the Merkel diagram for 
this system. The phase equilibrium compositions 
are also shown as dotted lines and the inscribed 


figures indicate the bubble point temperatures at 


50 mm Hg total pressure, read from a large scale 
plot of t r 
Table 3. 


saturated 


Enthalpy data of the 
and 


vapour in the 


SU pe rheated region : reference 
enthalpies are zero for liquid hydrogen peroxide and 


water at O 


Saturated Region Superheated Region 


Enthalpy 
(¢ AL Ib) 


Temp. 
(©) HO, 


506 


120 
421-4 
160 128-4 


576 
554 
534 
4558 
$45 
$255 


wl 


NOTATION 


constants in the equation 
heat 


cal g mole 


capacity of the saturated liquid, 
heat capacity of the liquid, cal/g ¢ 
heat capacity of the vapour in the ideal 
gaseous state, cal/g mole K 
heat 
mole of hydrogen peroxide, cal/g mole 
the 
Ib re spectively 

the 


Ib respectively 


EXCESS capacity of the mixture per 


enthalpy of Ib mol 

enthalpy of 
and 

integral heat of solution 
H,O, and C.H.LU, 
vapour pressure in mm Hg of 


liquid 
and ¢ 
vapour, ib mok 


C.H.U. 
Ib Hf , respective ly 


hvdrogen 
peroxide and water respectively 

total pressure in mm Hy 

gas constant 1-087 cal/g mole hk 
temperature 

temperature K. 


molal volume of the liquid, cm®/g mole 


mole fraction and weight fraction of 


hydrogen peroxide respectively 
excess chemical potential, cal g mole 
activity coefficient 


8 second Virial coeflicient cm’ /g mole 


Subseripts #, w refer to hydrogen peroxide and water 
respectively. A prime indicates either weight fraction or 


per pound of the substance 
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Enthalpy-concentration diagram for the system hydrogen peroxide-water at 50 mm Hg total pressure 
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Abstract As incidental to our main studies on mass transfer, flooding rate data have been 


obtained in a packed liquid-liquid extraction tower with no solute transfer, with 2 in. Raschig 


rings, } in. copper rings and 6 mm glass beads as the different packing materials, using water as 


the continuous phase and each of the following solvents, toluene, Pegasol, kerosene and methvl 


isobutyl ketone, which do not preferentially wet the packing, respectively as the dispersed phase 


Visual observations, aided by photographic studies supplemented the measurements of the 


flooding velocities. The flow rates at the flooding point have been satisfactorily fitted on the 


square root plot as a preliminary correlation, and have been well correlated by the method of 


Dent and Prarr. It was observed that flooding rates decreased with decreasing size of the packing 


for any given system, and that higher flooding rates are favoured by high density difference and 


low interfacial tension of the liquid systems 


Résumé— A loceasion de leurs principales études sur le transfert massique, les auteurs ont 


obtenu des résultats expérimentaux sur lengorgement produit dans une colonne a garnissage 


pour Pextraction liquide-liquide, sans transfert de solute. Les garnissages suivants ont été utilisés 


anneaux Raschig de 0.95 om, anneaux de cuivre de 0-64 em et billes de verre de 6 mm. La phase 


continue était de eau, alors que la phase dispersée contenait un des solvants suivants: toluéne. 


Pegasol, keroséne, meéthyléthyleetone, qui ne mouillent pas préférentiellement le matériau de 


remplissage. Les mesures de vitesses dengorgement ont été réalivées de maniére visuelle ainsi 


que par des études photographiques 


Une corrélation préliminaire a permis aux auteurs de reproduire d'une maniére satisfaisante 


les débits a Tengorgement par la méthode des racines carrées. Loutilisation de la méthode de 


Deut et Prarr a permis de corréler les résultats avec suceés, Les auteurs ont observe que les vitesses 


engorgement diminuent avec les dimensions du garnissage pour tout systéme, et que les grandes 


Vitesses Cengorgement sont favorisces par une grande différence de densité et une faible tension 


interfaciale des svstemes 


Zusammenfassung -Im Zuge unserer Untersuchungen tiber den Stoffaustausch wurden 
Flutungsgeschwindigkeiten in einer FallkOrpersiiule zur Flissig-Flissig-Extraktion ohne Cber 


tragung von Gelostem ermittelt Als Packungsmaterial dienten Raschigringe von 9-5 mm Kup 


ferringe von 64mm und Glasperlen von 6mm. Die geschlossene Phase war Wasser und die 


disperse Phase waren folgende Lisungsmittel, Toluol, Pegasol, Kerosen und Methvyl-Jsobutvl 


Keton, welche das Packungsmaterial vorzugsweise nicht benetzten. Visuelle Beobachtungen 


unterstitzt durch photographische Aufnahmen, ergiinzten die Messungen der Flutungsge 


schwindigkeiten. Der Mengenstrom am Flutungspunkt wurde befriedigend durch eine quad 


ratische Beziehung in vorliutiger Form wiedergegeben und konnte gut nach der Methode von 


Dew. und Prarr dargestellt werden. Es wurde festgestellt, dass die Flutungsgeschwindigkeiten 


abfielen mit fallender Packungsyrosse fiir ein gegebenes System und dass hOhere Flutungsyeschwin 


digkeiten begiinstigt werden durch hohe Dichtedifferenz und niedrige Grenztfliche nspannuny 


der flissigen Systeme 


INTRODUCTION and Exon [5] described their investigations of 
Mvucu of the early work on flooding rates in packed the design and operation of spray and packed 


liquid-liquid columns was done incidentally to columns, pointing out the importance of the 


other investigations, and in this category came proper design of the entrance sections and 
the reports of Rusuron [1], Apre. and [2], distributor. 
Suerwoop et al, [3], and Row et al. [4]. BLANpING During recent years, flooding rates have been 
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studied in packed liquid extraction towers by 
several investigators [6-11], and of these special 
mention should be made of and Prarr, [9] 
who have made a rather extensive study of flooding 
rates in columns of 3in. and 6 in. inside dia., 
using a wide variety of liquid systems and 
The their 


data by means of a theoretical equation : 


yackings. above authors correlated 
£ 


0-885 ‘ Un) 


PC 


The value of n was found to be 1/4. The 
constant C has slightly different values for each 
type of packing used, and the values recommended 
for design are 0-68, 0-80, and 0-88 for Raschig 
rings, Lessing rings and Berl saddles respectively. 
Each of the flooding correlations so far developed 
for packed liquid-liquid columns is applicable 
over a particular range of the different operating 
variables, and as pointed by Treypat [12], it is 
necessary to exercise due caution in utilizing any 
of the flooding correlations to cases where the 


packing sizes, system properties and also column 


Table 1. 


sizes are outside the range for which they are 
applicable. 
APPARATUS 

The extraction equipment consisted essentially 
of an extraction tower and the necessary acces- 
sories to maintain a steady flow of the two liquid 
phases. The design of the packed column was 
based on the recommendations of BLANDING and 
ELGIN [5], and the apparatus used has been fully 
described elsewhere [13,14]. To facilitate visual 
the 
column proper was constructed of Pyrex glass 
1-88 in. 


observation of column in operation, the 


tubing, inside dia., and two column 


lengths of 60 and 34 in. respectively were used in 
the experimental work. A dispersed phase dis- 


5 


distributor with 21 °; in. nozzles was used in 


the flooding rate studies. 


MATERIALS 


Of the dispersed solvents, kerosene and methyl! 
isobutyl ketone were supplied by the Burmah 
Shell Oil Co., Pegasol by the Standard Vacuum 
Oil Co., and toluene by the Bengal Chemical and 
Pharmaceutical Works, Calcutta. The physical 
solvents are recorded in 


properties of these 


Table 1. Tap water of negligible acidity, drawn 


(a) Properties of packing materials 


Average size of units 


1D. 


(in.) 


Packing Material 


in. Raschig 
rings 

} in, Copper 
rings 


Porcelain O-27% 


Copper 0-175 


6mm Glass beads Glass 


6 mm dia, sphere 


No. of 
Length units 
(in.) Calculated 


per ft 3 Measured 


0-384 25,680 0-610 0-667 


0-635 
0-381 


215-4 
198-0 


0-692 
0-396 


81,060 
163,000 


0-250 


(6) Properties of solvents at 30°C 


Solvent Density at 30°C (g/em*®) 
O-778 
O-785 
0-859 
0-800 


Kerosene 

Pegasol 

Toluene 

Methy! isobutyl ketone 


Viscosity in c.p. Interfacial tension (dyn/cm) 
1-08 

0-854 

0-556 

0-520 
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from the Andhra University mains supply was 
used in the experimental work. 

Three types of packing materials were used : 
(1) 2 in. 
Raschig rings, manufactured and supplied by 
the Maurice A. Knight Co. of Akron. (2) Lin, 
copper rings, and (3) 6mm glass beads. The 


non-porous and unglazed porcelain 


column was filled with water, and the packing 
material dropped in from the top, a few pieces at 
a time, slowly and at random to the desired 
height. The packing was not shaken or tamped. 
The properties of the packing materials are listed 
in Table 1. 

With 


materials, the 


different 
studied 


each of the three packing 


liquid systems were : 


1. toluene—water, 2. Pegasol—water, 3. kerosene 
water, and 4, methyl isobutyl ketone- water, using 
water as the continuous phase and the organic 


solvent as the dispersed phase. 


PROCEDURE 


Kach of the liquid phases was first saturated 
with the other, as this is particularly necessary 
when the solvent and water have appreciable 
mutual solubility, as in the case of water and 
methyl tsobutyl ketone in the present investiga- 
tion. The 
admitted into the column and then set at the 


continuous water phase was first 


desired flow rate. The dispersed solvent phase 
was then slowly introduced, and its flow rate was 
gradually increased in 5 to 10 steps until the 
flooding point, as indicated by the first appearance 
of a thin layer of the dispersed phase underneath 
the packing support, was reached. In a few runs, 
when the dispersed phase rate was high, the 
procedure was reversed by first establishing the 
rate of the dispersed phase and then gradually 
increasing the continuous phase rate up to the 
flooding point. The flooding results obtained by 
the two methods agreed quite well. 

In all the runs extra care was taken to maintain 
the liquid-liquid interface in the upper end 
section at the same level, and during each run, 
the flow rates were not only checked, but also 
the exact rates of the solvent and water phases 
were obtained by timing the exit streams for a 
time and measuring their 


known interval of 


volumes. 
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CoLUMN BEHAVIOUR 


A good number of photographs were taken of 
the column during operation with a view to 
interpreting the flooding data in terms of drop 
size, holdup and coalescence of the dispersed 
phase. In general, it was observed with each of 
the packings and liquid systems used that the 
streams of the dispersed phase dropk ts worked 
their way up the column, deforming and following 
a torturous path while passing up through the 
interstices of the packing. With increased flow 
rate of the dispersed phase its holdup was 
markedly increased, and the drops became very 
closely packed up in the column. As the flooding 
rate was closely approached the drops were 
observed to coalesce, forming slugs which gradu- 
ally tended to travel down the column. The 
flooding rate was taken at the first appearance 
of a thin layer of the dispersed phase below 
the packing support. As compared with # in. 
Raschig rings, the column was rather densely 


packed with 4 in. 


i rings, and more particularly 


with 6mm beads. With the latter packing the 
drops showed a greater tendency to coalesec nce, 
and for any liquid system the flooding rates were 
lower than those for ring packings. 

It was observed that smaller and finer drops 
were formed with solvents like methyl isobutyl 
ketone, and comparatively larger drops, showing 
increased tendency to coalesce, with kerosene and 
Pegasol. No attempt was made to determine 
either pressure drop or holdup at the flooding 
point. 


RESULTS AND CORRELATION 


Altogether over 50 runs have been taken on 
flooding rates with the different liquid systems 
and packings, and the data are presented in 
Table 2. Maximum flooding rates were obtained 
with methyl 


2 with the 
isobutyl ketone-water system, and the least 


Raschig rings 
with kerosene-water and Pegasol—water systems. 
The larger drops, and their increased tendency to 
coalesce observed with the last two systems, as 
compared with the relatively smaller drops of the 
ketone—water system, may perhaps be due to the 
fact that the interfacial tension (c) of the ketone 

water system is only 10 dyn/cm, which is much 
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Table 2. Experimental flooding data: 1-88 in. i.d. column, 21 °, in. nozzle distributor. 


Continuous phase : water. Dispersed phase : solvent 


Series No. Run No. System 


M.1.B.K. water 


Kerosene water 


Toluene water 


Toluene water 


Kerosene water 


Pegasol water 


water 


Kerosene water 


Pegasol- water 


Toluene-water 


Packing Vaterial (ft /hr) Up (ft/hr) 


2 in. Raschig rings 19-0 
64-0 
103-4 


, in Copper rines 


G mm Glass beads 


wea 


smaller than those of kerosene—-water 
dyn/em) and Pegasol-water 29 dyn cm). 
In general, flooding rates decreased with de- 
creasing size of the packing. These observations 
are in quite good agreement with those of the 
earlier investigators. The flooding rate data 
(Table 2) show clearly that higher flooding rates 


are favoured by high density difference and low 
interfacial tension of the system. 

\ general type of plot which has been utilized 
by previous investigators in preliminary cor 
relation of flooding data is to plot the square root 
of flow rate of the dispersed phase against that of 


the continuous phase at flooding. On such a plot 


: 
I 
2 
3 
27°5 
5 11-6 143-0 
105-6 20-4 
I! 7 64-9 38-5 
11 18-2 115°5 
12 30-3 SS-2 
ltl 13 : 
15 29-7 17:1 
7 59-4 
9 I\ Is 16-3 
55 32-5 13-0 
38/59 19-3 56-2 
21 T60 i 
14-6 30-3 
25 10-6 26-5 
\ 25 23-7 21-5 
12-7 33-6 
27 20-3 30-9 
28 36-4 12-1 
VI 26-5 71-6 
19-6 46-3 ai 
ol 15-4 84:3 
32 34-2 
33 79-3 
17-4 56-2 
36 10-8 28-1 
37 22-6 16-0 
4-9 190 
110) 21 
16-5 35 
25-9 26 
15-4 
26°5 29 
39-1 16-0 
7 20-8 23:1 
22-6 28-7 | 
25-9 12-7 
16-5 
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_raschig | copper mm gicss 
rings rings beods 


+ 


MLBK -—woter « 


Kerosene 


Pegosol- water | 


2 
Vea 
ge? 


Fic. 1. 


a series of almost parallel straight lines was 
The 
slope of all these lines was found to be close to — 1, 
1-0 to 1-1. An attempt has 
been made to fit the data by the well-established 


obtained. one for each set of conditions. 
the range being 


correlation of Dei. and Prart, and such a plot is 
represented in Fig. 1. The flooding data are quite 
satisfactorily fitted by the dimensionless equation 
and the computed values of the slope, », and the 
and 0.66. 


be 


found to 


constant C are 
These values are in quite good agreement with 
the 4 and 0-68 obtained by the 
above authors for Raschig rings. The 
fit of the flooding data on De. and Prart’s plot 
the different 
constant C in the 


equation is almost identical and constant at 0-66, 


values of - 
excellent 
with 


of correlation shows that, 


packings used, the packing 


probably because of the closer size range of the 
packings. In view of the excellent correlation of 
the flooding rate data by the method of Dex. and 
Pratt no other correlation has been attempted. 


CONCLUSIONS 


Flooding rate data have been obtained in the 
packed extraction tower, as incidental to our 
main studies on mass transfer, with @ in. Raschig 
rings, } in copper rings and 6 mm glass beads, 
using water as the continuous phase and each of 
the solvents toluene, and 


Pegasol, kerosene 


~ 200 


~woter 


oluene — water x I 


Sp 


Flooding correlation of Deu. and Prari 


methyl isobutyl ketone which do not preferenti- 
ally wet the packing, as the dispersed phase. 
Visual 
recordings of the column in operation supple- 
the 
Maximum 
in. 


methyl isobutyl ketone-water system, and the 


observations aided by photographic 


mented measurement of flooding rates. 


flooding rates were obtained with 


Raschig rings as packing material and 


least rates with kerosene-water and Pegasol- 
water systems. 

The flooding data have been fitted satisfactorily 
on the square root plot as a preliminary cor- 
relation, and have been well correlated by the 
method of Dett and Prarr with a value of — } 
for n and 0-66 for C. It was observed that flooding 
velocities decreased with decreasing size of the 
packing for any given system, and that low 
interfacial tension of the system and high density 
difference tended to give higher flooding rates by 
favouring the formation of smaller drops. 


NOTATION 


interfacial contact area per unit tower 
volume ft? /ft®. 
constants in Det. and Prart’s equation, 
acceleration due to gravity. 
flow rate of continuous water phase in 
ft® /(ft*) (hr). 

= flow rate of dispersed solvent phase in 
ft® /(ft*) (hr.) 


10 
| | ae 
£8). 
0 20 30 60 MM 400 600 1000 3000 6000 
= 9 
1958/ 
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void fraction of the packing material, Ap = density difference of the phases in lb ft.3 
ft? viscosity of the continuous phase _ in 
density of the continuous phase in Ib /ft*. Ib /(ft) (hr). 

density of the dispersed phase in Ib /ft*. oor y = interfacial tension in dyn/cm. 
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Abstract 
transfer equipment, handling nor 


iturated 


A new graphical method is presented for the design of simultaneous heat and mass 


is-Vapour mixtures 


The condition of the non-saturated gas-vapour mixture is followed from point to point ina 


semi-logarit hmix 


temperature. For this purpose 


d log 
dt 


plot of the partial pressure of the non-diffusing gas component against the 
the following equation has been derived, 


(>) 


The application of this equation ts shown for a simplitied cooler-condenser example, in which 


the right hand side bracket groups of the equation were taken equal to unity. 


Résumé 


L/auteur présente une nouvelle méthode graphique pour le caleul d'un appareillage 


de transtert simultané de masse et de chaleur, opérant sur des mélanges non saturés gaz-vapeur. 


L’état du melange non sature gaz-vapeur est suivi point par point sur un graphique semi- 


logarithmique de la pression partielle du composant gazeux non-diffusant, par rapport a la 


temperature 


(“ \ log p, 
it AM 


Dans ce but Pauteur a posé léquation suivante 


(ae —1 
out 


Il montre lapplication de cette Cqguation pour un exemple de condenseur simplifié, pour lequel 


les groupes entre parenthéses de la 2éme partie de léquation sont pris égaux a |. 


Zusammenfassung Fiir die Berechnung 


von 


Apparaten mit gleichzeitiger Wirme- und 


Stoffibertragung wird eine neue graphische Methode mitgeteilt, die nicht-gesittigte Gas-Dampf- 


Gemische behandelt. 


Die Bedingung des nicht-gesittigten Gas-Damp-Gemisches wird punktweise in einem halb- 


logarithmischen Diagramm verfolgt 
ponente tiber der 


abgeleitet worden 


d log p, A log p, 
| dt Ml 


Temperatur aufgetragen ist. 


in dem der Teildruck der nicht diffundierenden Gaskom- 


Fiir diesen Zweck ist folgende Gleichung 


(<—"). (4) 


/h 


Die Anwendung dieser Gleichung wird fiir ein vereinfachtes Beispiel eines Kiihler-Kond- 


ensators gezcigt, in welchem die in 


Gleichung gleich eins gesetzt werden 


INTRODUCTION 
Recent y, the author [1] proposed to use a dia- 


gas 


gram of the logarithm of the non—diffusing 
against the temperature for the determination of 
the gas-liquid interface conditions in the design 
of cooler condensers. The graphical design method 
was later improved [2], and it was also shown that 


the method is suitable for the calculation of 


Klammern stehenden 


Ausdriicke auf der rechten Seite der 


equipment in which gas and liquid are in direc 
contact, e.g. gas saturators, gas coolers, water 
heaters, dryers, ete. 

The use of a diagram of the above type enables 
the elimination of all trial and error calculations 
in the determination of the gas liquid interface 
conditions. 

When a non-saturated gas-vapour mixture is 
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subjected to simultaneous and mass-trans- 
fer, it will usually be necessary to follow the 
condition of the mixture from step to step. This 
may be done by using the equation proposed by 


CoLBuRN [3], 

dt!, Al Pos Se 

Equation (1) can be derived from the equations 
for the condensation of binary vapour mixtures, 
Drew [4]. 
derivation of an equation of this type was given 


given by CoLBuRN and Another 
by the author [5]. In the same article, equation 
(1) was used as the basis for a graphical design 
procedure for cooler condensers. The application 
of equation (1) for the design of a direct contact 
cooler for hot gases was published later [6]. 

These design calculations were carried out by a 
graphical point to point procedure in a diagram 
of partial vapour pressure against temperature. 

Thus, in the general case of a non-saturated 
gas-vapour mixture, subjected to simultaneous 
heat and mass transfer, the design computations 
would involve the following two steps: 


(a) The graphical determination of the inter- 
face conditions in a diagram of log p, 
against 

The graphical step by step determination 
of the condition of the gas vapour mixture 
in a diagram of p against ¢. 


Consequently, two different diagrams would be 
required, 

In order to avoid this complicated procedure, 
the following method was developed, which 
allows to carry out both operations in the same 
diagram, namely log p, vs. t. 


DIFFERENTIAL EQUATION 


and mass transfer can be 


The 


written, respectively, 


h, AtdA 


rate of heat 


- Gas dt, 
and, 


k,' {In(P — p,) —In(P — p,)}da 


From equations (2) and (3) follows, 


A ln dH 
(4) 
At dt, 


The ratio of the gas side mass, and heat transfer 
coefficients, k,'/h, may be derived from any of 


Here, 


the correlation proposed by and 


the correlations available in the literature. 


BURN [7] will be used, or, 
k,,’ M, 


h, cM 


m 


(5) 


In case of packed towers, somewhat different 
correlations of heat and mass transfer are usually 
preferred. 

The humid heat capacity, s, is equal to, 


e=c(l i) (6) 
The average molecular weight, M,, can be ex- 
pressed by, 


P.) M, 
P 


M Po M, 


The humidity, H, is given by the equation, 


p, M. 

H 

(P p,) M, p, M, 
(P — p,) M, 


From equations (7) and (9) follows, 


H)(P 
m P 


} M, 


M (10) 


When equations (5) and (6) are substituted in 
equation (4), equation (11) is obtained, 


Al 1+ H)M 
np, ( Jn (11) 
Al M,, Se Ih dt, 

Substitution of equation (10) in equation (11) 

delivers 


12 
dt 


Alnp, P 
“At (P — p,) (se, 


Differentiation of equation (8) gives the relation, 


d 
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Substitution of equation (13) in equation (12), 
gives after rearranging, 


dp, Aln p, 
(P Pr) dt, At Sc 


Equation (14) can now be written, 


d\n Pa A In Pe. (ery (15) 


dt, Al Sc Sh 


If the correction for the sensible heat transferred 
by the diffusing vapours is taken into account, 
the rate of sensible heat transfer at the gas side 
of the gas film can be expressed by an equation, 
recently derived by the author [6], or, 

(16) 


by sdt, h, At - 


— 3) 


Introduction of the correction factor in the above 


derivation will give the relation, 
e l 

(17) 

Dh ec 

to the ratio of the heat transfer coetlicient for 


2 


din py Aln Pe. (<") 
dt, At Se 


In this equation, « is equal 


the sensible heat transferred by the diffusing 
vapours and the sensible heat transfer coctflicient, 
h,. Before the numerical value of ¢ can be caleu 
lated, the 
determined. 


interface conditions must first be 
The determination of the interface 
the 


sensible heat transferred by the diffusing vapours, 


conditions, when making allowance for 


was discussed in some more detail in a previous 


publication, to which reference may be made [2]. 


Once the interface conditions have been estab- 
lished, the correction factor, (e* 1) (e e*), can 
be calculated, and equation (17) can be used in a 
straightforward manner and no further trial and 
error work is involved. It should be borne in 
mind, that the correction factor for the sensible 
heat of the diffusing vapours, (e* 1) /(e e*), 
has been theoretically derived, but needs more 


experimental verification. 


GRAPHICAL PROCEDUR! 


For the point to point calculation of the relation 
between the gas temperature, ¢,, and the water 


temperature, ¢,, reference will be made to Fig. 1. 


Partial pressure, 


Temperature 


i}. 1. Graphical Determination of Change in Gas Phase 


The curved line represents the partial pressure 
of the 
temperature in a 


saturated non-diffusing gas against the 


semi-logarithmic diagram. 
Point A, represents the known condition of the 
inlet gas-vapour mixture, p, 
Point C, the 
Pe = Poy. and t = ty. A method for the graphical 
the 


been described elsewhere [1, 2). 


Poy andt=t 


represents interface conditions, 


determination of interface conditions has 
The slope of the straight line A, C, will now 
A log p,/ Al. Application of 
equation (15) will further enable the calculation 
of d log Per dt,. 


equal to d log Pe dt, is drawn through point A,, 


be equal to 
Another straight line with a slope 
see Fig. 1. The direction of this line will indicate 
the initial trend in the change of the non-saturated 
gas vapour mixture, in respect to the logarithm 
of the partial pressure of the non-diffusing gas 
and the temperature. Although this trend is 
strictly valid for point A, only, the line may be 
followed over a relatively short distance, e.g., to 
point A,. The co-ordinates of point A,, p 
and t, 


£ 
g» can now be read from the diagram. 
Then, the temperature and the partial vapour 
pressure at point 4, are known, and the enthalpy 
of the gas mixture can be calculated. The differ- 
ence in enthalpy between the points A, and A, 
is equal to the total heat transferred. 

In case of direct contact operations, all this 
heat is exchanged with the liquid. 
cooler condensers, some of the heat will be in the 


In case of 


condensate. This amount can usually be neglect- 
ed, especially if the latent heat of the vapour is 
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relatively large. Thus, the change in the liquid 
temperature ¢,, can be calculated from the en- 
thalpy change, by application of an equation of 
the type, L ce, dt, = G, di,, after integration. 

After calculation of the liquid temperature, 
tig, in this way, the interface conditions for point 
A, are again obtained graphically. The procedure 
is repeated step by step until the outlet of the 
apparatus is reached. The typical curved line 
for the condition of a non-saturated gas-vapour 
mixture in a cooler condenser is shown in Fig. 1, 
by the curve A, — Ag. 

Sometimes, the gas may still be non—saturated 
when leaving the apparatus. Then, the exact 
partial vapour pressure and temperature of the 
exit gas will not be known until the above step 
by step procedure has been completed. 

The enthalpy of the exit gas vapour mixture 
will be equal to the enthalpy originally assumed 
for the calculation of the ratio L Gy. However, 
the temperature and partial vapour pressure of 
the gas vapour mixture may be different from the 
conditions originally assumed, Obviously, a trial 
and error procedure for the entire design calcula- 
tions would be required to let the gas leave at 
the temperature originally fixed. In most cases 
it may prove more practical to design the equip- 
ment for a certain gas exit enthalpy, rather than 
for exact gas exit temperature and partial vapour 
pressure, in order to avoid this complication. 

In direct contact equipment, the liquid flow 
rate, L, may vary between inlet and outlet due to 
condensation or evaporation. Usually, it will be 
accurate enough for engineering design purposes, 
to use an average value of the liquid flow rate, L. 
When the variations in the liquid flow rate are 
relatively large, the value of L can be calculated 
from point to point from the variation in the 
gas humidity. 


EXAMPLE 


The application of the new design method for 
non-saturated gas-vapour mixtures will be illus- 
trated by the following example: An air—water 
vapour mixture enters a tubular cooler condenser 
at 382°F, and 1 atm total pressure. The partial 
water vapour pressure is equal to 0-8 atm. The 
air—water mixture leaves the cooler 


vapour 
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The graphical determination of changes in the gas phase due to simultaneous heat and mass transfer 


condenser saturated at 86°F. Cooling water 
enters at 80°F and leaves at 120°F, and flows in 
countercurrent with the air-water vapour mix- 
ture. For the present purpose, the following 


assumptions are made: 


ho = 200 B.t.u/hr ft? °F 
h, (at gas inlet) 25 B.t.u/hr ft? °F 
A = 978 B.t.u Ib 
Pr Se 1-0 
Specific heat : air, average 
air, 382 


water vapour, average 


= 0-245 B.t.u 1b. °F 
0-25 B.t.u Ib °F 
0-47 B.t.u Ib °F 


water vapour, 382 °F 0-48 B.t.u lb °F 


First, the ratio LG, is calculated as follows : 
Humidity of entering air-water vapour mixture, 


H, (0-8 « 18) (0-2 x 29) = 2-483 lb ‘Ib. 


Enthalpy of entering air-water vapour mixture, 


i 0-245 « 350 2-483 (1071-0 — 0-47 350) 


3153-5 B.t.u/lb of dry air. 


Enthalpy of leaving air water vapour mixture, 
saturated at 86°F, 


i 42.7 B.t.u/lb of dry air. 


Total heat transferred, 3153-5 42-7 
8110-8 B.t.u lb of dry air. 


Ratio L/G, = 3110.8 (40 = 77-77 


Next, the interface conditions are determined 
by the graphical method given recently by the 
author [2]. 


for all details. 


R = hy/h, = 200/25 = 8.0 

C.. (0-8 0-48 « 18) + (0-2 
8-36 B.t.u/lb mole. °F 

,/C,,) (Pr/Sc)** = 18 

2106 

ty = (8-0 x 120 + 382)/9-0 
(Rt, + t,)/(R + 1) 


Reference is made to this publication 


0-25 29) = 


F = (L 978 8-36 


- 149-1 °F 


In Fig. 2, the curved line represents the satura- 
tion line in a semi-logarithmic plot of p, against f. 
Point B, has the co-ordinates tp = 149-1 °F and 
Pp, = 02atm. A straight line with a slope equal 
to (R+1)/2.303 F =9-0 (2-303 x 2106) =0-001856, 
is drawn through point B,. This line intersects the 
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Fic. 2. Graphical Procedure Applied to Exanmpk 


saturation curve in point C,. representing the 
interface conditions at the gas inlet. The co 


ordinates of point C, are read from the diagram, or, 


Dei 0-2455 atm, and f¢, 198-2 F. Thus the 
heat flux at this point is equal to hy (t t.) 
200 78.2 15,640 B.t.u ft®. 


Point A, represents the condition of the non- 
saturated air at the inlet of the cooler condenser. 
The known co-ordinates of point A, are, 
382°F. When a straight 
line is drawn through the points A, and C, the 
A log p,/At 

0-000484. From equation (17) follows, for 
(Pr/Sc)** 1-0, (e* 

10, dlog p,,/dt 
line A, C, may be followed over a relatively short 


Pe 0-2 atm, and t. 
slope of this line will be equal to 


1) 1-0, and Imm dh 
0-000484. Thus, the 


distance, until point A,. 

In this particular case, the two lines r present- 
ing Alogp,/At, and dlog p,./dt, respectively, 
have the same slope in Fig. 2. In general, the 
slopes of the two lines through point A, will be 
different, see Fig. 1. When the slope of the two 
lines is different, it will be obvious that point A, 
is located on the line for d log p,,/dt. 

The co-ordinates of point A, are read from the 
0-2124 atm, and ¢, 832°F. 
0-7876 atm, H 


diagram, or Ps 
Therefore at point Ag, p, 
(0-7876 » 18)/(0-2124 « 29) = 2-30 lb/lb, and the 
enthalpy is i = 0-245 « 300 + 2-80 (1071 +0-47 
300) = 2861-6 B.t.u, lb of dry air. Heat trans- 
ferred in the first increment is then Q, 
3153-5 — 2861-6 = 291-9 B.tu/lb of dry air. 
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Neglecting the heat of the condensate. the tem- 
perature of the cooling water at point A, is 
calculated from L (120 291-9. For 
L Gy 77-77. follows & 116-3 °F, 

For the design of the cooler condenser, the above 
calculation is repeated for various points. When 
necessary, new values for the heat and mass 
transfer coeflicients may be calculated at inter- 
mediate points. After calculation of the heat 
flux at various intermediate points, the condenser 
tube area is obtained by graphical integration of 


the equation, 


d Q 


4 jaa | he (t; t,) 


The graphical integration is carried out by plotting 
the reciprocal heat flux, 1/ hg (t; — t,), against the 
total heat transferred, Q@. The area under the 
integration curve so obtained, is equal to the 
area A. 

The design calculations for direct contact 


equipment are carried out along the same general 


lines. 
NOTATION 
surface urea, ft* 
C, €y specific heat at constant pressure of gas vapour 


mixture and vapour, respectively, B.t.u/Ib °F 
r specific heat of liquid, B.t.u/Ib °F 
. specific heat at constant pressure of gas vapour 
mixture, B.t.u/lb mole °F 
diffusion coefficient, ft® /hr 
“ mass velocity of gas-vapour mixture, Ib /hr ft? 
G flow rate of vapour free gas, Ib /hr 
i! = humidity, lb/lb of vapour-free gas 
h combined heat transfer coefficient for cooling 
water, tube wall, scale, and condensate, 
B.t.u/hr ft® °F, In case of direct contact, 
between gas and liquid, hy reduces to h,, the 
liquid side heat transfer coefficient, B.t.u/hr 


ft? °F 
h.= gas side heat transfer coefficient, B.t.u/hr 
ft? °F 


Jm heat and mass transfer factor, respectively, 
dimensionless 
k = thermal conductivity of gas-vapour mixture, 
B.t.u/hr ft °F 
k gas side mass transfer coefficient, |b /hr ft? atm 
k, /hr ft? 
L = liquid flow rate, lb/hr 
L,,, L; = latent heat of evaporation, respectively at 
main body of gas-vapour mixture, and at gas 
liquid interface, B.t.u/lb mole 
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M..M,, 


Ap 
A In p, 


Pei Por 


Pgs 


- molecular weight of vapour and of gas-vapour Pot 


mixture, respectively, lb /lb mole 

total pressure, atm 

Prandtl number, /k, dimensionless 

partial vapour 
and at the 

respectively, atm 

Po — Pp 

In 

partial pressure of non-diffusing gas component, 


pressure in the gas-—vapour 


mixture, gas liquid interface, 
atm 


ln Per 


respectively, at interface and in main body of 
the gas-vapour mixture, atm 
partial pressure of the saturated non-diffusing 


gas component, atm 
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Shorter Communications 


Zur Berechnung des Volumens adsorbierter Schichten 


Das Volumen ciner adsorbierten Schicht wird meist nach 
F \ berechnet, wobei F die Oberfliche 

\ die Dicke der adsor- 
Forme! gilt 


der Formel A 
des adsorbierenden Koérpers und 
darstellt. 
\ klein gegeniiber dem Kriimmungs- 
Abb. 1 stelle 


beispielsweise der schraflierte Teil den Querschnitt eines 


bierten Schicht Diese jedoch nur, 
solange die Dicke 


radius der betreffenden Oberfliiche ist. In 


Abb.1 Adsorbierte Schicht an einem Korper mit poly 


vonfOrmigem Querschnitt 


stabchenfOrmigen Korpers dar, um den sica eine adsor- 
Wie man sich leicht iiberlegt, 
Querschnitt 


bierte Schicht gelagert hat 
ist der Punkte dargestellte 
Schicht nicht gleich U A (fl 
\ A? 
sprechenden Rechtecken treten niimlich noch die finf 
Kreisflache 


und 


durch dieser 


Linge der Um inyslinie) 
sondern gleich Zu den fiinf den Seiten ent- 
Kreissektoren hinzu, die zusammen eine volk 
L die Lange 


Volumen der an den Enden absorbierten Schicht vernach- 


ergeben. Ist des Stibchens, wird das 


lissigt, so betriigt das am zvylindrischen Teil absorbierte. 
Volumen: 
FA 


L(UA A® 


a 


Wie 
Enden nicht mehr vernachlissigen diirfen, d.h. wenn wir 


steht es nun aber, wenn wir die Wirkung der 


zum allgemeinen dreidimensionalen Problem itibergehen 
Glicklicherweise haben auch diesen Fall die Mathematiker 
Von Sremer wurde nimiich ein 


Kikor 
nach dem 


bereits vorausygedacht. 


Satz *fusseren ParallelkOrper von 
dazu z.B. [2), 


Schicht 


uber den 


pern abgeleitet [1] und vergl. 


sich das Volumen der adsorbierten darstellen 


lisst durch: 


M + (4/3) A 


a 
Als Kikérper wird dabei ein Korper bezeichnet, der durch 
Ebenen oder konvexe Flichen begrenzt ist, also nirgends 
Ecken 
* Integral der 


eingebuchtet ist und nirgends einspringende 
Der Ausdruck M 
mittleren Kriimmung bezeichnet. 
fiir die Kugel (Radius R&), den Zylinder (HOhe und fiir 
das Oktaeder (R Radius der umschriebenen Kugel) 
beispielsweise die Werte 12.56 R #R 
mh, baw. 10,45 . R. Werte fiir zahlreiche andere Kérper 
finden sich in den Tabellen auf S. 36/37 Buch 
Natiirlich nur 
der 


aufweist. wird als das 


Fir M ergeben sich 


im von 


die Formel dann 
Abstand der einzelnen Koérper so 


Schichten 


HADWIGER. lisst sich 


anwenden, wenn 
gross ist, dass sich die an ihnen absorbierten 


nicht tiberdecken 


kidg. Technische Hochschule I’. GRASSMANN 
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kK. A. Guecexurim: Thermodynamics. (3rd Ed.) 
North Holland Publishing Co., 1957. xxii 176 pp. 


UNLIKE the second edition, which was substantially the 
same text as the first, the third edition of GUGGENHEIM’s 
valuable book has recently been published with two 
extensive modifications. 

The more important of the two is a rearrangement of 
the material on solutions and mixtures. In the original 
text, systems of two components were discussed in a 
separate chapter from the chapter in which systems 
containing more than two were considered, and this 
resulted in a good deal of repetition This has been 
avoided in the present edition which contains a new 
chapter on “ mixtures ” where all components are treated 
on the same basis 

As a result of this change the number of chapters in 
the book would have been less by one, but for the second 
major alteration, which is an entirely new chapter on 
“ Onsager’s Reciprocal Relations. Here the author 
gives an outline of what has become known as the thermo- 
dynamics of irreversible processes, or the thermodynamics 
of the steady state. This chapter occupies only 10 pages, 
with a result that the applications described are rather 
limited and do not include non-isothermal systems. 

Two other new features of the book are an “Introduction 
on Notation” and the use of Struue’s values (1955) of the 
physical constants in place of those of Bircre (1941) 
The most interesting change here is the value of the ice 


wint which is quoted 273-15 K in place of 273-16 K. 


K. G. Denpicu 


H. Hausen: Handbuch der KAaltetechnik. Vol. 5, 
The Liquefaction and Separation of Gas Mixtures 
(In German). 


Tue industrial application of low temperatures has 
developed during the past 50 years; and particularly 
during the past 20 years the development has been very 
rapid. Throughout most of this period Hausen has been 
in intimate contact with the subject, to which he has 
made many notable contributions. Thus the new textbook 
devoted to the low temperature industry, and forming 
part of the massive series on refrigeration edited by 
Rudolf Plank, speaks with authority and from rich 
experience, 

The book is divided into two parts, the first and larger 
part dealing with the theory of gas liquefaction and 
separation, the second part chiefly with technical aspects 
of low temperature processes. 

The first 50 pages are devoted to a review of the laws 
of thermodynamics and the properties of real gases. Equal 
space is then given to the theory of gas liquefaction ; 


1838 


all the common liquefaction cycles are discussed, including 
the new Philips gas expansion machine. The thermodyna- 
mics of gas separation and then treated, firstly by outlining 
the laws governing the behaviour of mixtures, and 
secondly by describing the properties of a rectification 
column and its associated design methods. 

Air separation is dealt with at length, all the common 
industrial cycles being reviewed. Tonnage oxygen plant 
features largely but not to the exclusion of plant for the 
production of high purity gases and the recovery of the 
rare gases. The separation of all other industrial gas 
mixtures is dismissed in 15 pages. The first section of the 
book concludes with the theory of heat exchangers and 
regenerators, subjects in which the author has done much 
pioneering work. 

The second part of the book includes very useful sections 
dealing with the practical design or selection of equipment 
for all the major components of a low temperature plant. 
Much useful detail is given in a field which has not previously 
been treated to any extent in a textbook. Some cost 
data are also provided. A short section is devoted to the 
physical properties of fluids at low temperatures, and the 
book concludes with a review of the applications of the 
products of low temperature plant. 

Two criticisms may be made of the work. On the one 
hand, in attempting to be completely comprehensive the 
author has included material, particularly at the beginning 
and end of the book, which could without great loss have 
been omitted, whilst treating some of the centre sections 
too briefly. On the other hand the author is biased 
strongly in the selection of material by his close association 
with the Linde Company, and this has resulted in a dispro- 
portionate concentration on air separation as compared 
to other technical gas mixtures. However, in view of the 
amount of useful and well ordered information contained 
within this volume it is likely to become a standard work 


of reference in low temperature technology. 


G. G. HAsELDEN 


Chemical Engineering Practice. Vol. 3: Solid 
Systems. Edited by H. W. Cremer and Treror Davies. 
Butterworths Scientific Publications, London, 1957. 
vi 534 pp., 95s. 

Tuts journal is concerned exclusively with the application 
of the sciences to chemical engineering problems. No 
doubt most of its readers, when not actually absorbed 
in its pages, have some contact with the rough-and- 
tumble side of the profession, in which it is more important 
to make something work than to know why it does so 
Nevertheless, we may suspect that they tend to be drawn 
from the “ Egghead” wing of the profession, and as 
such they are likely, paradoxically, to find the volume 
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Book reviews 


under review more stimulating than its two predecessors. 
Volume 3 really gets down to brass tacks ; here we are 
in the world of practical engineers at last. The volume 


is full of machines —for crushing, grinding, classifying, 
which are the product of experience, 


Science 


mixing and so on 
intuition and remarkable ingenuity. seems to 
have had very little to do with the development of most 
of them ; field 
even caught up in judging from the 
attached to many of the 


in fact, far from having led the it has not 


most cases, very 
meagre shreds of “ theory ”’ 
topics. It is true that the theoretical devek pment has 
gone further than might be supposed from reading some 
of the sections — that on cevclones, for instance. In general, 
however, the techniques described pose a large number of 
unsolved and interesting problems, which should serve 
to stimulate the theoretician. In many cases the general 
principles involved are clear, but there is no quantitative 
theory which can be used for that accurate prediction of 
For 


performance which forms the basis of good design. 


instance, although a number of tentative ipproa hes 


have been made to the theory of mixtures and mixing, 
they have probably had no influence at all on the methods 
of specifying the design and size of a mixer for a given 
job — methods which are about as scientific as winetasting. 
Dr. H. Heywood, in his interesting article on the Pring iples 
idea of the 


which a scientific attack can be opened on an intractable 


of Crushing and Grinding, gives an wav in 


problem. 


In general, however, the theoretical treatment is 


perfunctory, and the book consists mainly of descriptions 


of machines, their applications and mode of action, and 
illustrated, 


which are (for the most 


is copiously with diagrammatic drawings 


part) simple enough ty 
The 


crushing, grinding and pulverising ; 


‘ \pose 


the mode of action clearly. topics covered are 
size analysis, s¢ rechiny, 
sedimentation ; the 
and the 


I believe the volume will be a useful preliminary 


grading, classifving and mixing 


storage and handling of solids: cleaning of 
vases. 
reference book for those in industry who are unfamiliar 
with these fields, and who want to know the general types 
of equipment which are available. However this may be, 
wholesome 
more academic chemical engineer Apart 
fact that 


inventiveness is at least as important as scientific analysis 


it cannot be too highly recommended as 
reading for the 
from its intrinsic interest, it drives home the 


in chemical engineering. P. V. Dancxwerts 


Modern Chemistry for the Engineer and Scientist. 
Edited by G. Ross Roperrson. 


different 
of high 


book 
The chapters are, without exception 
but the 
style and objective are marked. The title 
of the book suggest that the text might form a suitable 
for the 
to modern chemistry. This is only true in the sense that 


Tuts consists of 19 chapters cach by 


authors. 
interest, differences in 


technical and scientific 


ind preface 


introduction non-chemical scientist or engineer 


such a reader would discover some of the topics which 
make up the subject but this could be achieved in perhaps 
one quarter of the present text, much of which is taken 
In this 
respect the book is probably better suited to the applied 


up with specific examples rather than concepts. 


chemist, chemical engineer or even the specialist chemist. 

It is not possible to classify the chapters unequivocally 
in terms of objectives, but those on thermodynamics, 
contact catalysis, photochemistry, rubber elasticity, metal 
mechanisms 
The 
last of these is exceptional in demanding specialised 


creep, clay minerals and organic reaction 


are largely devoted to reviews of fundamentals. 


knowledge and vo« abulary in the reader. All are accom- 
panied by good references and bibliographies. 

\ second group of chapters is devoted to measurement 
techniques including isotopic tracer methods, column 
chromatography, and methods for the study of fast 
reaction kinetics. These chapters are uniformly interesting 
and well presented although that on chromatography 
relates almost entirely to specialised branches of organic 
chemistry. 

\ third group of chapters provides reviews of the 
production and properties of silicones and carbon fluorine 
compounds, the petroleum chemical and food technology 
industries and of chemical relationships with the earth. 
These reviews accord well with the professed objectives 
of the book each giving a well rounded and informative 
account of its topic. That on petroleum chemicals is 
exceptionally good : one regrets the paucity of its list of 
references. 

The last four chapters of the book are devoted to 
biochemical and biomolecular topics, the biochemistry of 
insecticides, chemistry and physiological action in the 
living organism, and the structural chemistry of proteins. 
They will tax the scientific vocabulary of the average 
chemical engineer without direct reward, although the 
survey of synthesis in the living organism gives a most 
useful bac kground and a proper sense of modesty. 

In summary the text is worth-while reading and 
probably a worth-while purchase for the chemical engineer 


with a chemical bent. Most of the chapters repay several 


readings. N. L. 


FRANKLIN 


and T. Pierorp : Nuclear Chemical 


McGraw-Hill, London, 1957. xiv 597 


M. Benepicr 
Engineering. 
pp.. Tis. 6d. 


Tuts book is intended for the student although it includes 


a number of graphs and tables and comprehensive 


bibliography which are a useful source of reference. 
It presents with extreme clarity an account of some of 
the process aspects of nuclear engineering, (¢ hapters on 
isotope separation and one on properties of irradiated fuel 
is particularly notable. 

After 


reactions, there is a detailed account of fuel eveles and 


a chapter of introduction and one on nuclear 


irradiation schemes. The algebraic manipulation of this 
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Book 


section conceals the practical uncertainties in effective 
Much 
of the discussion is based on reactivity arguments, their 


neutron cross-sections and reactor performance. 


relevance to chemical engineering is not clear. Subsequent 
chapters on uranium production and on special nuclear 
materials are impaired by the absence of information 
on fabrication techniques, fuel element fabrication being 
dismissed in nine lines. This makes it impossible to define 


adequately the properties, other than nuclear purity, 
which are being sought in the production techniques. 
An adequate chapter on solvent extraction is followed 
by an excellent discussion of irradiated fuel properties 
which defines the requirements of irradiated fuel processing 
indication of the theoretical basis 
The 


discussion on irradiated fuel processing is disappointing 


and which gives an 


for solvent extraction processing. subsequent 
to the specialist. There is no discussion of effluent disposal 
problems, and plant design concepts, plant control, process 
equipment, plant inspection and decontamination escape 
comment. Even criticality control in enriched uranium 
or plutonium receives only passing mention The ch ipter 
provides a simple and clear discussion of the flowsheets 
of a number of solvent processes but excludes most of the 
detailed 1955 


Geneva Conference 


more information declassified since the 


The last third of the book Is devoted to isotope st par i 


tion. The treatment is related periodically to the more 


familiar concepts of distillation theory, an advantage to 
the chemical Classification 


engineer. precludes any 


discussion of gaseous diffusion equipment technology, 
treated 


principles 


and transient behaviour of cascades is not 


within these limitations the treatment of the 
of isotope separation is excellent. 


reviews 


The 
specialist alike. 


book can be recommended to the student and 
The latter will regret that classification 
and timing have prevented the authors from treating a 
number of topics which would have benefited from the 
clear analysis and style which they apply in this book. 
N. L. 


FRANKLIN 


Kdited by K. 
200 pp., 80s 


Chemical 


Reaction Engineering. 


Pergamon Press, London, 1957. 


Tuts book is a collection of fourteen papers, mostly in 
English, presented at an international symposium held in 
Holland in 1957. 


Chemical reaction engineering, which has as its ultimate 


aim the successful design of commercial reactors from 


first principles, is still in its infancy. Important in this 


field is not only the study of fluid mechanics, heat and 


mass transfer, and chemical reaction kinetics, but also 


the interrelation of these subjects. The problems involved 
little 


are complex and as yet progress has been 


towards their solution. 


The symposium has been divided into five parts 


introductory papers, transport phenomena in _hetero- 


geneous reactions, non-uniform concentration distributions 
reactor efliciency and stability, and reactor development 
Of necessity, much of the material is in the nature of a 
review, but some interesting new work is also presented. 


This volume will undoubtedly prove valuable to all 


interested in this field. 


G. A. Rarer: 
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